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ABSTRACT 


An aggregate production planning system has been 
designed and developed for use in a variety of small manu- 
facturing industries. The central analytical component 
of this system is formulated within the framework of the 
Search Decision Rule for solving the operations planning 
problem. The production/cost models used in the system 
are designed to reflect the typical operating environment 
encountered in most manufacturing industries. The system 
performance was assessed using a hypothetical planning 
problem. Further empirical testing is recommended in 


order to assess the practical capabilities of the system. 
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CHAPTER 1 


INTRODUCTION 


1.1 Purpose 


This research is directed at developing an aggregate 
production planning system as part of a total corporate 
planning model that can be used by small manufacturers. 

It is hoped that this research will offer the small manu- 
facturer a planning model that can improve his competitive 
position in todays business world. The production planning 
model developed herein uses a search routine ER treetacts 
with the objective function to generate near optimum values 
for each period (month) for (1) direct work force levels, 
(2) indirect work force levels and (3) production volumes 
by product line. With values for these three parameters 
reasonably complete cost data for the production planning 


purposes is generated and output by the system. 


1.2 Background Information 


Considerable work has been done with respect to the 
development of aggregate planning models. These models 
fall into the general classification of the Linear Decision 
Rule (LDR), linear programming methods, heuristic programn- 
ing models and computer search methods in the form of the 


Search Decision Rule (SDR). 
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The development of the Linear Decision Rule by Holt, 
Modigliani, Muth and Simon [20,21,22] in 1955 represented 
a major step in the development of a planning model that 
appeared to have some real value in a practical sense. The 
model was tested in a paint factory using actual past data 
over a period of six years. The analyses indicated a mini- 
mum cost saving of 8.5 percent by the use of the LDR over 
actual management performance. The LDR model formulation 
involves the development of a quadratic cost function for 
the company in question. The cost function consists of such 
costs as the regular payroll, hiring, layoff, overtime, 
inventory holding, backordering and machine setup costs. 
Then, using this quadratic cost function two linear decision 
rules are derived for computing work force levels and pro- 
duction rate for the upcoming period based on the forecasts 
of aggregate sales over some planning horizon. These two 
linear decision rules are optimum for the model. The LDR 
was later abtenaty to a multi-item production situation by 
Bergstrom and Smith [1] who refer to it as the "Multi-Item 
Decision Rule (MDR)". The Linear Decision Rule has been 
used aS a comparative test model By hand other individuals 
over the past twenty years. 

Among the linear programming models work by Bowman 
[3] and by McGarrah [29, pp. 124-127] are worthy of mention. 
Bowman proposed the distribution model of linear programming 
“as a mode for aggregate production planning in 1956. The 


model was built around the objective of allocating units of 
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productive capacity such that the total production and 
storage costs were minimized and sales demand was met with- 
in the constraints of available capacity. McGarrah pro- 
posed a simplex model for aggregate production planning 
covering such aspects as production change costs, inventory 
holding and backordering costs. The major disadvantage 
with these models is that they are based on a planning hori- 
zon of only one period. The simplex models which have 
larger horizon times have also been developed by McGarrah 
[29, pp. 127-129] and Hanssmann and Hess [18]. The Hanssman- 
Hess model determines the work force level and the production 
rate as independent variables with the regular payroll, 
hiring, layoff, overtime, inventory, and backordering costs 
as dependent variables. The model allows for a presetting 
of the planning horizon time. 

The linear programming methods and the Linear Decision 
Rule produce mathematically optimum decisions provided the 
basic cost structure of the model accurately represents the 
dynamic operating characteristics of the firm. These models, 
in fact, are unable to realistically duplicate the actual 
operating conditions because of their very basic require- 
ments of continuities and specific forms for the cost func- 
tions. 

Under the heuristic programming methods, there are 
two rather different proposals for handling the aggregate 
production planning problems. The first, Bowman's manage- 


ment coefficients model [4], uses past average managerial 
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performance to derive coefficients for decision rules. 
Bowman hypothesizes that management's decision-making 
criteria is in no way different from that used in analytical 
models except that managerial behaviour tends to be highly 
variable about optimality rather than being off-centre. 
Thus, he presumes that overall managerial decisions can be 
improved considerably over a period of time through consis- 
tent use of decision rules. The industrial application: of 
this technique is limited due to the fact that the manager 
may not be fully aware of (the components, the exact nature, 
and the dynamic characteristics) the firm's actual cost 
structure and the lack of this knowledge will lead to the 
biasing of decisions. The second heuristic model is called 
Parametric Production Planning (PPP) and was developed by 
Jones [24]. The PPP model is based on a simulation approach 
built around the postulated existence of two specific linear 
feedback rules with coefficient values obtained from a 
four dimensional search of possible parameter values derived 
from a linear and/or non-linear model of the firm's cost 
structure. This model does not place any limitations on 
the mathematical form Of¥thescosttiuicti ons, paspihenemane 
with the Linear Decision Rule and linear programming. In 
the paint factory application the PPP model yielded decisions 
which were almost the same as those given by the optimum 
Linear Decision Rule. 

Finally, the Search Decision Rule (SDR), developed by 


Taubert [42,43], utilizes a modified Hooke-Jeeves computer 
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search method. Under the Search Decision Rule, the plann- 
ing problem is first defined as a multi-stage decision 
system and then solved for a global optimum by using a 
computer search technique as a heuristic optimum-seeking 
method. The search method used in SDR, the "Pattern Search", 
was developed by Hooke and Jeeves [23] and later modified 

by Weisman, Wood, and Rivlin [50], and by Taubert [43]. 

This method does not require the cost/criterion functions 

to be either linear or quadratic but can take on any form 
desirable to duplicate reality. The SDR was tested by 
Taubert in a variety of planning situations. In the paint 
factory application the SDR decisions were only 0.1 percent 
above the analytically determined optimum computed by the 
LDR. Other successful applications of SDR include an exist- 
ing research company, where SDR during a test with actual 
past data over a 5.5 year period produced highly realistic 
decisions. The test results indicated that the SDR decisions 
would have produced a potential 11.9 percent cost saving 
over the actual management decisions. In a recent simula- 
tion test conducted by Lee and Khumawala [26] using both 
perfect ad imperfect forecast data SDR has outperformed 
Bowman's Management Coefficients Model, the Linear Decision 
Rule and Jones' Parametric Production Planning Model. Thus, 
SDR seems to have significant potential for the solution 

of more realistic planning models. Therefore, the SDR 
approach has been selected for use in this research project 


for developing an aggregate production planning system for 
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small manufacturing industries. 


testeScope 


The aggregate production planning system designed 
herein is structured around the Search Decision Rule. The 
computational model for the system is designed to handle 
deterministic data only. The system programs consist of: 
(1) a main routine which contains the master control 
strategy; (2) a planning cost subroutine for computing the 
costs. and the criterion functions: and (3) a search routine, 
the “subroutine PATS" taken from reference [40]. The pro- 
grams are designed to work out near optimum values of the 
independent decision variables (direct and indirect work 
force levels and production volumes by product line) for 
each period. The system programs have built-in routines to 
compute various planning details based on the optimum values 
of the independent decision variables. These details 
include the amount of overtime, the level of plant capacity 
utilization, and the period/month ending inventory and 
shortages levels for each product line. The planning cost 
subroutine is designed to compute various operating costs 
such as regular payroll, overtime, hiring and training, 
employee termination, shift premiums, inventory holding, 
backordering, and production fluctuation costs. 

The system has been tested using a hypothetical pro- 
duction planning situation representing a realistic planning 


environment. The system performance was assessed on the 
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basis of this preliminary testing and conclusions were 
drawn and listed in the report. The next logical step is 
to test the system empirically in actual factory planning 
Situations in order to make desirable refinements. Sugges- 
tions have been made which should be helpful in refining 


the model for use in a practical environment. 
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CHAPTER 2 


THE AGGREGATE PRODUCTION PLANNING SYSTEM 


2.1 System Overview 


The aggregate production planning system, presented 
herein, is a time-shared type computer-based system. The 
System generates all pertinent planning and cost information 
necessary for scheduling a suitable production plan to meet 
the sales projections. Ideally, for the most effective 
operation, the system should be administered and made avail- 
able for use to small manufacturers through a service 
organization. 

The entire system is sub-divided into functional sub- 
systems in order to proceed with the development and the 
ultimate design of the total system in an organized manner. 
These sub-systems are: 

A. Input information management system. 

B. Optimization system. 

C. Planning data generating system. 


Each sub-system is described in detail below. 


A. Input Information Management System 


The input information management system controls all 
the input information flows into and within the aggregate 


production planning system. The input information consists 
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of management assumptions about the company's goals and 

objectives including such information as the sales forecast 
data, operating constraints, production cost data and plann- 
ing horizons. The system function begins by processing the 


raw data to prepare a workable data base. 


B. Optimization System 


The optimization system is designed within the basic 
framework of the SDR strategy. Under the SDR approach the 
problem is defined as a multi-stage decision model and 
solved for a global optimum through use of a computer search 
technique by optimizing each stage over multi-period plann- 
ing horizons. Where, a stage refers to any particular point 
in time when the decisions are to be made, it is usually, 
but not necessarily, monthly. 

The optimization system consists of: (1) a planning 
cost model which defines the cost structure of the firm; 

(2) a set of constraints to define the boundaries of the 
solution space; and (3) an optimization technique (the 
computer search routine). The system program contains, 

(1) a master control routine (included in the main routine) 
to operate the system according to the SDR stragegy; (2) a 
planning cost subroutine to compute various costs, planning 
details and the value of the objective function during each 
iteration; and (3) a search routine, "the subroutine PATS" 
taken from reference [40]. The subroutine PATS seeks opti- 


mum values for the independent variables (production levels 
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for each product line and the sizes of direct and indirect 
work forces) for each period in the entire planning span 
through an iterative process directed by the master control 
program. It does so by minimizing the value of the objective 
function which is the sum total of the variable production 
costs over the planning horizon and the penalty charges for 
violating one or more of the operating constraints. The 
variable production costs include such costs as the work 
force payroll, overtime, employee hiring and training and 
termination, shift premiums, inventory carrying, backorder- 
ing and the production change costs. The optimization of 
each stage's decisions is accomplished through continuous 
interaction between the search routine and the planning cost 
subroutine. The entire optimization process and the com- 
plete search procedure are described in detail in Section 


4.3. 


C. Planning Data Generating System 


The planning data generating system receives optimal 
values of the independent decision variables from the 
optimization system. The system routines then compute and 
output, for management, all the necessary planning details, 
such as the amount of overtime, level of the plant's capa- 
city utilization, inventory and shortages levels, work 
force, inventory holding, backordering and other operating 
costs olisitedsingSection 243): 


The results of each run are reviewed by management 
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and if desired, the system is rerun with an alternate set 
of assumptions. This process is repeated until a solution 
acceptable to management is obtained. 

The structure of the computational algorithm designed 
for the system is; 1)1ustrated inPFiguréesé2d. te This figure 
also depicts how the computational algorithm fits into the 


framework of the sub-systems. 


2.2 Input Information Requirements and Data Formats 


The entire input information required to initiate 
a run.on the aggregate production planning model has been 
classified into different data set categories from the 
standpoint of programming. Included in this section is a 
complete description of the total anime Be edema ti on require- 
ments and the details of the input data format specifica- 
tions (i.e., how the input data is to be punched on the 
cards). 

Before discussing the input information requirements 
some of the terms used in this section are defined. 

Stage. In the aggregate planning system, using the 
search decision rule the entire planning problem is defined 
in terms of a multi-stage decision model, where each stage 
represents a point in time (the beginning of each successive 
month/period) when the decisions are to be made concerning 
the operations of the firm. 

Planning Horizon. The planning horizon is the time 


span over which the operation of the system (firm) is to 
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Figure 2.1: Flow Chart of the SDR - Computational 
Model for the Aggregate Production 
Planning System 
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be optimized during each planning stage. 


ence | "LIiStvotothev input Data 


In the list given below, each input data-variable 
is described along with the corresponding notation used 
to represent it, in the same sequence (under each respec- 
tive data set class) as it has to be punched on the input 


data cards. 


Rae oMaster Controlavariabpies 

(1) NOY the length of the entire planning span 
in years. This variable must be a whole 
number. Therefore, if its value is in 
between two integers it should be rounded 
off to the next higher whole (integer) 
number. 

(2) NPY the number of periods per year, must be 
a whole number. 

(3) NH the number of periods in each planning 
horizon, must be a whole number. 

(4) NHM - the number of periods/months which would 
serve as an upper limit on the length 
of the planning horizon. It is used for 
convenience in handling the input data. 

(5) JJ the total number of product lines or 
groups being manufactured by the firm. 


(6) NHP the number of regular hours of work in 
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JPROPX 
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each Period/month. 
corresponds to a technique developed 

DY Tk. Sikes C40] in which an arbitrary 
Factor (JPROPX Zl) 4lscusedmtoamaditee ip 
the month/period ending inventory level 
close to the desired minimum level. 

It does so by penalizing the inventory 
deviations from the specified level for 
the periods close to the end Of the 
Planning horizon. The higher the value 
Of this variable the more effective act 
TSeeelits ava Vue can be set anywhere between 
T-and'50 (it ds quite effective for 
values between 19 ANG 15) tee JPROPX. can 
be made totally ineffective by setting 
Up ee coal Be 

a controlling Variable, whose Value jis 
either one or zero, branches the choice 
Of approach for estimating the indirect 


work force requirements in the model. 


Regression method is used when the value 


is one, Extrapolation is used when the 
value is zero. 

an upper limit on the maximum number of 
Shitts, the Plant can be Operated, 

a print-out control Character equal to 


one or zero, (suppresses a large fraction 
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(99 percent) of the detailed output of 
the analysis when the value is set equal 
to one and allows the full details to 

be printed out when its value is equal 
LOEZEGOt 

an approximate upper limit on the total 
number of search iterations (for optimi- 
zation) during each stage. Its value is 
usually set at 3000 (see Appendix I) for 


more details. 


Data Set -"Typeni 


This set contains information regarding the initial 


conditions at the start of the first planning stage. 


(1) 


(2) 


(5) 


OWD 


OW! 


OED 


OE} 


OTDPI 


the beginning direct work force level at 
the start of the month/period one. 

the beginning indirect work force level 
at the start of the month/period one. 
the average employee efficiency level of 


the direct work force at the beginning 


'of the period one. 


the average employee efficiency level of 
the indirect work force at the start of 
the period one. 

the level of the direct work force over- 
time, expressed as a fraction of the 


work force level, used during the period 
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prior to the start of the plan. 


(6) OTIPI same as above, for the indirect work 
force. 
bh) MSO the number of shifts, the plant was ope- 


rating in prior to period/month one. 


©. Data Set = ly per 


This set consists of all the information which is 
different for each product line and, hence, will be entered 


into the data file as separate sub-sets (one per product 


line). 

(1) OP the production level during the period/ 
month prior to period/month one. 

(2) OVI the beginning inventory level at the 
start of month/period one. 

(33) S the sales forecast for a number of periods/ 
month = NOY x NPY + NHM - 1 (see the 
master control data set for the defini- 
tions of these variables). 

(4) VIL the lower limits on the levels of the 
month/period ending inventory stocks for 
each year of the plan. 

5) ENO the upper limits on the levels of month/ 
period ending inventory stocks, for each 
year in the plan. 

(6jey cCR the cost of carrying inventories per 


dollar value of the total inventories. 
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ICR is computed by summing up the follow- 

ing rates: 

RMAR = the management's desired 
after-tax rate of return on 
its investments. 

ISR = the cost of inventory losses 
through breakage and pilfer- 
age per dollar value of the 
total inventories. 


ISOR 


the cost of inventory losses 

through obsolescence and 

decomposition of the product 
per dollar value of the total 
inventories. 

IITR = the insurance and tax costs 
per dollar value of the total 
inventories. 

IOFR = the overhead costs of carry- 
ing inventories per dollar 
value of the total inventories. 

the desired upper limit on the level of 

inventory shortages (expressed as a frac- 

GON OT | thie: Ss anics:) . 

the cost (per unit) of processing the 

customer requests for better delivery 

periods for the sales orders postponed 


because of a stock shortage. 
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where 

F = the (expected) total annual 


costs of processing these 
requests’. 
M = the average number of orders 
(per year) met by shortages. 
AOS = the average order size, number 
of units per order. 
the cost of goodwill (per unit) that will 
be host.1f a distributorucerminares his 


business with the manufacturer. 


CGW Se yas) 

where, 

g = the management's desired after- 
tax rate of return on the costs. 

ASC = the average dollar value of 
sales per distributor measured 
in terms of the costs. 

AOS = the average size of a sales 
order (i.e., the number of 
units per order). 

(p/a) = the present value factor 

= ; {] - Coane 


sale ve nee 
a — oy 


zvebyo to 1sdaun BpRTSYS One 


Tree eee ee ) ‘et Seal 
-2e0styore vd tem (rae eq 
a i 
os er awekh«. wdehtan VAD es = . 2 
Saini U i » | “4 a j ’» SOS 7SV5 oat : 6 
te 7 +3. na 3 
‘ Do Voq e@eyrnu 79 - ; 
: eit “a 
a a P _ , r - . re jhe q + ~ ‘ ’ " : 
1ifw Fens (LIT Wap i pirweRoey Te Jeo2 ond . i 


bnitmyes Vosudtaseth s TF taps 3d 


sud >6TURSM) ed Hitw ee nfevd 


: a 7 , 
(o\q).2. ee) ne ae 
' CUA 7, F . 
; ; ef. 
pt 
.ototw ‘ 
- 
ssifasb 2‘insmepansm say = p 


y i 


: 
S20. snt no axnuteyt te stey xs2 ; ri 


4 
Qo sufev wsifeb epsreya ett = 224 
| a a rh 
bsvuzcam yosludirteth veaq: gefa2 Pes 


.2e2e00 aft Yo me ead 


“oM - ™ & 

eefé2 6 to exte s e1OVE ods Serne . A 
i. Pt. Be ~) 

7 iy al 

on 
(yeb49 194 ettny, an ‘) 


rotost oulsy ingore 5 \ 


to todmuna on? ney abv 


(10) 


(alee) 


(12) 


(13) 


TPTU 


16 


where 

i = the management's desired 
after-tax rate of return on 
its investments. 

the probability that the requests for 

improvements in delivery periods will be 

received if the sales deliveries were 

postponed because of the stock shortages. 

the probability that an order will be 

cancelled if it is met with a condition 

of stockout. 

the probability that the distributor will 

terminate business with the manufacturer 

(for this particular item only) if the 

deliveries are postponed because of stock 

shortages. 

the average production time (man-hours) 

per unit of the product line. This should 


be determined as: 


k 
TPTU = _—+— Tee ital UA Eig mex ori) 
pr sudiiew 
i=1 
where 
j = 1,2,3,...,k, the i*" product- 


jtem in the product line. 


SJ(i) = the average yearly sales level 
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for the product-item (i). 
TPTU(i) = the production time (man-hours) 
per unit for the product-item 
(i) (its value is equal to the 
inverse of the production 
rate). 
r4)o AL the average volume of the inspection and 
quality control work-load, man-hours per 
unit of the product line. Its value may 
be set = 0 if NWI = 1 in the master 
control data. 
(15) CDMU the average cost of direct materials per 
unit of the product line. For a more com- 


plete definition refer to Section 3.9. 


Doe wDatagset - Typess 


This set includes the information which is required 


separately for each individual year in the entire planning 


period. 

Gly AED the total yearly depreciation costs. 

(2) ARM the annual expenses for repairs and 
maintenance of buildings and office 
equipment. 

(3) ARN the yearly rental costs. 

(4) AUT the utility expenses for each year. 

(5) AOS the yearly expenses of the office supplies. 


(6) APT the yearly insurance and property tax 


(7) 


(11) 


ees) 


(13) 


(14) 


(15) 


ADD 


ASW 


ASO 


AIS 


COER 


CPER 


MDHU1 


MDHU2 


MDHL1 
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expenses, 

the annual engineering and drafting 
expenses. 

the yearly costs of scrap, rework and 
in-house material losses. 

the yearly costs of outside services such 
as contract services and travel expenses. 
the yearly indirect overhead expenses such 
as personnel support and training depart- 
ment expenses. 

the rate for the costs of operation, 
repairs and maintenance of the plant 
equipment and machines, dollars per pro- 
duction man-hours. 

the rate for the power and electricity 
costs, dollars per production man-hour. 
the maximum direct work force requirements 
(man-hours per period), for each year, to 
operate the plant during the first shift 
without any overtime 

the maximum direct work force requirements 
(man-hours per period), for each year, to 
operate the plant during the second or 
each subsequent shift without any overtime. 
the minimum direct work force requirements 
(man-hours per period), for each year, to 


operate the plant during the first shift 
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without overtime. 

(16) MDHL2 the minimum direct work force requirement 
(man-hours per period), for each year, to 
operate the plant for the second or each 
subsequent shift without overtime. 

G7 jie WI the maximum indirect work force require- 

| ments (number of employees), for each 
year, to operate the plant during the 
first shift without overtime. In case it 
Lsefdi FFiCUhG ‘Vol estimate its: value, ‘seit 
it = 0 and the model will extrapolate 


the value internally. 


Pevauaeoet: = ype met 


This data set consists of the miscellaneous type of 
information. 

The variables ‘umnbiers "WW2))y 72°38) 9 and 10° (listed 
below) become inoperative when NWI = 1 in the master 
control data. Hence, their values may be set = 0 when 
NWI = 1. Similarly, variables numbers 4, 5 and 6 may be 
set = 0 when NWI = O. | 
(1) XI the minimum indirect work force require- 

. ments (the number of employees) for operat- 
ing the plant during the first shift with- 
out overtime. 

(2) X13 the minimum number of indirect work force 


employees required to maintain the plant 
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(operations) during the first shift for 
a zero direct work force level. — 

(3) X14 the minimum number of indirect work force 
employees required during the second or 
each subsequent shift to maintain the 
plant (operations) for a zero direct 


work force level during that shift. 


(4) X the regression coefficient, see Appendix 3. 
(5) Y the regression coefficient, see Appendix 3. 
(6) Z the regression coefficient, see Appendix 3. 
jm was the number of supervisors required per 


employee of the direct work force. 

(8) AM the average volume of plant repairs and 
maintenance work load (man-hours) per 
production man-hour. 

(9) AP the volume of planning, scheduling, rout- 
ing and inventory and stores work in man- 
hours per production man-hour. 

ClO)” “AC the number of truckers, shipping and time 
clerks required per employee of the 
direct work force. 

(11) PP the ratio of the number of production 
planning and scheduling employees to. the 
total number of employees of the indirect 
work force. 

G12) .PD4 the average proportion of hirings in the 


direct work force from skilled or experi- 
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(18) 
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EHD] 
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NHTD1 
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enced type of employees, referred to as 
class (a). 

the average proportion of the direct 
work force employees hired from unskilled 
or semi-skilled type of employees, 
referred to as class (b). 

the average initial efficiency, at the 
start of the training period, of the 
direct work force employees from class 
(a). Assuming 1.0 as the efficiency level 
for an average trained employee. 

the average initial efficiency, at the 
beginning of the training period, of the 
direct work force employees hired from 
class (b). Assuming an efficiency level 
equal to 1.0 for an average trained 
employee. 

the average duration of the training/ 
learning period in hours for the direct 
work force employees from class (a). 

the average duration of the training/ 
learning period in hours for the direct 
work force employees from class (b). 

the average proportion of each hour the 
supervisor spends in directly training 
each new employee of the direct work 


force from class (a). 
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(20) 


(23) 


(24) 


(25) 


(26) 


(27) 


STD2 


RD 


CRSD 


CTFHD 


PHDU 


eatty 


PI2 


NHTI1 


NHTI2 


ee 


the average proportion of each hour the 
supervisor spends in directly training 
each new employee of the direct work 

force from class (b). 

the effective learning rate constant (for 
the combined plant operations) for the 
direct work force employees. A method 

for estimating this learning rate constant 
is described in the Appendix 4. 

the average recruitment and selection costs 
for hiring a direct work force employee. 
the fixed costs of training for the direct 
work force per trainee, per hour. 

the upper limit on the rate of hiring in 
the direct work force in any period/month. 
the average proportion of hiring in the 
indirect work force from skilled/experi- 
enced employees, referred to hereafter 

as class (a). 

the average proportion of the employees 
hired in the indirect work force from 
unskilled or sémneakeaied employees, to 

be referred as class (b). 

the average duration of the training/ 
learning period in hours for the indirect 
work force employees from class (a). 


the average duration of the training/ 
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learning period in hours for the indirect 
work force employees from class (b). 

the average initial efficiency, at the 
start of the training period, of the 
employees hired in the indirect work 

force from class (a). Assuming an effici- 
ency level equal to 1.0 for an average 
trained employee. 

the average initial efficiency, at the 
beginning of the training period, of the 
indirect work force employees from class 
(b). Assuming 1.0 as the efficiency level 
for an average trained employee. 
the (average) proportion of each hour the 
supervisor spends in directly training 
each new indirect work force employee 

from class (a). 

the (average) proportion of each hour the 
supervisor spends in directly training 
each new indirect work force employee 


fvomecdasse lb) 4 


the recruitment and selection costs for 


hiring a new employee in the indirect 
work force. 

the fixed costs of training for the in- 
direct work force employees per trainee, 


per hour. 
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(34) PHIU the upper limit on the rate of hiring 
in the indirect work force in any period/ 
month. 

(35) RI the average (effective) learning rate 
constant for the indirect work force 
employees. Its value can be determined 
in the same way as explained earlier for 
variable number 20 (RD), or it may be 
set at approximately 0.75 to 1.00. 

(36) PTDU the upper limit on the rate of termina- 
tions of the direct workforce employees 


in any period. 


(37) 2 (PITRE the maximum acceptable rate of employee 
turnover. 
(38) see bli the upper limit on the rate of termina- 


tions of the indirect work force employees 
in any period/month. 

(39) ES the expected value by which the average 
employee efficiency level would drop at 
the beginning of the period of a new shift 
program (starting or shutting down an 
additional shift). The efficiency falls 
because some of the work force has to be 
redistributed in the process of the pro- 
gram change. The estimate is to be based 
upon an efficiency level of 1.0 for an 


average trained employee. 
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the expected value by which an employee's 
efficiency level would drop after some of 
the operations have been reorganized for 
changing some of the production output 
levels. The estimate is to be based upon 
an efficiency level of 1.0 for an average 
trained employee. 

the time in hours that would be required 
by an average employee to recover the 
efficiency lost in the process of a shift 
program change. 

the time in hours that would be required 
by an average employee to recover the 
efficiency lost in the process of changing 
the production schedules (changes in pro- 
duction levels). 

the expected standby expenses of idle 
facility (the entire plant) for a period 
of one year. 

the upper limit on the (percent) level 

of overtime, for both types of work force, 
during any period. 

the upper limit on the (percent) level 

of employee idle time, during any period, 
for both types of work force. 

the rate for overtime (employee wage) 


premium. 
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the rate for shift (employee wage) 
premium. 

a factor to be applied to the minimum 
value amongst sales forecast to determine 
the lower limits for the production levels. 
Its value can be set anywhere between 0 
and 0.5 (where 0 < PMN < 1). 

the ratio of the total annual indirect 
material costs to the total annual direct 
material costs. 

the average hourly wage rate of the direct 
work force employees, including fringe 
benefits and all other types of allowances. 
the average hourly wage rate of the in- 
direct work force employess, including 
fringe benefits and all other types of 
allowances. 

the average hourly salary of a supervisor 
in the direct work force. 

the average hourly salary of the super- 
visor in the indirect work force. 

the rate for the administration costs, 

to be applied to the costs of production. 
the rate for the marketing and distribu- 
tion costs to be applied to the costs of 
production. 


the management's desired minimum after- 
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tax rate of return on its investments. 


C57) 246 this variable is the same as RMAR 
above. 
(58) PMX a factor between 1.3 and 2.0, to be 


applied to the maximum value among the 
sales projections to determine the upper 
limits for production levels during any 


period. 


Fs Data set - Typemo 


This data set contains the values of the penalty 
cost factors which are used during the evaluation of the 
objective function. These factors are applied to the cost 
of operations, to charge the penalties for the violation 
of the environmental constraints on the plant operations. 
The values of these factors are to be set between zero and 
one. Any of these factors, if the value is set equal to 
zero, makes the corresponding constraint totally ineffective. 
The penalty cost factors are: 

(1) F(1) for exceeding the upper limits on the 
levels of the inventory stocks. 

(2) F(2) for violating the constraints regarding 
the levels of stock shortages (expressed 
as fraction of the sales). 

(3) F(3) for the violation of the upper limit on 
the level of overtime in the direct work 
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for the violation of the upper limit on 
the level of employee idle time in the 
direct work force. 

for exceeding the upper limit on the level 
of overtime in the indirect work force. 
for exceeding the upper limit on the level 
of employee idle time. 

for exceeding the upper limit on the rate 
of hirings in the direct work force in 

any period. 

for exceeding the upper limit on the rate 
of employee terminations in the direct 
work force in any period. 

for exceeding the upper limit on the rate 
of hirings in the indirect work force 
during any period. 

for exceeding the upper limit on the rate 
of employee termination in the indirect 


work force during any period. 


The values of these factors used in the sample 


problem, shown in the data (Chapter 5), can be used as 


reference points for initiating the first run with the 


problem. 


But, in the subsequent runs these values should 


be set through trial and error procedures by increasing the 


values of those factors for which the corresponding con- 


straints are being violated. 
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G. Data Set - Type 6 


This data set consists of the values for all of the 
independent decision variables in the first planning hori- 
zon. These values are approximate estimates based upon 
past experience and form a starting solution base (the 
initial reference point) for intiating the optimization 
process In the model..ouLnisndata,.inge fach,wismeeguired only 
for a number of periods equal to the length of the planning 
horizon. But, it is necessary that these values be pro- 
vided for (NHM) periods (see master control data set for 
the definition of NHM) for convenience in changing the 
length of the planning horizon during different runs. The 


data in this set consists of the following information: 


(1) WD the direct work force level for (NHM) 
periods. 

(2y)7 611 the indirect work force level for (NHM) 
periods. 

(3) R the production levels for all the product 


lines for (NHM) periods, in separate sets 


(one set per product line). 


2.282  InputaDatal Eormats 


A. General Formats 
There are two types of cards in the data file: 
(a) Indicator cards, which specify the type of data which 
| isco Tol low. 


(b) Data cards, which contain the actual data values. 
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The indicator cards contain two integer numbers to 
specify the type of data cards that follow. These cards 


have the following format: 


ee ee Field-1 Field-2 
ems [sp ee 


where JA = the data set type number 


JF 


the number of (data) variables in this parti- 
cular data set. 

All data cards follow the same general formats, except 
the master control card (the first data card; the format 
rs-shownine rigqure 223). Each cardeis idivided into six 
fields as shown in Figure 2.2, one field for each numeric 
value. The number of fields assigned to each variable vary 
with the type of data variable. Each is explained in 


detail below. 


Bb. Notes onsinput Data 


(a) All numeric values appearing in the fields (1 through 
6) of the data cards (except the master control card) 
are defined by ten characters; these characters can 
include a decimal point. Specification of a sign is 
optional; if no sign is specified the plus (+) sign 


is implied. The position of the sign and the decimal 
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point is not fixed in these fields. However, for 
speed and the avoidance of scaling errors it is 
strongly recommended that the decimal point be placed 
in some fixed column in all cards and that there be 
no blanks between the sign and the most significant 
digit. 

The number may be followed by a decimal expon- 
ent, written as the letter E followed by a signed or 
unsigned one or two-digit integer. It is essential 
that there be no blanks after this integer number. 

(b) All the values appearing on the master control card 
and the indicator cards have to be integer numbers 
and must be placed as far to the right (in the fields) 
as possible so that there are no blank spaces left 
within the fields at the end of these numbers. 

ic) All the values of variables expressed as percents 
must be divided by 100 before entering them on the 
data cards. For instance, 56.3% will be entered as 
045.633. 

(d) The input information will be punched on the cards 
in the same order in which it has been listed in the 


input data list. 


C..) Description. of. -the.Data .Deck 


The organization of the data deck and the manner in 
which each input information will appear on the data cards 


is described in detail in this section. But, before pre- 


i ai, | 
no beaghe 6 Soa tte pach) oe! 
isttweeze et 31 negadnt dtpth-owl Yo ano bok “a ) 

xadmuq tspstah etat vests ednstd on od syed wit oa _ 
bis2 fortno? ietzsm oda no gotrseqqs, 2ouTey odd cy (a) 
ateduun vepetnt sd of syst ebise Vos 69 that wit iat 1A ! 
(2bfet? edt nt) digi sft) ot tot 26 basetg ad feum bis 
ste! seosge Anstd. on 2x6 s1eds Jails oe efdtezoq a6! - i 
— .etedmun ozedd to bas add se eblbett sitdomtdshe <= 4 ie 
atnsoteg 26 bazzavqxs eaidetvav to esulav sft, ie ae 
sit ng won? ontrigan welled Oat vd bobivtb od teum . 
as bovsdns od [ftw 26.08 a a on ee “ 
abve9 oft no bédanug sd ‘ote nots seneer dieil oat a 
ont oF pital), aon lane Senta ab. nine oat at oH 


33 


senting these details a few important points must be 

explained: 

(a) The statement, "one variable per card" to follow, 
means the variable has more than one value in the 
data and each value will appear one per field on the 
cand. 

(b) If the number of values of the variable(s) (being 
punched) exceed the number of fields available on 
the card, use more than one card. This applies to 
all the statements such as: "one variable per field 
or per card" in the explanations to follow. 

(c) Almost all the input variables have been referred to 
herein by their numbers (for identification purposes) 
which are as per the input data list (in Section 
28281)% 

All cards in the data deck will be punched and 
arranged as explained below: 

(a) Master control card - to have one variable per field. 

(b) Cards for data set - type 1. 

(i) Indicator card, to have a number 1 in column 
(3) and a number 7 in column (6). 

(ii) Data cards will have one variable per field for 
all the variables in this set. 

(c) Cards for data set - type 2. These cards will con- 
tain one set of cards per product line, where each 
set consists of: 


(i) Indicator card, to have a number 2 in column 
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(3) and a number 15 in columns (5) and (6). 

Data cards: 

- One variable per card for variables (1) and 
(ay. 

- One value per field for variable (3) (sales 
forecast data). 

- One variable per card for variables (4) and 
(5). 

- One variable per field for variables 6 to 15. 

for data set - type 3. 

Indicator card, to have a number 3 in column 

(3) and a number 17 in columns (5) and (6). 

Data cards will have one variable per card for 

all the 17 variables in this data set. 

for data set - type 4. 

Indicator card, to have a number 4 in column 

(3) and a number 58 in columns (5) and (6). 

Data cards will have one variable per field 

for all the 58 variables in this set. 

for data set - type 5. 

Indicator card to have a number 5 in column 

(3) and a number 10 in columns (5) and (6). 

Data cards will have one variable per field 

for all the 10 variables in this set. 

for data set - type 6. 

The indicator card, to have a number 6 in 


column (3) and a number 3 in column (6). 
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(ii) Data cards will have one variable per card for 
the variables (1) and (2) and one card per 
product line for the variable number (3). 
The listing of the input data deck, for the sample 
problem used herein to test the model (given in Appendix 2), 


can be used as a specimen for reference. 


2.3 Output (Planning) Information and Data Reporting 
Formats 


The aggregate production planning system generates 
all pertinent planning information necessary for decision 
making, which is reported on the output data formats shown 
iMermkOUres 2c. 4.97.5, .2.60, 2.7 dande2.6. tnese data report- 
ing formats are self-explanatory, hence, do not need to be 
explained here. 

The output information essentially consists of the 
following: 

(1) Aggregate manpower plan - includes both the direct and 
indirect work force levels for each period/month (in 
the plan) along with the respective period's capacity 
Utilization factor to indicate the actual utilization 
of the total (installed) capacity available through 
all the shifts the plant is (currently) operating. 

(2) Aggregate production plan - for each product line, 
consists of production and ending inventory levels, 
inventory shortage levels (expressed as a percent of 


the sales) for each period/month in the plan. 
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Work force performance analysis - which includes 
information on average employee's work-efficiency, 
overtime and employee idle time (the value will be 
negative in this case) by period/month. 

Production work force costs - gives a detailed break- 
down of the various manpower costs by period/month 
and by year. 

Total cost of plant operations - provides a complete 
analysis of the production and operating costs by 


period/month and by year. 
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AGGREGATE MANPOWER PLAN 


Direct Indirect | Capacity 
Period WoLetot Work Force Work Force Utilization 
Year /Month Ses Level Level Factor (3%) 


Figure 2.4: Reporting Format for the 
Aggregate Manpower Plan 
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CHAPTER 3 


BASIC FUNCTIONS USED IN THE MODEL 


All the basic production/cost functions used in the 
aggregate production planning model have been described 
inethis chapter. 7 Each function is explained herein; tn 
brief details, so as to provide the reader with a clear 
picture of the background, mathematical analysis and the 


underlying assumptions used in developing these functions. 


Sle EMO LOVEE Efficiency During Learning Periods 


The model for predicting employee efficiency (per- 
formance) during the learning period has been developed 
from the classic learning curve theory. This theory as 
applied to the learning of industrial jobs involving skills 
states that the production time declines by some constant 
percentage every time the number of units made is doubled. 


It is represented mathematically by a two parameter function: 
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where f(v) = the time required to make the vt" unit 
a = the time required to make the first uni 
b = a negative constant 
R = the learning rate constant, (0 < R “04 ) 
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3.1.1. Employee Efficiency 


An employee's efficiency at any point in time during 
the learning period can be defined as the ratio of the time 
required by an average experienced employee to complete 
job tasks to the time required by the employee under train- 
ing. The efficiency of an experienced employee is, there- 
fore, defined as 1.00. The efficiency of an inexperienced 
employee will take on values between O and 1.00 depending 
upon the employee's experience with the job. For the sake 
of analysis, all new employees are assumed to reach peak 
efficiency after a certain job experience, supposedly, the 
same period in each case (say after making x number of 


units or F(x) hours). Then, the employee efficiency (E) 


will be 
p= fH) 2 (x) 
f(v) vA 
for v = | 
dee oyt 
so 
eee eh 9.) 


eliminating v from equations [4] and [5] 
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For E = 1.00 and v = x the equation [6] will give 
Ve SC Bate [9] 


where F(x) = the time (in hours) to make x number of units 
or the total job experience required to 
attain the peak efficiency or in other words 


the total training period. 


3.1.2 Average Employee Efficiency 


The average employee efficiency (EA) over a certain 


length of the period, from F(v), to F(v), is given by 
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F(v)y 
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(B+1) (B+1) 

E {F(v) - F(v) } 

ER. = = . Ab head #2te [10] 
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The employee efficiency function is used in the aggre- 
gate production planning model to determine the average 
efficiency of the employees during the training/learning 
periods. This efficiency is an aggregate value for the 
total plant and is computed using the effective values for 


the training period and for the learning rate constant. 


3.2 Overhead Costs 

Overhead costs are those costs that cannot be clearly 
associated with particular operations or products and there- 
fore must be prorated by some method. There are many methods 
whereby overhead costs may be allocated. One should adopt 
a system that gives the desired accuracy with a minimum of 
effort on a continuing basis. In this model the overhead 
costs have been distributed on a man~-hour (equivalent of 


total production) basis: 


(or 


, et 
WF4#8), 
\ : 
: 


mnboveet \eniniss? sat entrud zaoyofqne ane ne sles 


< 5 
rl 4 
: 
7. 
2 
~ 
, | 
P 
S} : 
a 
‘ 
. 
tae! ’ 
os 
=» % 
, 
| 
é 
° 
) 
~ . 
i 4 
a 5 : 
d ' : a 
& 1 ‘ a 
Rs ) 

- nas ae —— 
‘y fT iran 
. ; 
hs ok 
psanveve 

. 
; ‘ 
va ~ te 


oo ets7 nateyeot + sit wt ba owas ware 


fonns. teAd eteod Sons ors 23900 bent v0 | 


vode anf ibeteootts se {om ebzo0 ‘bes 


. 
—- a a 
a aan: 


fe . j 
(va - «f v # Aue on : - 


a2u af noktomut yonetol cy 9 99Y 


gteb oF Tebom ontansta nots: vbang 9: BE 


fev sispeipes os ef vanatottte ett Ros” 
evtiostte oat? petew bssugqme > 2F one raat fi 


a ue 


,¥e 


et202. baad avd 


a 
ibovg ¥o-enortersqe rb Toot are wit ts fstace 
27 6 etait yond om omg2 yd benerorg pri 270" 


6 Asiw “oe 296 bentest ote. av 
eno : fis a ws , 
en? foebom ght! 1 vaband ntunt nas en ee 
i ry 


a 


ae 
vtups) ywod-em & 10 parva b bud oved 
4 n 


47 


_ Actual factory overhead expenses 
Overhead rate = Man-hours equivalent of the 


or 


where 


items: 


(1) 


(7) 
(8) 


total production 


pb TAF O 
OHR MHD 
OHR = the overhead rate 
AFOSR ipesactlal factory overheadmexpenses 
MHD = the direct labour man-hours equivalent 


of production. 


The overhead expenses consist of the following cost 


Depreciation expenses: 

(a) Land and buildings. 

(b) Plant equipment, machines and vehicles. 
(c) Office equipment. 

Rental expenses: 

(a) Machines and equipment. 

(b) Production space or facility. 

Repair and maintenance of buildings and office equip- 
ment. 

Utilities: water, heating and lighting, etc. 
Office supplies. 

Property taxes and insurance expenses. 
Engineering, design and development expenses. 


Scrap, rework and in-house material losses. 
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(9) Costs of outside services: contract services and 
travelling, etc. 
(10) Indirect overhead expenses: 
(a) Personnel support. 
(b) Training office expenses. 
(11) Operating expenses including repair and maintenance 
of plant equipment, machines and vehicles. 
(12) Indirect employees salaries and wages: 
(a) Inspectors. 
(b) Supervisors. 
(c) Repair and maintenance crews. 
(d) Clerical personnel: schedulers, time office 
clerks, truckers and shipping clerks, etc. 
(13) Indirect material costs. 


(14) Power and electricity expenses. 


3434-Costeof.+Rroduction al ime 


The cost of production time, used in the aggregate 
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production planning model, is defined as the fixed cost of 


operations per man-hour of the total work force and is 


computed using the following mathematical model: 


CFOH 
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where CPTH = the cost of production time per man-hour 
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CFOH = the fixed cost of operations, includes 
all the cost items from one to ten listed 
under overhead costs in Section 3.2 

WD = the direct work force level 

WI = the indirect work force level 

NHP = the number of regular hours of work per 


period. 


3.4 Direct Work Force Model 

The direct work force consists of all the production 
workers who are directly engaged in the manufacture of the 
product, i.e., it includes all workers such as, machine 
operators, fabricators and assemblers, etc. 

In production type industries, the direct work force 
requirements are usually estimated by various time study 
methods. This involves dividing operations into their 
basic elements, applying time factors to these elements and 
finally arriving at the total time for each operation. The 
total operation time consists of: the setup time (the time 
for preparing and setting up the machine for the next pro- 
duction run) and the cycle time or the run time (material 
handling time, machine time, and any delay times). The 
total production time, after accounting for man and machine 
inefficiencies, is then used to ascertain the total direct 
work force requirements using the following mathematical 


model: 
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J 
m clal 
WD” = yup * ) P x TPTU, 


where WD = the direct work force requirements 
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per unit of the product line (j) 

NHP = the number of regular hours of work, per 
period 

j = a subscript for product Tines:, 
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3.5 Indirect Work Force Model 
The indirect work force, in a production-type indus- 

try, consists of all categories of employees which are 

supplemental to the direct work force, such as: 

(1) Supervisory staff. 

(2,) Inspectors and quality control personnel. 

(3) Maintenance crews. 

(4) Planning, scheduling and routing personnel. 

(5) Inventory and stores personnel. 

(6) Other clerical start, Such as: truckers and shipping 
clerks, time office and other personnel directly 
engaged in services concerning production. 

Material handlers are also sometimes classified as 
indirect employees for the convenience of costing. But, 


in this model this category of employees is assumed to be a 
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direct work force class. 
Two methods have been built into the model to esti- 
mate the indirect work force requirements. Either method 


can be used depending upon the availability of the data. 


3.5.1 Regression Method 


In this method, the indirect work force is estimated 
as the function of the direct work force, the form of the 


regression equation is: 


Hine CWDS 


where WI” 


the indirect work force requirements 
WD = the direct work force level 
X = a constant equal to the minimum level of 
the indirect work force required to main- 
tain the plant operations at a zero direct 
work force level 
Y and Z are the regression coefficients and will 
be determined from historical data. 
The fortran code for developing this model is given 


in Appendix 4. 


3.5.2 Extrapolation Method 


This method should be used when sufficient historical 


data is not available to develop the regression model. 
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Here, the indirect work force is expressed as a linear 


function of direct work force and production levels. The 


Slope coefficients for this equation can be set using 


intuition/judgment estimates. The form of the mathematical 
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indirect work force requirements 
direct work force level 
actual direct work force requirements 


planned production level for the pro- 
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WI™ = (AS 
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where WI* = the 
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minimum indirect work force level 


required to run the plant for a single 
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required to maintain the plant for the 
second shift or subsequent shifts. 

AS = the number of supervisors required per 
employee of the direct work force 

Al, = the volume of inspection and quality control 
tasks in man-hours per unit of the product 
line (j) 

AM = the volume of maintenance work in man-hours 
per production man-hour 

AP = the volume of planning, scheduling, routing, 
inventory and stores work load in man-hours 
per production man-hour 

AC =» the number of all the clerical personnel 
of category 6 required per employee of the 
direct work force 

NHP = the total number of regular hours of work 


per period, 


3.6 Hiring and Training Cost Model 
The cost of hiring and training an employee consists 


of the following cost components. 


3.6.1 The Recruitment and Selection Costs 

(a) Advertising costs, including the fees paid to employ- 
ment agencies. 

(b) The costs of making selection decisions. 


(i) Salary and wage expenses of the interviewing 
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stath. 
(ii) The costs of materials and equipment used 
during assessment or tests. 
(iii) Employment office overheads. 
(iv) Other expenses, such as: physical examination 
fees, medical and clinical expenses and the 


costs of verifying the references. 


The Training Costs 


Ane =TIXed costs Of training; such lass sche CcOSst of 
training materials, books and reviews, etc.; and the 
training department overhead. 

The costs of productivity losses resulting from 
inefficiencies of the new employees, during the 
learning periods. These productivity losses are a 
maximum in the beginning of the training period and 
will gradually decrease as the employee's efficiency 
increases towards the end of the period. The cost 
of these losses is estimated as the cost of actual 


production time lost due to inefficiencies and is 


equal to 
CPTH(1 - EHA) 
where CPTH = the cost of the production time per 
man-hour 
EHA = 


the average efficiency of the employee 
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during the training period. 

(ici) The supervision costs, i.e., the cost of the time the 
Supervisor spends in directly training the employee 
which would not have to be spent if the employee were 
fully trained to perform his job. The amount of time 
spent will decrease as the efficiency increases. The 


actual amount of time spent for an hour is 
ST(1 - EHA) 


where ST = the average proportion of each hour 
the supervisor spends in directly 
training the new employee 
EHA = the average efficiency of the employee 
during the training period. 
Hence the total cost of hiring and training an 


employee is: 


CHTA Recruitment and selection costs + training costs 


CRSEPILGEFAMESC1 eS BAR) PPeCCPRHME S XPST) Pex NAT 


where CHTA the cost of hiring and training an employee 
CRS = the recruitment and selection costs per 


employee hired 


CPTH the cost of the production time per man-hour 
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CTFH = the fixed costs of training per trainee, 
per hour 

S = the hourly salary of the employee's super- 
visor 

ST = the average proportion of each hour the 
Supervisor spends in directly training the 
new employee during the entire training 
period 

EHA = the average efficiency of the employees 
during the training period 

NHT = the expected duration of the total training 


period in hours. 


3.7 Employee Termination Cost Model 


Termination costs within most companies are usually 
much more significant than they appear to be on the surface. 


These costs basically consist of the following. 


oo 721 Costlor Hiring aidalraininguan.emp lovee 


An employee is a valuable asset to the company. 
Therefore, termination of an employee means a net loss 


equal to the cost of hiring and training the employee. 


Sajne COS tt OT sexecessiver and rrequent lermina tions 


Very frequent and excessive termination rates create 
a feeling of job insecurity among the employees who in turn 


will find other employment and leave the company. Asa 
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consequence, the hiring and training costs increase and 
the plant production suffers in quality and output declines. 
This cost has been estimated empirically as the hiring and 
training cost times the amount by which the average rate 
of terminations over the past six periods/months exceeds a 
predetermined level. 


Hence, the total cost of terminating an employee is 


CEl”>= niring and training, cost 
+ cost of excessive terminations 


PLM - Wy 
= CHITA {1+ Tirayent: Woot 
(t-1) C 


where 


PLM = PTRW - PTR for WPTRWTSePTR 

= 0 otherwise 

PTRW = the average rate of terminations over the 
past six periods/months 

PTR = the maximum acceptable level of labour 
turnover 

CET = the cost of terminating an employee 

W, = the workforce level at the end of the 


current period 
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Wi) = the work force level at the end of the 


previous period. 


3.0 WorkeForce Costs 


The work force costs consist of the following costs. 


3.8.1 Regular Payroll Costs 


These costs are purely the salary and wage expenses 


of all the employees: 


CEDURE RX SNP eX W, 


where CLP = the regular payroll costs 
NHP = the number of regular hours of work per 
period/month 
W = the work force level during the period (t) 


ALR 


the average hourly wage rate per employee, 
inclusive of fringe benefits and all the 


allowances. 


3.8.¢c Overtime: Costs 

In this model it has been assumed that the overtime 
is used only to meet (temporary) shortages of work force, 
not other kinds such as the shortage of (installed) produc- 


tion capacity. Thus the overtime cost will be: 
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* 
COT = ROT x ALR x NHP x et - LI fe 02 w* > W 
= 0 otherwise 
where COT = the overtime costs 
ROT = the overtime rate premium 
ALR = the average hourly (gross) wage rate 
NHP = the number of regular hours of work per 
month/period 
W, = the work force level during the period 
We = the actual work force requirements for the 
period 
E. = the average efficiency of the employees 


during the period (t). 


3.8.3 Shitt. Premium Costs 


The shift premiums are paid to those employees who 


only work in shifts, but here, for convenience in analyzing 


the costs, it has been assumed that all employees wil] 


receive the shift premiums. The shift premium costs are: 


CSPo= (PSP = 1) x CLP for Si; > | 
a) FON yin See = ] 
where PSP = the rate for shift premiums, (PSP > 1) 
Ci Pema yetne 
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Si = the number of shifts the plant operates 


during the period (t). 


3.8.4 Hiring and Training Costs 


The employees hired normally have different skill/ 
experience backgrounds. Thus, the job learning/training 
time will, most of the time, be different in each case. 
Therefore, in order to maintain uniformity in the analysis, 
the skill or experience background of the new employees 
can be broadly classified into the following two categories: 
(a) The skilled or experienced class. 

(b) The unskilled or semi-skilled class. 
The hiring and training costs can be calculated as 


given below: 


CHT = CHTA{W, - We, _,)} for W, > We4_y) 


= 0 | otherwise 
where CHT = the total hiring and training costs 
CHTA = PHh-xeGHTAhetePHAayxICHEA2 
PH] = the expected proportion of hirings from 


the category (a) 

PH2 = the expected proportion of hirings from 
the category (b) 

CHTA] 


the hiring and training costs per employee 


from the category (a) (refer to Section 3.4) 
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the hiring and training costs per employee 
from the category (b) (refer to Section 3.4) 
the work force level at the end of the 
period (t) 

the work force level at the beginning of 


the period (t). 


3.8.5 Employee Termination Costs 


These costs are determined as follows: 


Cee a Cen(Ny ee Ty he) eeetoye Nam mh 


= 0 


where CETW = 
GEd 


otherwise 


the total employee termination costs 

the employee termination costs per 
employee (refer to Section 3.7 for a more 
complete explanation) 

the work force level at the end of the 
period (t) 

the work force level at the beginning of 


the period (t). 


The (total) work force cost is obtained by summing 


all the cost components listed above. Thus, the work force 


COSt (CW) 1s 


CW = Regular payroll costs + Overtime costs 


+ Shift. Premium costs + Hiring and training costs 
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+ Employee termination costs. 


3.9 Direct Material Costs 

Direct Materials are all those materials (used in 
manufacturing) which become a part of the final product 
and are involved in the operation in such a way that the 
material cost can be estimated. The major component of 
the direct material cost can be derived by completing a 
material take-off from the product blueprint(s). To this 
material cost one must add scrap material costs which are 
represented by two components: (1) scrap due to waste 
material, such as, slugs that are discarded from punching 
a hole in the component being manufactured and (2) compon- 
ents that have to be scrapped due to such factors as work- 
manship and machine misalignment. 


The mathematical model for the total direct material 


costs is: 
J 
CDMT, = iat Peg CDMU; 
P 
CD : P 

at MU X SAG 5 

where CDMU, Se eT Re 
SA. 
CDMT, = the total direct material costs for opera- 


tions during period (t) 
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(ep) 
a] 
= 
SS 
i 


the direct material costs per unit of the 


th th 


p product item in the j product line 


the expected direct material cost per 


QO 
=] 
= 
= 
Cu. 
It 


unit for the product line (j) 


SA: = the annual sales level for the pth pro- 


a 
duct item in the yeh 


product line 
PS. = the production level for the product line 


(j) during period (t). 


3.10 Indirect Material Costs 

Indirect materials are those materials which are 
critical to the operation (for manufacture) but do not 
become a part of the final product. This may include cutting 
lubricants, welding rods, perishable tools and clerical 
Supplies S*'Pecalse thdse Seosts avekdifficulll*to assess, it 
is charged to the product cost through an overhead distri- 


bution based upon direct material cost: 


CIM = (eam) x CDM 


where CIM = the cost of indirect materials 
CDM = the cost of direct materials 
CIMA = the total yearly indirect material pets 
CDMA = the total yearly direct material costs. 
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3.11 Inventory Carrying Costs 


Inventory carrying costs are the costs of maintaining 
inventory stocks in a production and/or distribution system. 
These costs basically consist of the following cost com- 
ponents: 

(als) Opportunity cost of capital invested in inventories. 
This cost is equal to the company's desired after-tax 
rate of return on its invested capital. 

C2) Shrinkage losses, the loss through breakage and pil- 
ferage at the storage site. 

(3) Loss through aging due to decomposition and/or decay 
of the products on shelves (with time) or decline in 
value of the products in storage due to obsolescence. 

(4) The cost of insurance and taxes. 

(Sy) Overhead or fixed costs. These include: 

(a) Depreciation and/or the rent of the storage 

Pacqeb Yu 

(b) Operating expenses: 

(i) Operating and maintaining the equipment. 
(ii) Utilities expenses. 
(fi7).,Sataries; and wages of the staff. 

All of these costs vary directly with the amount of 
investment in inventories, therefore, can be estimated.as 
rates per dollar value of the inventory stock. Hence, the 


total inventory cost function will be: 
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where ICR = the rate for inventory carrying costs as 


a percent of the total inventory value 


RMAR = the company's desired after-tax rate of 
return on its invested capital 

ISR = the rate for the shrinkage losses as a 
percent of the total inventory value 

ISOR = the rate of loss through decomposition and 
obsolescence as a percent of total inven- 
tory value 

IITR = rate of insurance and property tax expenses 
as a percent of the total inventory value 

IOFR = rate for overhead costs of inventories as 


a percent of total inventory value. 


In the criterion function, the inventory carrying 
costs serve to penalize the undesirable stock surpluses. 
Therefore, for use in the objective function, the inventory 
carrying costs are computed for the surplus stock levels 
only, i.e., for the number of units above the safety or 
minimum reserve stock level. But, for all other costing 
purposes the inventory carrying costs are computed corres- 


ponding to the actual stock levels. 


3.lc Inventory Snortages@osus 


Inventory shortage costs are the costs incurred as 
a consequence of stock out. That is, when the demand can- 
not be fully and immediately satisifed due to a stock short- 


age. The shortage costs consists of the following cost 
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components. 


312.1 oThe Cosititof Processing ‘the Requests for Improve- 
ments in Delivery Dates 


As a result of inventory shortages the deliveries 
for all the orders falling under shortages are postponed 
and the customers and/or distributors are informed about 
the expected delivery periods. Consequently, most of the 
customers and/or distributors will ask for better delivery 
dates from the manufacturer. Then, depending on the produc- 
tion schedule, the manufacturer is either able to give an 
improved delivery date to the satisfaction of some distri- 
butors or is forced to deny the request. But, this will 
incur costs (labour and communication charges, etc.) to 
process these requests whether or not the manufacturer is 
able to give an improved service. These costs of process- 


ing the requests (per unit short) are: 


where CPR = the average cost to process a request for 


a better delivery date 


M = the average number of stock-outs per year 
AOS = the average order size, number of units 

per order 
F = the expected annual costs of processing the 


requests for better delivery periods. 
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3312.2  The?Cost-ofelkost Sales. 

Some sales may be lest on orders met by the stock out 
condition, because the customers and/or distributors to 
whom the late deliveries are not acceptable may cancel their 
orders for this particular period (for this particular item 
only). This loss will be equal to the amount of after-tax 


profit foregone on these sales. 


3.12.3 Goodwill Losses 

Some customers/distributors who have been given the 
delayed delivery dates may terminate their business with 
the manufacturer, supposedly for the item short only. This 
means a loss of all the future sales, expected through 
these clients. This loss can be estimated to be equal to 
the present worth of these sales. For the sake of analysis 
it is assumed that the distributor who terminates his busi- 
ness with the manufacturer for this item, generally speak- 
ing, would resume business after a certain length of time 
(assumed to be the same in each case, say (q) years). 

LiaP 


P. and P are the probabilities associated 


WY 
with the three cases discussed above, then the expected 


COSt Of GCacn UnNiL SNOYUTIWTL) Des 


CINSA = Ps x Cost of processing the requests for better 
delivery periods + P. x Costs of “ost sales 


+ P, x Cost of goodwill losses per unit 
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: F g x PQI x ASC 
UPA TROS) OTe tay xe CPU) tenet AOS ) 


where CINSA 


the cost of shortages per unit short 
ie = the expected yearly cost of processing 
the requests for better delivery periods 


(includes labour, communication charges, 


etc.) 
M = the average number of stock outs per year 
AOS = average order size, number of units per 
order 
g = the rate of (after-tax) profit earned on 


the cost of sales 


CEE = the cost of the product (per unit) before 
Sales 
ASC = the average yearly sales per distributor 


(to be estimated in terms of the cost of 
sales) 


PQI = present value factor. 


The inventory shortage costs are subject to estimates 
of probabilities that the actual demand quantities will 
exceed stock quantities during the future periods. These 
probabilities are quite difficult to estimate and would 
complicate the analysis. Therefore, as an alternative, the 
method used in the model for estimating this cost presumes 
that the loss is sustained whenever the actual stock level 


drops below the safety or minimum reserve stock level. The 
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number of units below the safety stock level multiplied 
by the dollar contribution (loss) per unit would give a 


reasonable estimate of the cost of stock shortages. 


3.13 Shift Start-Up and Shut-Down Costs 


Start-up or shut-down of an extra shift imposes addi- 
tional costs which are incurred at the beginning of the 
period of program change. These costs basically consist 
of the following cost components: 

(1) The cost of additional planning work: to reschedule 
the man-machine loading; to make adjustments in the 
production runs; and to work out the new manpower 
requirements (for hiring/layoffs and overtime). | 

K25) The cost of productivity lost due to work-inefficiencies 
resulting from the redistribution of the work force 
to suit the new shift arrangement. This cost is esti- 
mated as the cost of the useful production time lost 
in the inefficiencies. 

The following assumptions have been used in developing 
the shift start-up and shut-down costs model: 

() The minimum production capacity (man-hours of direct 
work force) requirements for starting or continuing 
the additional shift are assumed to be greater than 
that available through use of a maximum amount of 
overtime. 

(2) The number of employees redistributed for implement- 


ing the new shift program is equal to half the minimum 
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number of employees required to start an additional 
shift. 

(3) All the additional planning tasks, resulting from the 
start-up or the shut-down of an additional shift, 
Supposedly would make half a days work for the plann- 
ing and scheduling employees whose number is equal 
to the number of total planning and scheduling 
employees times the (minimum) proportion by which 
the existing production capacity will be increased 
or decreased. 

The mathematical model for estimating these costs is 


given below: 


CSS = Cost of extra planning work 

TaCOSUsOh DROCUCGTVILy alos t 

is MHDL2 

=) 4.0 .eAURT x PRY x Wi =x (WHDU ) 

+ 0.25(1 - ES) x FD(S) x CPTH x (MAL 7 MBI 2) 

where GS$ = the additional shift start-up or shut- 
down cost 
ALRI = the average hourly wage rate for the 


indirect work force employees 
ee = the ratio of the planning and scheduling 


employees to the total indirect employee 
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work force 

ES = the level of work-efficiency at the 
beginning of the period of the program 
change (i.e., after a reorganization of 


the operations) 


FD(S) the time in hours required to regain the 
peak work-efficiency (100%) 


CPTH 


the cost of the production time per 
man-hour 

NHP = the number of regular hours of work per 
period/month, per employee 


MHDL2 


the minimum production capacity (man- 
hours of the direct work force) require- 
ments for the start-up or shut-down of 
an additional shift 


MHIL2 


the minimum requirements for the indirect 
work force (man-hours) in the second and 
each subsequent shift 


MHDU1 the maximum production capacity (man- 


hours of direct work force) available in 
tharkinrse srpvet 
WI] = the indirect work force requirements for 


Che Gf1 esit Tshifematmiull production capacity. 


Rie pee Costs of *Fluctuations “in eProduction Level’s 
The costs of fluctuations in the production levels 


may be classified into the following two groups: 
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(1) The expenses are those incurred only once at the 


beginning of the period of the new program. These 


include: 


(a) 
(b) 
(c) 
(d) 


(2y) The 
the 
(a) 
(b) 
(c) 


Hiring and training costs. 

Employee termination costs. 

Shift start-up and shut-down costs. 

The cost of productivity lost due to the work- 
inefficiencies resulting from the redistribution 
of the work force in order to adjust to (a) the 
changes in the work force levels (hirings and 
terminations) and/or (b) the changes made in 

the set-up of the production lines or the work 
centres. This cost will be estimated as 
explained in the previous section. 

The cost of additional planning work resulting 
from the changes in the production programs. 
This includes the tasks such as: rescheduling 
the man-machine loading and the production runs; 


and the manpower planning. 


expenses that may be imposed for the duration of 


period of the new program. This includes: 


The overtime expenses. 
The costs of shift premiums. 


The standby expenses of idle facilities. 


All of the cost items except 1({c), 1(d) and 1(e) have 


been considered in the aggregate production planning model 


separately for the sake of convenience in the analysis. 
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Therefore, the costs of fluctuations in production levels 

Williconsist of the cost: wtemset (c)<> 1(d) Candeive) only: 

In order to facilitate the mathematical analysis of 
this model the following set of assumptions has been used: 
As) The total number of employees redistributed (as a 

part of the process of reorganizing the operations): 

(ayy. ine caserog hiningseor terminations. isu equal’ to 

half the number of employees hired or terminated; 

(b) in case of a change in the production level, 

their proportion is equal to the proportion by 
which the production level has changed. 

(2s) The additional planning tasks are assumed to be three 
hours of work for the planning and scheduling 
employees whose number is equal to the number of 
the total planning and scheduling employees times the 
proportion by which the production level has changed. 
The mathematical model for these costs is as given 


below: 


CCP, = Costs of inefficiencies resulting from changes in 
production levels + costs due to additional planning 


tasks + shift start-up and shut-down costs. 


= 0.5{0.5(AWD, + AWI,) + PMHD,(WD(, 4) + WI¢4_4))} 
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|WD, - WD¢4_y)! 
De teden cal), 


) |MHD, - MHD, 


| 
j=] t-1)j 


arsupsienipt for periods/months, t=1,2,3,... 


AMSUDS CRIP terol PrOdUGtMINnNeS wij 1,2 poees « 
the costs of fluctuations in production 
levels 

the direct work force level 

the indirect work force level 

the production of the product line (j), 
man-hours of direct work force 

the cost of the production time per man- 
hour 

the shift start-up and shut-down costs 
per shift started or shut-down 

the number of additional shifts started 
or shut-down 

the number of regular hours of work per 
employee, per period 

the average hourly wage rate for the 
indirect work force employees 

the ratio of the production planning and 


scheduling employees to the total indirect 
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employee work force 

the expected work-efficiency level in the 
beginning of the period of the new program 
(i.e., after reorganizing the operations) 
the expected time in hours required to 
regain the peak work-efficiency level 
(100%) after a change in the production 


program. 


3-15 Unbt Manutactunings€ost 


The manufacturing cost consists of the direct material 


costs, the direct labour costs and the overhead costs. 


The mathematical model is as given below: 


CPU, 


where CPU, 
CDMU. 

: J 
TeIU. 

J 

COHR 


CDWH 


CDMU, sy TPTU.(CDWH + COHR) 


the manufacturing cost per unit 

the direct material costs per unit 

the production time, man-hours per unit 
the overhead rate, dollars per production 
man-hour 

the labour rate, dollars per production 
man-hour 

a subseript tor the product: limes, 


Pa ese a Us 


of beytupes 2% 
level cone tat eo 
-fottauborg odj nt seuss 5 peo cs . i 2 
Se ee offs et ‘Paes: nan a ON Te aay) ; - 
So pee ih ays oe # o _ O° tae 7 


ms 


‘nip tay a A Coy t 


-ed202 dahisereis sunishaaieeiaal nC a a 
oe, twoked, peytp, 28, 2h Isbom feottsmedism bi 


ae] > a) Tae) OY e (2% oe ae ech f Lah ai | : 
PVRHOO weer i i tam 
HEF Mint? Se a ead beeen me 


sau v9q $403 patawtostunsm aay uaa | ened 
““ghau 199 23200 fetvedem soarth and <a 
ttnv veq 24vod-nem ,smit natisuborg odd 
nottouborg San avst toh os r bach odd . 
notsoubowg req ite swodal co he! 
v¢ ; . 

Senet ee 
revel Gir eanae << tune th 
bad pnianalo haan aaa a 


a shin a> as nay 


a 
Cad 
= 
a 
~ 
s 


| : 4" 
- 


76— 


3.16 Total Cost of Operations 


The total cost of operations is obtained by summing 
the following set of costs: 

Ol) Direct work force costs - are the labour costs associ- 
ated with specific machine operators. 

G2) Direct material costs - are the costs of all those 
materials which are directly consumed in the manufac- 
ture of the product (including scrap). 

(3) Overhead costs - includes such costs as rents, deprecia- 
tion, taxes, maintenance, indirect materials, indirect 
work force, power and electricity, heat, light, water 
and utilities. 

(4) Inventory carrying costs - are the costs of maintain- 
ing (carrying) the inventory stocks. 

(5) Inventory shortage costs - are the costs of expected 
losses that will be incurred when an order is met by 
a condition of stock shortage. These are intangible 
costs and are assessed through probabilistic estimates 
(very approximate). Even then, these costs are of 
great significance in the analysis. 

(6) Costs of fluctuations in production levels - are the 
costs of productivity lost due to the work-ineffici- 
encies resulting from the redistribution of the work 
force and reorganization of operations necessary for 
a new production program. 

Cre, Marketing and distribution costs - consist of such 


jtems as salaries and commissions for sales personnel, 
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advertising, entertainment and travel expenses. 

(8) Administration costs - consist of such cost items as 
salaries and allowances of management personnel, 

wages of clerical and accounting staff, office rent 

and supplies and utilities. 

The last two costs (numbers 7 and 8) will be estimated 
through fixed rates applied to the total production costs 
(the sum total of the first six costs). The costs of opera- 
tions, computed, herein, do not include any profits. But, 
the gross profits can be added separately, to arrive at the 
cost of sales, through a rate applied to the cost of opera- 
tions. It may be mentioned here that the system routines 
are designed to calculate the cost of operations only, not 


the cost of sales. 
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CHAPTER 4 


OPTIMIZATION MODEL 


The optimization of the decisions in SDR involves 
the act of selecting the "best" solution from among a set 
of feasible alternatives. The entire optimization model 
has been sub-divided into: a planning cost model - to com- 
pute the production costs and the criterion function; a 
een! of constraints to define the decision space; and an 
optimization technique - a computer search routine to 
select the best feasible solution. Each is discussed below 


in details. 


4.1 Planning Cost Model 


The planning cost model determines the cost of opera- 
tions over the planning horizon and, hence, the objective 
function for evaluating the performance of the system 
(firm). The flow chart shown in Figure 4.1 provides a 
detailed description of the planning cost model. The 
Fortran equations for the model are contained in the docu- 
mented listing of the subroutine FCT], given in Appendix 2. 

The model, virtually operates the system (firm) for 
each month over the planning horizon based upon the decisions 
supplied by the pattern search routine. This operation 


Consists of costing out the impact. of the decisions im each 
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ENTRY | @ 


INITIALIZE 


COMPUTE (FHESDIRECT 
WORK FORCE REQUIREMENTS 


DETERMINE THE NUMBER OF SHTETS 
REQUIRED TOSMEET THE SSALES 
COMPUTE THE INDIRECT 
WORK FORCE REQUIREMENTS ; 
COMPUTE OAVERAGE TERE ICTENCYIUEVELS 
POROTHESBOTHOTYPES! OF) WORK. FORCES 
COMPUTE DIRECT AND INDIRECT 
WORK FORCE PAYROLL COSTS 


COMPUTE DIRECT AND INDIRECT 
WORK FORCE SHIFT PREMIUM COSTS 


IS THERE ANY (EMPLOYEE) YES 
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NO 


COMPUTE DIRECT AND INDIRECT WORK FORCE 
HIRING AND TRAINING AND TERMINATION COSTS 


YESS IS RATE OF HIRING/ 
TERMINATION WITHIN LIMITS? 


NO 


COMPUTE SPENABEETGOSTS 


Figure 4.1: Flow Chart of the Planning Cost Model 
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Figure 4.1: (Continued) 
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month and computing the various cost penalties for violat- 
ing the constraints on the production parameters. The 
critical operating costs are then totaled (for all the 
months in the planning horizon) and the penalty costs 
added to get the value of the objective function which is 
then returned to the search routine where it is assessed 
to establish the next move for a new trial solution. After 
the search terminates the planning cost model computes all 
the planning data required for reporting the results of 
this period's decisions, 
The objective function consists of the following cost 
components: 
(az) Work force costs: 
(a) The direct work force costs. 
(b) The indirect work force costs. 
(2) Inventory costs: 
(a) The inventory carrying costs. 
(b) The inventory shortages costs. 
(3) Costs of fluctuations in production levels. 
(4) Penalty costs - charged for violating the constraints 
on the production parameters (discussed below). 
These are charged in such a way that the amount of 
the penalty cost is proportional to the amount by 


which the constraint is exceeded. 


4.2 Set of Constraints 


The constraints used in the aggregate production 
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planning model are classified as follows: 


(1) 


The range constraints for the decision vectors (pro- 

duction, direct work force and indirect work force). 

These constraints are used within the search routine 

(the subroutine PATS) to define the boundaries of the 

decision space for the problem. The search routine 

has built-in provisions to ensure that these con- 
straints are not violated. 

The constraints imposed by the operating environment 

of the firm, includes: 

(a) The maximum limit on the amount of overtime as 
a percent of the work force level, for both 
types of work forces. 

(b) The maximum limit on the amount of idle time 
as a percent of the work force level, for both 
types of work forces. 

(c) The upper limit on the rate of hiring in any 
period, for both types of work forces. 

(d) The upper limit on the rate of employee termina- 
tion during any period, for both types of work 
forces. 

(e) The upper limits on the levels of the month end- 
ing inventory stocks of each product line for 
each year in the total planning period. 

(f) The upper limits on the levels of stock short- 
ages (expressed as a percent of sales). for each 


product line. 
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These constraints control the independent decision 
variables indirectly. A heavy cost penalty is charged to 
the objective function costs each time any restriction is 
exceeded. This provides a check on the activity levels of 
the variables. The size of the penalty cost for each con- 
straint is determined by the value of its corresponding 
penalty cost factor. The values of these factors can be 
set anywhere between zero and one in the input data (data 
set-type 5). These constraints, if desired, can be made 
inoperative by setting the value of the corresponding 
penalty cost factor equal to zero. Thus, in effect, these 


constraints are optional. 


4.3 Optimization Method 
Under the SDR formulation of the aggregate production 


planning model the optimization function is performed 
through a direct computer search method using a "Pattern 
Search Technique". The optimization method involves defin- 
ing the problem as a multi-stage decision system. A stage 
represents the cost structure of the firm (in the form of 

a cost subroutine) at a particular point in time when the 
decisions are to be made, which is usually, but not necess- 
arily, monthly. The decisions for each stage are optimized 
with respect to overlapping n-period planning horizon. 

This process progresses gradually, through all the N- 
stages (periods/months) of the system, optimizing one stage 


at a time, in succession. Figure 4.2, taken from reference 
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[43], illustrates the complete SDR optimization process 


for a simple multi-stage decision model. 


The 
S = 


The 


iteration 


notations used in Figure 4.2 are defined below: 

a (i) component state vector, describes the state 
of the system at the beginning or at the end of 
the period of stage transformation and serves to 
carry information into or out of the stage-system. 
a (k) “component vector, consists of those factors 
which affect the stage return (vy) and the state 
vector (S). It contains both the demand fore- 
casts and the production/cost parameters. 

a (j) component vector, transmits the selected 
(feasible) trial values of the independent vari- 
ables from the search routine to stage systems 

for evaluating the objective function (the stage 
returns). 

the stage return - a scaler variable, measures 

the utility of the stage as a single valued func- 
tion of the input state, parameter and decision 
vectors. "The*sum or “the returns (or-all>tne-(n) 
stages in the planning horizon forms the criterion 
function. 

optimization’ process’ 1s. iterative*®in nature. Each 


involves, formulating a feasible trial solution 


by specifying sehective Values to*the decision’ variabies 


in the search routine and evaluating the criterion function 


(the stage-returns) for the planning horizon in the planning 
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cost subroutine. The value of the criterion function is 
then compared to the best previous value. If an improve- 
ment is observed the value (the trial solution) is accepted; 
if not it is rejected. Accordingly, the search proceeds 
in a systematic and sequential manner, directed by the 
search logic (Pattern Search) looking for better solution 
points in the solution space defined by the production/cost 
parameters. This iterative process continues, for each 
period/month, for a finite set of iterations and the search 
is terminated when the improvement observed is quite 
negligible, less than a specified value; or when a pre- 
determined number of evaluations have been made. 

Following the optimization of the first period/month 
the decisions are reviewed and implemented and a search 
for the optimization of the next period's decisions is 
started with Das D3, Daseees Dd, as the initial starting 
decision vectors. This process is repeated for all sub- 
sequent periods in the entire planning span. Figure 4.3, 
taken from reference [43], explains the information flow 
in the SDR monthly updating cycle. It will be observed 
here that the search for the first period's optimization 
will start with trial values (supplied as input in the 
data) for all the independent decision variables in the 
planning horizon. But, for the optimization of each sub- 
sequent period the routine needs to add trial values only 
for last period in the planning horizon. The demand fore- 


cast and the production/cost data will be required for the 
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(N + n - 1) periods, where N is the total number of periods 
in the entire planning span and (n) is the number of 
periods in a planning horizon. 

Pattern Search Code. The pattern search code used 
in the SDR is a modified version of the original code 
developed by Hooke and Jeeves [23]. The strategy is com- 
pletely heuristic in nature and is not based on any specific 
model of the response surface. The search consists of 
looking for better solution points in the n-dimensional 
decision space. A move from one point to another in the 
decision space is called a "pattern move" and is based, in 
magnitude and direction, on the information about the local 
characteristics of the response surface supplied by the 
"exploratory search". A master control strategy within the 
search code guides the pattern move and exploratory search 
logic in coordinating their actions and conducts appropri- 
ate tests for convergence and repeated pattern move failures. 

The exploratory search logic is explained in detail 
in Figure 4.4 taken from reference [43]. Vector X contains 
n coordinate values of the base point and vector D contains 
the corresponding step sizes computed during the last 
exploratory search. In the search process the behavior 
of the function is inferred from the success or failure 
of the criterion function evaluations as the value of each 
coordinate is modified. This success/failure information 
is used to modify the values of the step sizes for future 


exploratory searches and to formulate a pattern for the 
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Figure 4.4: Descriptive Flow Chart of the Exploratory Search Logic 
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the Pattern Move Logic 
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next move in the probable direction of success. 

Figure 4.5, taken from reference [43], shows the 
details of the pattern move logic. The pattern move logic 
serves to move the search to a new temporary base point by 
changing the coordinate values by an amount equal to the 
difference between the new base point and the previous 
base point. If the local exploratory search at the new 
temporary base point indicates a success, the base point 
is accepted and the search proceeds. In case of anticipated 
success for future moves the built-in acceleration proce- 
dures speeds up the search process by increasing the step 
sizes to cover larger distances during each move. If the 
pattern move is a failure, the temporary base point is dis- 
carded and the pattern is destroyed. The search returns 
to the previous successful base point and attempts to esta- 
blish a new pattern. If this fails, the step sizes are 
reduced and another exploratory search is conducted. This 
process continues until a new pattern is established and 
the search moves on, or the step sizes fall below a pre- 
set minimum, and the search terminates. 

The Fortran version of the pattern search code, used 
here, has some built-in procedures for independently con- 
trolling the step size of each variable as well as some 
rather sophisticated search termination logic. A fully 
documented listing of the pattern search program is given 


in the Appendix. 
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CHAPTER 5 


TESTING THE MODEL 


The aggregate production planning model was tested 
on a hypothetical problem in order to evaluate its per- 
formance. The sample problem was designed to represent 
a realistic factory environment so that the results would 
enable one to assess the operational capabilities of the 
model. The problem used for testing is described in sec- 
tions. leand the results are contained ingsection S.200Ff 


this chapter. 


5.1 Sample Problem Formulation 


The sample problem is a hypothetical, but neverthe- 
less realistic model of a manufacturing firm making two 
Phines*®oreproducts inga@batch production type facility. The 
firm presently operates two shifts but has an installed 
capacity for a three shift operation. The overtime is 
used to meet the temporary extra capacity requirements 
resulting from the fluctuations in sales. The company hires 
new emptoyees for both the direct and indirect types of 
work force from two categories based upon the skill back- 
ground: (1) the skilled and experienced class; and (2) semi- 
skilled or unskilled class. The new employees undergo a 


formal training program in order to familiarize themselves 
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with the job and improve their performance on the job. 

A planning span of two years is selected with fluctu- 
ating sales demand. Decisions are made by period (each 
period consisting of 400 hours, 2,000 hours per year) and 
are expressed in terms of direct work force size, indirect 
work force size and production levels for both the product 
lines. A planning horizon of four periods is used for 
optimizing each period's decisions. 

The input data used here for testing, listed below, 
is estimated to be reasonably representative of a practical 
Situation. In the list, each (data) value is written 
against the corresponding variable name and in the same 


order as is described in Section 2.2.1. 


3. bo leginpurddata 


A. Master Control Data 


(3) ~ NOY 2 (7)  “JRROPX 1 
(2) NPY 5 (8) NWI 0 
(3) NH 4 (9) MX 3 
(4) NHM 6 (10) KPRESS 1 
(5) Jd 2 CODE 3,000 
(6)  NHP 400 


B. Data Set - Type 1 


Beginning work force levels. 
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Operating data during the period prior to the 


start of the plan. 


(2) eV0ED 0.95 
(aye OE 1 0297 
(5)  OTDPI 0.10 
(6 pee OnLPT 0.10 
(7) anSO 2 


C. (Dataesete- Type 2 


Previous period's production levels 
Goclor tortie start..of ‘the plan) 


(gl) n0P 


Beginning inventory levels 
(number of units). 


(23) OVI 
3) Sales forecast (number of units) 
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Limits on the levels of inventory 


stocks. 
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Other production/cost data regard- 
ing the product lines. 
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Data set = Typegs 


Data which is different for each year. 
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Data Set - Type 4 


Data regarding the operating environment. 
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RentDatatSetre Type 5 


Penalty cost factors. 


Ch) enF (a) 0.10 
f2vet Fi2) 0.07 
(3)cluF(3) 0.015 
(4) F(4) 0.010 
(5) FS) 0.015 
(6) F«6) 0.015 
(7) E(#) 0.001 
(8)  F(8) 0.002 
(9) —F(99 0.005 
(10) F(10) . 0.010 


G.| DatalSet eflype 6 


Initial starting solution. 


(1) WD 108 95 98 93 90 90 
(2) WI 7] 65 7A 68 68 68 
(3) R(1) 348 447 44] 400 400 400 

R(2) 830 620 750 750 800 800 


5.2 Test Results 

The aggregate production planning model was thoroughly 
tested using hand calculated test cases prior to actual 
SDR optimization. At the end of this testing phase the 
model was operated over a ten period planning span with the 
sample problem. Several test runs were made to assess the 


performance of the model under different operating environ- 
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ments. The results for a test run with the sample data 

are presentedirvin, tab esicSi ih, thicZ yc. 37 Th, whe land oi36, at 

the end of this chapter. The operating behaviour of the 

model observed during the test runs led to the following 
conclusions: 

(21) Several runs are required to obtain a satisfactory 
solution to the problem. Because, each time the 
initial starting solution is replaced with the solu- 
tion from the previous run, the model results in a 
better solution. The values of the penalty cost 
factors in the input data set - type 5, also need 
to be reset, sometimes during each run. This is a 
trial and error procedure and involves increasing the 
values of the factors for which the corresponding 
constraints have been violated. 

(2) the independent decision variables are quite inter- 
dependent in the criterion function cost model. Their 
interrelationship imposes additional (indirect) 
constraints, which sometimes suppresses a solution 
move even if it was in the right direction. This 
problem can be overcome by changing the values of 
some of the penalty cost coefficients as explained 
above. 

(:3)) The quality of the solution and the speed of the 
optimizing process improves considerably with a 
reduction in the number of penalty cost constraints 


(listed in Section 4.2). The model has a built-in 


es 


ef? Fo +woty dod nieve 


wo Tot ont of bef zav 


va a s . 7 


one cee t ° 
Mg nite 
o on UJ \p - 


etetise 6 tase of bot tupon. ons anu Ts4ava 
git 


vzov Tsbom ont coy euotvang ais on a nore 


F cc 1) ' Ae 
j} dosed ,82¥eos a mel dorq sf2 0 
igtw boeosTqat ef Pout oe entire te fotgtm 


[snoq afd To esvtev sat seh: Ne 
ie gavt = see sited sug> “edt abe F 
.oun does pnfvud comtveno® «2204 od oe 
zovfovat bas stubs9019. TONS ne taint 
| i 
10> aft dotdw wot e102 387, ont to ul PBv 
i , ts 

.batsfoty: need oved asn te: an 

7 ee ——) 
4g eeldsirev aotztosb tn pennant 
S202 wnetsonu? noties t49 ond af 30 obns 


5 


4) Se eszoqmt athenotantorss 
ttoovtbh dapts edt at 26w- ant aga py 
sit pateneda yd emoa19vO od 0 donq 
ax stnatatttees ¢z00 wt! jonsg att see 


| ine, ie “tees oem ods. ar 
aa i= 
bsoqe oft bine norsutos “ pre) aut, 


dttw yldsvebPeneo aovorgm sae ste ttq 


2intent2zno> 


i” 


oD vttensg 7 


9d im 9 pie aap 


o 5 oe oe hel a 
est fabom ent 6 nok 4082 af 


Te i y a ; ’ ie i 1 r, 
a aed 


101 


provision of using all or none of these constraints. 
Therefore, one should include in the problem only 


those constraints that are necessary and important. 
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VABC ESS ai 


AGGREGATE MANPOWER PLAN 


Direct Indirect Capacity 
Period No. of work Force Work Force Utilization 
Year /Month Shifts Level Level Factor (%) 
1 0 2 8] 58 -- 
1 1 2 84 67 97202 
1 2 2 84. 67 Siso72 
1 3 2 84 67 94.33 
1 a Z 84 69 95.74 
1 5 2 84 71 109.33 
Z 6 2 84 7] 108.70 
2 7 2 84 68 94.64 
2 8 2 84 70 109 G22 
2 9 2 84 70 TOSS 
2 10 2 84 70 107.87 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


The aggregate production planning system designed 
herein is an attempt at improving the overall effectiveness 
of planning in small manufacturing companies. The aim is 
to enable small business to be more competitive. 

The aggregate production planning system, in broad 
terms, iS a time sharing type of computer-based management 
system. The entire system consists of three sub-systems: 
(i) an information management system - which collects and 
prepares all the input information; (ii) an optimization 
system - which determines optimal planning decisions; ane 
(iii) a planning data generating system - which develops 
and prints out, in specified formats, all the necessary 
planning information pase upon the "optimal" decisions. 

The computational model for the system was formulated 
within the framework of the SDR technique. SDR in the 
application is used in a similar manner to any other 
analytical method (e.g., Linear programming, Linear Decision 
Rule, etc.) for solving the operations planning problem. 
Under the SDR approach the problem is defined in terms of 
a multi-stage decision model which is then solved using a 
computer search technique (a direct search method for 


optimization). The decisions for the entire planning period 
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are determined stage by stage, each stage being "optimized" 
over a multi-period, overlapping planning horizon. 

The SDR technique has certain definite advantages 
over other analytical methods. Unlike other methods SDR 
places fewer constraints on the model builder. The cost 
functions need not be linear or quadratic and continuous 
but can be in a form which seems to duplicate reality. 
Realistic constraints on such factors as: normal capacity, 
overtime, hiring and termination rates, inventories and 
storage capacity can be handled by the model. In other 
words, under SDR the model can duplicate a reasonably com- 
plete phenomenon of actual factory operations, where the 
optimization process involves operating the model for each 
decision stage to evaluate the cost impact of the period's 
decisions over a certain number of periods (planning hori- 
zon). This process is repeated within the search to obtain 
the "best" (near optimal) decisions for the period and 
covers the entire planning span stage by stage. 

Use of this aggregate production planning system 
requires some background knowledge of the subroutine Pats - 
the search routine used in SDR. Because, with the present 
form of the subroutine Pats the quality of the solution 
and the speed of optimization depends somewhat upon the 
values of the step functions DEL, BET and ALP. These vari- 
ables are used for search and acceleration procedures in 
the search routine and determine the size and nature of 


the step by which the value of a particular independent 
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variable will be incremented during the formulation of a 
trial solution. Therefore, it becomes somewhat mandatory 
to offer the system for service through some service 
Organization where the officials can make themselves 
familiar with the necessary background required for using 
the system. This factor, as it presently exists, is a weak 
spot within SDR. Future study of the model will be necess- 
ary to overcome this deficiency. 

The test runs with the sample problem alone are not 
sufficient for a complete assessment of the practical cap- 
abilities of the model. The model should be thoroughly 
tested in an actual factory planning situation within two 
or three small manufacturing companies. These tests should 
involve testing over some past period so as to compare the 
test results with actual (past) management decisions. This 
testing would probably reveal the areas needing improvements 
and modification. After necessary refinements a somewhat 
larger cross-section of industry should be approached and 
further testing of the system should be initiated. UIti- 
mately, when fully satisfied with the performance the 
system may be introduced to interested companies through 


some central service organization. 


6.1 .TopicsntorshutunelResearch 


The aggregate production planning system can be 
expanded to include the following extensions in order to 


make the system more sophisticated and versatile. 
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Stochastic extension. The work carried out herein 

is based on the deterministic, or "perfect" fore- 
casts and cost models. A logical next step is to 
extend the model to a stochastic decision model to 
handle indeterministic forecasts and cost variables. 
Disaggregation extension. The aggregate production 
planning system, at present, handles all the produc- 
tion planning decision details to the extent of pro- 
duct lines at the overall plant level. The system 
can be extended to a disaggregate planning system to 
deal with the decision details at the level of pro- 
ducts and departments/sections. 

Extension to include sensitivity analysis. The pre- 
sent design of the system does not have any provision 
for doing sensitivity analysis regarding planning 
decisions. This extension would include designing 
routines for doing the sensitivity analysis with 
respect to various cost parameters such as labour 
COStS: spayrolis, overtime, Niring-and training) costs 
and inventory costs. | 

Extension for a conversational mode. Presently the 
system is meant for use in a batch type operating 
mode. It should be modified to operate in some type 
of conversational mode so that the user is able to 
interact directly with the system in his own language. 
By building an interactive/interrogative mode into 


the system it would be possible to enable the user 
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to establish a meaningful dialogue with the system. 
This addition would involve designing a special 
question-answer type of language. 
The aggregate production planning system can be 
extended to develop a total corporate planning system. 
This would involve, developing similar planning systems 
for all the divisions of an industrial organization such 
as marketing, production, finance, engineering and develop- 
ment and control (purchasing, accounting, etc.). Each 
System can be designed and tested separately to develop 


a fully integrated total corporate planning system. 
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APPENDIX 1] 


GLOSSARY OF THE VARIABLES USED IN 
THE COMPUTER PROGRAM 


The glossary is presented, herein, in two parts. 
The first part covers the main routine and the subroutines 
FCT), LEABFFE, OVEFF and WCOST; the second part, the sub- 


routine PATS only. 


A.1 Glossary for the Main Routine and Subroutines 
FCT], LEAREF ,2OVEPFVandawe0st 


A(T) a variable which is used for different purposes 
in different locations. in subroutine FCT! 1t 
represents the amount of penalty that will be 


u penalty cost con- 


charged each time the i 
Straint 1s Violated, 

In subrouctne EEAETT at=1s Used in: tne 
computations for determining the efficiency of 
the new employees during the training period 
and appears under the name "A" in the text 
(refer to Section 3.1 for further details). 

AA(J) a variable which is used in the subroutine 
LEAEFF in computing the efficiency of the new 


employees during the training period; appears 


under the name "a" in the text (refer to 
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Set tor 3B) 

ALR a variable in subroutine WCOST which represents 
the average (gross) labour rate for the 
employees. 

AMHDMX(K,Y) the maximum production capacity of the plant 
that will be available through (K) number of 
Silt Us during the yen year. 

B(J) a variable which is similar to A(J); appears 
under the name "B" in the text (refer to 
Selet won es) 

BB(J) Same as above, appears under the name "b" in 
theerexp tGréefier ito Section os 41 )e 

Bi) a variable name used for the variables in the 
data set - type 3, which represents the value 


h variable for the yeh year. 


of the L° 
Gis.) a variable name used for the variables (6) to 


(15) in the data set - type 2, represents the 
th th 


value of the L variable for the J product 
line. 

CASE(NP) the administrative expenses during the period/ 
month (NP). 

CASEK a variable in subroutine FCT! which means the 
same thing as CASE(NP). 

CCP(NP) the? costs tof luctuations’ in the production 
levels during the period (NP). 

CCPK a variable in subroutine FCTI, which means 


the same thing as CCP(NP). 
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CDM(NP) 


CDMK 


CDMT 


CDW(NP) 
CDWH 


CDWK 


CET(M) 


CETD(NP) 


CETDK 


CETI(NP) 


CERLK 


CFOH(Y) 


CHT(M) 
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the total cost of direct materials during the 
period (NP). 

a fie aie in subroutine FCT1l which is the 

same as CDM(NP). 

a variable, which represents the direct material 
costs for any particular product line during 

the period (NP). This variable is used in com- 
puting the value of CDMK and OCDM. 

the direct labour costs during the period (NP). 
the direct work force costs per production man- 
hours; used in the subroutine FCT]. 

a variable in subroutine FCT] which means the 
same things as CDW(NP). 

a variable in subroutine WCOST which represents 
the employee termination costs during the period. 


(M) of the npc 


planning horizon. 

the direct work force employee termination 
costs during the period (NP). 

a variable in subroutine FCT1 which means the 
same thing as CETD(NP). 

same as CETD(NP), for the indirect work force. 
same tasaCETDK} foritherindinectrwork fiorces 
the fixed portion of the overhead costs during 
the year (Y). 

the employee hiring and training costs during 


th 


the period (M) of the NP planning horizon; 


used in the subroutine WCOST. 
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the hiring and training cost for each employee 
hired; used in the subroutine WCOST. 

the expected cost of hiring and training an 
employee in the direct work force. 

same as above, for the employees from category 
(a) (refer to Section 3.8). 

same as above, for the category (b) (refer to 
Section 3.8). 

a variable which stores the values of CHTD(NP) 
in the subroutine FCTI and then transfers them 
to the main routine. 

same as CHTD(NP), for the indirect work force. 
same as CHIDA, for the indirect work force. 
same as CHTDAI, for the indirect work force. 
same as CHTDAI], for the indirect work force. 
same as CHTDK, for the indirect work force. 

a variable in subroutine FCT1l which is used in 
computing the inventory carrying costs. 

the standby costs of idle facility during any 
particular period; used in the subroutine 
EGind. 

the LonaTnades of indirect materials during 
the period (NP); used in the subroutine FCT]. 
a variable in the main routine which represents 
the inventory carrying costs during the period 
(NP). 


a variable in subroutine FCT] which means the 
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CINS(NP) 


CINSK 


CINSA(J) 


CIS 


-CIWK 


CLDO 


CLIO 
CLP(M) 
CLPD(NP) 
CLPDK 
CLPI(NP) 
CPPEK 


CME(NP) 


CMEK 


Per 


Same as above. 
the inventory shortage costs during the period 
(NP); geen in the main routine. 

a variable in subroutine FCT1 which means the 


same thing as CINS(NP). 


the cost of inventory shortages per unit (short) 


of therproduct) line: (J). 
a variable in subroutine FCT1 which is used in 
computing the value of CINSK. 

Same as CDWK, for the indirect work force. 
represents the direct labour costs, used in 
computing the manufacturing cost per unit for 
the period prior to the start of the plan. 
Same as above, for the indirect work force. 
the regular payroll costs of the employees 


during the wth ae 


period of the NP planning 
horizon; used in the subroutine WCOST. 

the regular payroll costs of the direct work 
force during the period (NP). 

a variable in subroutine FCT] which means the 
same thing as CLPD(NP). 

same as CLPD(NP), for the indirect work force. 
samemasvebROketiorethe indinect workstorces 
the marketing and distribution costs during 
the period (NP). 


a variable in subroutine FCT] which means the 


same as CME(NP). 
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COH(NP) 


COHK 


COT(M) 


COTD(NP) 


COTDK 


COTI (NP) 


COTIK 


COVI(J) 


CPRK 


CPTH 


Chis 


CPUK(JC) 


CSPD(NP) 
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the overhead costs during the period (NP). 
Same as above, used in the subroutine FCTl. 


the cost of employee overtime during the yeh 


period of the np’? 


planning horizon; used in 
the subroutine WCOST. 

the direct work force overtime costs during 
the period (NP). 

a variable which stores the values of COTD(NP) 
in the subroutine FCT1 and then transfers them 
to the main routine. 

same as COTD(NP), for the indirect work force. 
same as COTDK, for the indirect work force. 
the value of the beginning inventory stock, 

of the product line (J), at the start of the 


month/period one. 


a variable in subroutine FCT] which is used 


in computing the value of the objective function. 


the cost of production time per man-hour of 
the total work force (see Section 3.3 for a 
more complete explanation). 

the total production costs during the period 
NPy used imethe subroutine FCTT. 

the manufacturing cost per unit of the product 
linee(JG), during therpenioda( NP); sused «in 

the subroutine FCTl. 

the shift premium costs of the direct work 


force during the period (NP); used in the main 
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CSPI(NP) 
CSPDK 


CSPIK 
CSPO 


CSSA(T) 


CSSK 


CUMC 


CVIT(NP,J) 


CVITK(JC) 


CW] 


DHP 


DMHD 
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routine. 

Same as above, for the indirect work force. 

a variable in subroutine FCT1 which means the 
same thing as CSPD(NP). 

Same as above, for the indirect work force. 
the shift premium costs for the period prior 
to the start of the plan. 

the coefficient for the shift start-up and 
Shut-down cost equation which represents the 
value for the pth year. 

the shift start-up and shut-down costs; used 
in the subroutine FCT]. 

a localized variable in the main routine which 
is used to cumulate the cost of operations for 
the whole plan. 

the value of the ending inventories for the 
product line (J) during the period (NP). 

a variable in subroutine FCT1 which means the 
same thing as CVIT(NP,J). 

a variable in the subroutine OVEFF which is 
used in computing the overall efficiency of 
the work force. 

the amount by which the constraint on the 
(percent) level of hirings in any particular 
period is exceeded; used in the subroutine 
WCOST. 


a variable in subroutine FCTI which is used 
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for computing the value of PMHD (refer to 
Section 3644). 


DOP same as DHP, for overtime. 

DTP Same as above, for the employee termination. 
DUP Same as above, for the employee idle time. 
DWD a variable in subroutine FCT] that measures 


the (absolute) amount by which the direct 
work force level changes during the period (NP). 
DWI same as above, for the indirect work force. 
E(L) a variable which has different meanings in 
different locations. In main routine and sub- 
routine FCT] it represents the value of the 


, th 


variable of the data set - type 4. 

In subroutine OVEFF this variable means 
the average efficiency of a new employee dur- 
ing each period of the training program. 

Edel) a variable in subroutine LEAEFF which is used 
in computing the average efficiency of a new 
employee during the learning period. 

EDO(NP) the average efficiency of the direct work 
force employees during the period (NP). 

EDOK a variable which transfers the value of EDO(NP) 
from the subroutine FCT1 to the main routine. 

EH(J) the initial efficiency of the new employees, 
from the category (J) (refer to Section 3.6 
and 3.8), at the start of the training period; 


used in the subroutine-LEAEFF. 
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same as above, for the category (a). 
Same as above, for the category (b). 


the average efficiency of the new employees 


‘during the (entire) training period; used in 


the subroutine LEAEFF. 

same as ibove, for the category (a). 
same as above, for the category (b). 

the expected average efficiency of a new 


employee during the poh 


period of the train- 
ing program. 

the average efficiency of a new employee, 
hired from the category (a), during the train- 
ing period (see Sections 3.6 and 3.8 for a 
more complete explanation). 

Same as above, for the category (b). 

same as EHD(I), for the indirect work force. 
same as EHDA], for the indirect work force. 
same as EHDA2, for the indirect work force. 
the initial efficiency of the work force at 
the beginning of the period/month one. 

same as EDO(NP), for the indirect work force. 
a variable which transfers the values of 
EIO(NP) from the subroutine FCTl to the main 
routine. 

the average efficiency of the work force dur- 


h 


ing the me period of the planning horizon; 


used in the subroutines OVEFF and WCOST. 
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FMDK(J) 


FMHW(M) 


FMHD(M) 
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the expected average efficiency of a new 


employee during the yth 


period of the training 
program; used in the subroutine LEAEFF. 
a variable which represents the value of the 


, th 


variable of the data set - type 5. 

the capacity utilization factor; used in the 
subroutine FCTl. 

the capacity utilization factor during the 
period (NP); used in the main routine. 

a variable, which transfers the values of FCU 
from the subroutine FCT1 to the main routine. 
the man-hour equivalents of the production 
level of the product line (J) during the period 
(M) of the planning horizon; used in the sub- 
routine FCTI. 

a variable which stores the values of FMD(1,J) 
for future use, during the optimization of the 
next period; used in the subroutine FCTI. 

the work force requirements, in man-hours, 


during the ne} 


period of the planning horizon; 
used in the subroutine WCOQST. 

the man-hour equivalents of the total direct 
work force requirements during the period (M) 
of the planning horizon; used in the subroutine 
FG Tale. 


a variable in subroutine FCT! which stores the 


value of FMHD(1) for future use, during the 
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optimization of the next period. 

same as FMHD(M), for the indirect work force. 

a real variable in subroutine LEAEFF which is 
equivalent to NHP. 

a variable in the main routine which is used in 
computing the coefficient for the shift start- 
up and shut-down costs. 

a variable in subroutine LEAEFF which repre- 
sents the length of the training period, in 
hours, for the employees hired from the cate- 
gory (J). This variable appears under the name 
"NHT" in the text (refer to Sections 3.6 and 
3.8 for a more complete explanation). 

Same as above, for the category (a). 

Same as above, for the category (b). 

a variable in subroutine OVEFF which is used 

to compute the number of employees hired during 
the period (K). 

a highly localized variable which is used for 
different purposes in different locations, 
mainly for indexing in the program-loops. 

a variable which represents the maximum dura- 
tion of training programs in periods (rounded 
to the next whole number) for the new employees 
from each category. This variable is used in 
the subroutines LEAEFF and OVEFF. 


a counter that records the number of times 
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(in a planning horizon) the constraint for 

the upper limits on the levels of month-ending 
inventories has been violated; used in the 
subroutine FCT] as a penalty cost factor. 

same as above, for the lower limit on the 

level of inventory stockss® 

ar pena lityercostsfactor for the constraint on the 
level of inventory shortage which is initiated 
at a zero value and incremented, each time the 
constraint is violated, by a value equal to 

the amount by which the constraint has been 
exceeded. 

a multiplying factor to charge the penalty 

cost for violating the constraint on the 
(direct work force) overtime. The value of 

I1D is computed by summing up the percent levels 
by which the constraint is violated during each 
period of the planning horizon. 

same as above, for the employee idle time. 

same as above, for the rate of hiring. 

Same as above, for the rate of employee termina- 
tion. 

same as 11D, for the indirect work force. 

same as IIF, for the indirect work force. 

same as I1G, for the indirect work force. 

same as IIH, for the indirect work force. 


a highly localized variable which is used for 
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different purposes in different locations, 
mainly for indexing in the program loops. 

JA a highly localized variable which represents 
the data set type-number in the data reading 
routines. | 

JB same as above, for the item/variable-number 
within a data set. 

JC same as above, for product lines. 

JD same as above, for the year (number) corres- 
ponding to the particular data value. 

JF Same as above, for the total number of vari- 


ables in a data set. 


Jd the total number of product lines being manu- 
factured. 
JPROPX an arbitrary factor used to penalize inventory 


deviations from a specified minimum inventory 
level (refer to Section 2.2.1, the master 
control data set, for a more complete explana- 
tion) 

K a highly localized variable, used mainly for 
representing the month/period in the program 
loops. 

KFINAL a variable, when the value is equal to one 
computes the required planning data based on 
the results of the search (for optimal decisions). 

KJ the number of independent variables in each 


month/period (KJ = JJ + 2). 
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a highly localized variable, used as an index 
in some of the DO loops. 

Same as above. 

the number of periods that will be considered 
during optimization in the SDR decision model 
(KM = NPY x NOY + NHM - 1). 

a variable which, when set = 1, Suppresses a 
large fraction of the detailed output. 

a variable in subroutine WCOST and OVEFF which 
is used for indexing the DO loops. 

a variable, which is used to represent the 
period/month of the plan in various program 
loops. 

same as KQ. 

a highly localized variable, used in different 
locations for different purposes. 

same as III; used in the main routine and the 
subroutine FCT] for the direct work force. 
same as above, for the indirect work force. 
see the glossary for the subroutine PATS. 

a variable in subroutine WCOST which is used 
as an index in the program. 

a localized variable, which is used to repre- 
sent the period/month of the planning horizon 
in various program loops. 

the number of shifts the plant operates during 


the period (NP). 
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a variable which transfers the values of MS(NP) 
from the subroutine FCT] to the main routine. 
the number of shifts the plant operated during 
the period prior to the start of the plan. 

the maximum number of shifts the plant can 
operate. 

a (local) variable used in a DO loop. 

the number of independent variables in the 
search (N«=«NH x-Kd), 

see glossary for the subroutine PATS. 

the length of the planning horizon in months/ 
periods. | 

the number of regular hours of work per period/ 
month. 

the length of the entire planning span in 
periods/months. 

the length of the entire planning span in years. 
a counter which represents the period/stage 
being optimized. 

the number of periods/months per year. 

a variable used in subroutine OVEFF as a 
counter. 

a controlling variable whose value is either 
one or zero, which branches the choice of 
approach for estimating the indirect work force 
requirements in the model. Regression method 


is used when the value is one. Extrapolation 
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is used when the value is zero. 

a highly localized variable used in a DO loop. 
the direct material costs during the period 
prior to the start of the plan. 

Same as above, for the indirect materials. 

the manufacturing cost per unit of the product 
line (JC) during the period prior to the start 
of the plan. 

represents the maximum value among the values 
Of VOCPUCIG) WRIGH=CPP293 HTS. FIG 

a variable used in computing OCPU(JC). 
a variable in subroutine OVEFF which represents 
the average efficiency of the work force at 
the beginning of the period one. 

Same as above, for the direct work force; used 
in the main routine and the subroutine FCT]. 
Same as above, for the direct work force; used 
in the main routine and the subroutine FCT]. 
the man-hour equivalents of the production 
level for the product line (JC) during the 
period Upriorsto*the start *of the plan: 

Same as above, for the total production of 

all the product lines. 

the number of units of product line (JC) pro- 
duced during the period prior to the start of 
the plan. 


a variable in subroutine OVEFF which means the 
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same thing as OW. 

the level of the direct work force overtime 
(expressed as a fraction of the work force 
level) used during the period (NP). 

Same as above, for the indirect work force. 
the percent level of the direct work force 
overtime used during the period prior to the 
start of the plan. 

Same aS above, for the indirect work force. 

a variable which transfers the values of 
OTDP(NP) from the subroutine FCT1 to the main 
routine. 

a variable in subroutine WCOST which is used 
in computing the amount of work force overtime 
and idle time. 

same as OTDPK, for the indirect work. force. 
the employee overtime (expressed as a percent 


of the work force level) during the wth 


period 
of the planning horizon; used in the sub- 
routine WCOST, 

the beginning inventory level of the product 
line (JC) at the start of the period one. 

the work force level at the start of the 
period/month one; used in the subroutines 
WCOST and OVEFF. 


Same as above for the direct work force; used 


in the main routine and the subroutine FCTl. 
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Same as above, for the indirect work force; 
used in the main routine and the subroutine 
ECTS 

a variable in subroutine LEAEFF which repre- 
sents the expected proportion of the employees 
hired from the category (J) (refer to Section 
3.8 for a more complete explanation). 

same as above, for the category (a). 

same as above, for the category (b). 

a variable in subroutine WCOST which represents 
the value of the constraint (upper limit) on 
the rate of hiring during any particular 
period. The rate of hiring is defined as the 
ratio of the new employees hired to the start- 
ing work force level at the beginning of the 
period. | 

a variable in subroutine WCOST which is used 
in computing the employee termination costs 
(refer to Section 3.7 for a more complete 
explanation). 

a variable in subroutine FCT1 which is used 
in computing the cost of fluctuations in pro- 
duction levels, appears under the name “PMHD," 
inmethey text? (refer to Section 3.14). 

same as PHTU, for the overtime. 

a vector of penalty factors used in applying 


the JPROPX concept (refer to Section 2.2.1, 
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master control variables for an explanation 

of JPROPX). 

a localized variable, used in computing 
PROP(M). 

a localized variable in subroutine FCT1, used 
in applying the PROP(M) function. 

the maximum acceptable level of employee turn- 
over; used in the subroutine WCOST in computing 
PLM (refer to Section 3.7). 

the actual average level of the employee turn- 
over; used in the subroutine WCOST in computing 
PLM. 

Same as PHTU, for the employee terminations. 
Same as above, for the employee idle time. 

the production level for the product line (J) 
during the period (K). 

the expected value of the learning rate con- 
stant sbnede met oeSectiron: .3..4) )s. 

the overtime (wage) rate premium. 

a variable in subroutine OVEFF which means the 
same thing as W(K). 

the sales forecast for the product line (J) 
for the period (K). 

a dummy variable. 

the value of the objective function before 
adding. the, penadLtysicos:ts. tosiut.. 
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penalty costs, see the glossary for the sub- 
routine PATS for more details). 
T(M,J) the sales forecast for the product line (J) 


for the wth h 


period of the npt planning hori- 
zon. This variable means the same thing as 
SK S09" TTY + NPRUIF= SCKAI) 1. 

TCOP(NP) the total cost of operations during the 
period (NP). 

TCOPK a variable in subroutine FCT1 which means the 
same thing as TCOP(NP). 

TCPK the value of the objective function. 

Tepe the total variable production costs during the 
entire planning horizon. 

TNK(J) a variable in subroutine LEAEFF which is used 
in computing E(U5K)" 

TNP (J) same as above, 

UCP a variable in subroutine FCTI which means the 
same thing as CPUK(JC) and is used in computing 
the inventory carrying costs. 

UMHD (JD) a localized variable in the main routine which 
is used in computing the range for the direct 
work force decision vectors. 

UTDP(M) the level of the direct work force idle time 
(expressed as a fraction of the work force 
level) during the period M of the planning 


norrzon: USedvin tne’ subDrTOoULINe Port. 


UTIP(M) same as above, for the indirect work force. 
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the amount of employee idle time in man-hours 
during any particular period. 

the level of the employee idle time expressed 
as a fraction of the work force level, 

the ending inventory level of the product line 
(J) during the period (NP). 

a variable which stores the values of VI(1,J) 
in the subroutine FCT] and then transfers them 
to the main routine. 

the lower limit on the levels of inventory 
stocks of the product line (JC) during the 
year GY). 

the inventory level of the product line (J) 

at the start of the period one. 

the level of inventory shortage (expressed as 
a fraction of sales level) of the product Tine 
(J) during the period (NP). 

a variable in subroutine FCT1 which means the 
same thing as VISP(NP,J). 

the upper limit on the level of inventory 
stocks of the product line (J) during the 

year (Y). 

a variable in subroutine FCT1 which computes 
the number of inet abe which the inventory 
level of a product line during any period is 
below the minimum reserve stock level. 


the work force level during the period (K); 
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used in the subroutines OVEFF and WCOST. 

WD (NP) the direct and the indirect work force levels 

Bes during the period (NP). 

WDMAX the upper limit on the levels of the direct 
work force during any period. 

WDO a variable in subroutine FCT1; used in comput- 
ing the cost of fluctuations in the production 
levels. 

WH a localized variable in subroutine WCOST which 
computes the number of employees hired during 
any particular period. 

WO Same as above. 

WHP the level of the employees hired, expressed as 
a fraction of the work force level; used in 
the subroutine WCOST. 

WIF the indirect work force requirements (number 
of employees) during any particular period; 
used in the subroutine FCTI. 

WIFI | the number of inspection and quality control 
employees required during any particular 


period; used in the subroutine FCT1. 


WIMAX same as WDMAX, for the indirect work force. 
WIO0 same as WDO, for the indirect work force. 
WT the number of employees terminated during any 


particular period; used in the subroutine WCOST. 
WTP(LS) same as WHP, for employee terminations. 


X(K) the kth independent variable (see glossary for 
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the subroutine PATS). 

the upper and lower bounds on the acceptable 
vValuescofex(kKi. 

a variable in main routine, used for printing 
the output (final results). 

same as above. 

Same as above. 

Same as above. 

variables in main routine, used for printing 


out the final results. 


A.2 Glossary for the Subroutine PATS 


The material in this section is taken from reference 


[40]. 


ALP 


Dien 


D(1) 


alpha, the factor by which the step size, 
D(I), grows when a forward move is successful 
andeca = "2. "tntenattzed ataas, used in 

stm 292). 

beta, the multiplicative factor by which step 
size for an independent variable is reduced 
if forward and reverse move for that variable 
Taal 

the current value of the step size for the 
yeh independent variable. 

delta, the multiplier which is used to deter- 


mine the initral watue of Dil), the step size, 


in accordance with statement 180. 
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a local quantity used to determine whether the 
lower bound on the step size has been reached. 
DX is computed and used only between stms 480 
and 485, where step size reduction takes place. 
a quantity used to increment the value of GR 
in the adaptive logic. Dl is set in stms 

$02 (+ elipandaso4y TteTrstusedetin'’s tmss10; 

the factor by which the pattern move vector 

is multiplied to obtain the actual size of the 
pattern move. (GR is initialized at 2.2 and 
adjusted upward, usually by increments Ofte), 
in statements 510 through 783.) When GR 
meaehest3e5 tithis vesatrtoe2 27 

a highly localized variable used in an index 
in DO loops - see stms 180, 420, 786. 

counter to record passes through stm 802. 
counter to record passes through stm 803. 
counter to record passes through stm 804. 
counter to record passes through stm 801. 

a printout control character of little import- 
ance; initialized at 1 and left at that value. 
It is tested in stm 888; if the value is > ] 
it causes deletion of certain print lines? 

(In effect, it is not used unless initialized 
AtehERS> i) 

a subscript which defines which independent 


variable is now being studied -- as in X(K), 
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DOK) § sete. 

a counter for one plus the number of variables 
studied since the last test for a new base 
point. When KK reaches N + 1 atest for a 
new base point is made. (The foregoing sent- 
ence applies when the subroutine is in the 
full exploratory search mode, as signified 

by LT < 0 -- this is the usual mode.) KK is 
Set equals CO ml at SUMS: 1OU e994 40 “toc geatnG 
784. In the full exploratory search mode 

(LT < 0) KK is incremented at 330 and tested 
at 330 + 1 and 404; in the truncated search 
mode (LT > 1) KK is incremented at 778 and 
tested at 404. 

counter of the number of times we enter the 
adaptive logic preparatory to attempting a 
pattern move. (--has same value as LT7). 
Incremented at stm 510 + 1. 


ae printout contro chardadcver mn ich 1s) a nput 


on the first data card read by the main routine. 


KPRESS = 1 suppresses about 99% of the output 
generated by subroutine PATS: KPRESS = 0 
allows the full details to be output during 
PATS. 

a Master Monitor of Subroutine Status which 
tells where to GO TO next. The primary job 


of LA is to control traffic through the Bound- 
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ary Check and the objective function subroutine. 
The following values of LA correspond to the 


following destinations: 


Destination 


Value Stm 
of LA No. Task 
] 100 Initia lizetion 
Z eo2 Forward Exploratory Move - Normal 
3 463 Reverse Exploratory Move - Normal 
(Following Fwd Failure) 
4 580 Evaluate Base Point Following 
Pattern Move 
5 285 Forward Exploratory Move - Follow- 
ing Pattern Move 
6 466 Reverse Exploratory Move - Follow- 
ing Pattern Move (and Fwd Failure) 
7 510 Attempt Adaptive Pattern Move 
8 500 Terminate Search and Exit from 


Subroutine. 


an approximate upper. bound on the number of 
calls of the objective function during a 
particular call of subroutine PATS. (LIM is 
usually set at 3000, but a different value 
might be appropriate for some applications.) 
See final remark under NEVAL for identification 
of location where NEVAL is tested. 


set at 0 in fifth statement of subroutine and 
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kept there -- never referred to again except 
in statement 784 + 1 (see Weisman, Wood, and 
Rivlin for explanation). 

controls the choice between the standard full 
exploratory search mode (LT < 0) and the 
truncated search mode (LT > 1). Under trunc- 
ated mode the exploratory search is stopped 
as soon as any move produces an improvement 
in the last base point. This point is saved 
as the new base point, and a pattern move is 
made. The next exploratory search starts with 
the variable after the one which produced the 
last success. Under the full exploratory 
search mode, which is the mode this coding 
reflects, since LT = 0, an exploratory step 
is made with all N variables before a pattern 
move is attempted. LT is set in the fourth 
statement of the subroutine and does not appear 
again except in statements 320 and 400. 
number of times this subroutine has reached 
statusabA tsi2%. 

number of times this subroutine has reached 
Status sLAT=a3 40 

number of times this subroutine has reached 
status LA = 4. 

number of times this subroutine has reached 


status LA = 5*., 
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number of times this subroutine has reached 
statusslAreroe . 
number of times this subroutine has reached 
status LA = 7. 

*jn the event that a proposed forward 
or reverse move fails the boundary check 
the appropriate one of these counters 
will fail to be incremented. 

a status variable always equal to 1 or 2: 

L4 = 1 says we are in the process of making 
exploratory moves in normal fashion, i.e., we 
are searching for a pattern. L4 = 2 Says we 
are in the process of making exploratory moves 
following a pattern move. L4 is set at stms 
1802+) 1081190 ,1580)+) 1seandh Gineorrectt ys. to 
no effect) at 492. It is tested at stms 290, 
200, e@e60teS508 77942Zand- 492. ,At the time “its 
vValuevds loutputeby stm? SotplgeLe4t=vdrmeans 
"the pattern move we're in the process of 
attempting is the first attempt in an effort 
to establish a new series of pattern moves." 
L4 = 2 means "the pattern move we're in the 
process of attempting has been preceeded by 
one or more successful pattern moves in this 
series." 

one plus the number of independent variables 


which have experienced both forward and reverse 
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failures with the minimum step size since the 
last test for a new base point. The search is 
terminated when Ml] is > N + 1 during the 
exploratory search mode, i.e., when L4 = 1] 

and we are attempting to establish a new 
pattern. Even though M1 may be (redundantly) 
manipulated while the search is conducting a 
post-pattern-move exploration (when L4 = 2), 
it is never used, i.e., tested, except when 

we are attempting to establish a new pattern, 
namely when L4 = 1. (In reprogramming, this 
variable could be reduced by 1, i.e., initia- 
lizedeto 0 andttested for -> pNsjiwith improved 
clarity of interpretation; this also applies 
to M2s) Ml is initialized in stms 180 + 6, 
300 tu2dse3s52crahs and r440 ct ky incremented rat 
490; and tested at 429 + 1. 

one plus the number of independent variables 
which have experienced both forward and reverse 
failures since the last attempt to make a 
pattern move. The pattern is considered broken 
and the search is restored to exploratory mode 
(L4's value changes from 2 to 1) if all vari- 
ables fail following a single pattern move 
attempt, Ihe) Wood fall sartielesand ttheur 
source code says (--page 63, under M2) that 


the search is restored to exploratory mode 
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following five successive attempts to make a 
pattern move, and NPF counts these five fail- 
ures. But there's a bug in their code (and 
Taubert's) which causes (a) NPF never to get 
as high as 2, and (b) the only test of NPF 
(stm 353) never to be reached unless a new 
base point has been established after a 
successful pattern move. Such a successful 
move, is extremely unlikely when all adjust- 
ments to the attempted pattern move have 
failed. The net result is that the pattern 
is broken after the failure of a single attempt 
to make a pattern move -- probably a good 
nesubtetf athetbesthvaldesotf Ssinsstmiss3nts 
too high. In this area the source code could 
be condensed and clarified to an advantage, 
however, by modifications between stms 340 
and 353, and elimination or constructive use 
of NPF as a variable. M2 is (redundantly) 
manipulated when we are in the exploratory 
search mode !(L4e=°1), ebutSit is smever sincre= 
mented or tested except when we are in the 
exploration-after-pattern-move-attempt mode: 
Be si2. IM24hs THN ha huved aingstns 880" 7), 
300 + 1, 440, and 780; incremented at stm 493; 
tested at 493 + 1. 
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dimensionality of the space being searched. 

a counter of the number of evaluations of 
FCT1, the objective function, which have been 
made during this particular call of subroutine 
PATS. NEVAL is incremented inside subroutine 
FCT1; it is not tested every time it is incre- 
mented. LIM is an. approximate upper bound on 
NEVAL; the test is made in stm 7, just prior 
to attempting a pattern move. This is prefer- 
able to testing every time FCT1 is called. 
variable used to store the previous value of 
NEVAL, the "total number of evaluations of the 
objective function made thus far." NEVOLD 
Tssresetoatistm? Sl5a+ 2vandeas used inethe 
computation? ofVV’ inastmesos, 

counter for the number of successive pattern 
move (attempts) which are followed by failure 
of all individual steps that try to adjust the 
pattern move attempt. (Set at 0 in stms 

TO0U*s 5h 300 +76, 7 4208F 45 Oihicrementedt at 
Stmss780 +) 13 tested’ at-353..)" “See discussion 
under M2. 

a variable used to store a certain prior value 
OTPSNOSTICOTSEInitializedmatmetmei so: +e si 
reset at stm 815, and used in the computation 
of V at stm 398 (within the adaptive logic). 


variable used to store the previous value of V 
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Gresetringsstm 815 +.1; used in stm 782 to 
compare old and new values of V). 

a highly localized temporary storage variable 
used only between stms 530 + 3 and 530 + 5, 
where the pattern move is attempted. 

ass torage (matrix. for esitoring ithe “old, si vel; 
base point, values of the independent vari- 
ables. These values are initially set equal 
to the X(I) values and are updated in stm 

530 + 4. If an attempted pattern move fails, 
the old values of X(I) are recovered from 

Q(I) in stm 420. 

a variable used to store the value of the 
objective function at the most recent base 
point. SSCovseinitiabnzedGatwstm 180 th4, yre= 
set at stm 530 (just prior to attempting a 
pattern move), tested at statement 340 to 
determine whether a new base point has been 
established, and used to restore SP to the old 
value in stm 410 if the attempted pattern move 
fails. 

the value of the objective function returned by 
the#s ubrnoutrnel hes liei.es3 ethescunviati oneof 
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a redundant variable which is initialized at 


stm 180 + 2 and never referred to again for 
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any purpose. 

a variable which (except in the following cir- 
cumstance) is equal to the minimum value 
returned by the objective function thus far. 
The exception occurs when an exploratory search 
is being conducted immediately after a pattern 
move.o” Ine that: situation SP’ 1s set (stm 580+ 
2) equal to the value of the objective function 
which reflects the unadjusted pattern move, 
and thereafter, as the exploratory search to 
adjust the pattern move proceeds, SP is up- 
dated to reflect the minimum of this value and 
the best exploratory move to date. After the 
explorations are completed, this updated value 
of SP is compared (stm 340) with SC, the value 
of the objective function at the last base 
point. If this base point test is passed the 
(adjusted) pattern move is declared success- 
ful and a new base point is established. If 
the base point test is not passed the pattern 
move attempt is declared unsuccessful, SP 

is restored to its old value (stm 410), and 
local explorations are initiated about the 

old base point in an effort to establish a 

new pattern. (With reference to the source 
eoqe. SP appears in stms 180 + 5, 280, 300, 
3402-741 0-460 785-9530, ange oe0nt 2.) 
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a quantity used (stm 480 + 2) in obtaining 
lower bounds on step sizes for exploratory 
moves. (has been set at 1 x Te perhaps 
should be changed.) See also stm 340, where 
TOL is not used but perhaps should be. 
percentage improvement in value of objective 
function per call of objective function, Com- 
puted and used in adaptive logic - see stms 
398 through 782. 

the independent decision variables (of which 
there are N: 1 < I < N). In the aggregate 
production planning model the vector X(I) 
contains the decisions for each of the (NH) 
months in the planning horizon. In each month 
the first value is the direct work force size, 
the second value the indirect work force size 
and the subsequent values are the production 
levels for each product line. 

the upper and lower bounds on the acceptable 
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APPENDIX 2 


ANNOTATED LISTING OF THE COMPUTER PROGRAM, 
WITH SAMPLE DATA 


This appendix contains a listing of the computer 
program developed for the aggregate production planning 
system. The program, written in FORTRAN IV, was run on 
an Amdahl 470V/6 machine at the University of Alberta 
Computing Services installation. A listing of the sample 
data, used for testing the model, is given at the end of 


this appendix. 
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FORTRAN CODE FOR THE AGGKEGATE PRODUCTION PLANNING MODEL, 
SDR VERSION.» PKOGSAMMED IN MAY,1976 BY TEJINDER SINGH. 
THE SEARCH ROUTINE ( SUBROUTINE PATS, BELGW) USED IN THIS 
PROGKAM WAS TAKEN FROM KEFE RENCE (40). 

—--SEE GLOSSARY FOR THE DEFINITICN OP THE YARIABLES 
USED iN THE PRCGEBAMS, BELOW. 


THE MAIN ROUTINE BEGINS HERE 


THIS CARD CALLS FOR DOUBLE PRECISION REAL VARIASLES 


qaHaaeaanga aagaaga aes 


IMPLICIT REAL*8 ( A-H,O-Z) 


oO 


DIMENSTON S(930,10),28(50,10),V1(80,10),WD( 80), 
1VISP(S0,10),CCP( 80 ),CDM( 80 ),CDW{ 80), 

2CME( 80 ),CASE{ 80), COH( 50 ),CETDI 80), W171 80 ) 
3,CETI( 80 ),CHID( 80 ),CHTI( 80 ),CLPD(80),CLPiI( S80), 
4CINC(80),CiNS( 80 ),COTD( 80 ),COTI( 50 ),CSPD( 80), 
SCSPi(80),CViIT( 80,10), EDG( 80 ), E101 80 ), OTDP( 80), 
6OTIP( 30 ),FCUP( 80 9), TCOP( 80 ),0P110),2Z1016), 

7Z2( 6 )y9Z31 G6) yZ 40 6) 251 6) ZE61 6 IZ 716 D4 Z81 6 )D,Z9(6) 
8,Z116 ),COVI( 10), UMHD( 6 ),MS( 80) 


COMMON 1(412,10 ),8(18,6),C(13,10),8165),F( 15), 
LAMHDMX(4,6),X(150),XMAX({ 150), XMIN(150),0V1(10), 
2ViL( 10,6 ),V1U( 10,6), ERD( 30),EHI( 50), 
SCINSA( 10 ),CSSA(6 ),PROP{ 12), VISPK(10), 
4V1iK(10),5N,CFMDi 10), CFPOH(6),6CPU(10),CVATK( 10), 
SCCPK,CDMK, CDWK,yCMEK,COHK,CASER,CETDER ,CETI£KR.,CPTH, 
6OCHTDK ,CHUTI »,CLPDK,CLPIK, COTDK,COTIK, CINCK,CINSK, 
7TCSPDK,CSPIK, TCOPK,EDOK, EIOK ,OTDPEK,OTIPK ,OCPUM, 
SFCUPK,OFMHD ,OWD,OWi,OED,CEI,CHTDA,CHTIA,SJ,SIGMNA, 
t+MX yNEVAL,KFAINAL,KPRESS 
c 
C READ IN THE MASTER CONTROL CARD 
READ{( 5,306 )ACY,NPY¥,NH yNHM,JJ,NHP,JPROPK,NWI,MX,KPRESS 
WRITE( 6, 306 JNOY,NPY,NH,NHNM,JJ,NHP,JPROPXK,NWI,yMX,KRPRESS 
1,Lim 
NN=NOY*¥NPY¥ 
KM=NAtNH 
KN=NN+NHW-1 
Cc 
C READ IN THE DATA SET —- TYPE 1 


1A50G0e OF DMAALS HORT: 
-HOula S9080L9T Va 
iar #) @220 ( woeee 


PaGaL RAY SE 
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c 
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READ( 5,303)dA,JF 

WRITE (6,303) JA,JF 

READ (5,300) CWD,OW1,OCED,OEI,OTDPI,OTIPI 
READ( 5,301) MSO 

WRITE (6,500) OWD,OW1,OED,CEI,OTDPI,OTIPI 
WRITE (6,504) MSO 

JC=0 


C READ IN THE DATA SET -—- TYPE 2 


Cc 


10 JC=JCt1 
WRITE (6,502) JA,JF 
READ (5,301) oP( JC) 
READ( 5,301) CVI( JC) 
WRITE (6,501) OP(Jc) 
WRITE (6,501) OVI{JC) 
SALES FORECASTS 
READ( 5,300)(S(1,3C),i1=1,KN) 
UPPER AND LOWER LIMITS ON TRE LEVELS OF 
OF INVENTOCEIES FOR EACH YEAR 
TF (NOY.GT.2i1) GO TO 2 
READ (5,301) VIL(JC,1) 
READ (5,301) ViU(JC,i) 
WRITE (6,501) ViLiJC,1) 
WRITE (6,501) Vil(JC,1) 
Go To 3 
2 READ (5,900) (ViL( JC,1),1=1,NC0Y¥) 
READ (5,300) (V1U(JC,1),1=1,NOY) 
WRITE (6,500) (ViIL(JC,1),i1=1, NOY) 
WRITE (6,50C) (V1U(JC,1),i=1, NOY) 
COST BATA BWHICH IS DIFFERENT FOR GACH PRODUCT LINE 
3 JF=JF-5 
READ( 5,300) (C(J8B,JC),JB=1,J35F ) 
LP( JC.LT.~JJ) GO TO 10 


C READ IN THE DATA SET - TYPE 3 


READ(5,303) JA, JF 
WRITE (6,502) JAyJP 
DO 13 JH=1, JF 
IF (NOY.GT.1) GO To il 
READ (5,301) H(JB,1) 
WRITE (6,501) B(J35B,1) 
GO To 13 
11 READ (5,300) (8(JB,1),1=1,N0Y¥) 


1S5 


WRITE (6,500 ) (58(J6,1),1=1,N0Y) 
13 CONTINUE 


: © 

© READ IN THE DATA SET - TYPE 4 
READ( 5,303) JA,JF 
READ(S,300)( E(.J8 ),J5=1,JF) 
WRITE (6,500) (E(J6),JB=1, JF) 
B( 39 )=1-E( 39) 
£E( 40 )=1-EFi 40) 

c 


C READ IN THE DATA SET - TYPE 5 
READ( 5,303) JA,JF 
WRITE (6,502) JA,JF 
READ( 5,300)( F(JB),JB=1,JF) 
WRITE (6,500) (F(3B),3B=1,J5F) 


READ IN THE DATA SET - TYPE 6 (THE INITIAL STARTING 
SULUTION) 
READ (5,303) JA,JF 
WRITE (6,302) JA JSF 
if { NHM.GT.~1) GO To 1 
READ (5,901) WD(1) 
READ (5,301) Wii) 
WRITE (6,501) Wb i) 
WRITE( 6,501) WI(1) 
GO To 7 
i READ (5,300 }( ¥p(1),1=1, NHM) 
READ (5,900) (W111),i=1,NHMN) 
WRITE (6,500) @WD(1),i1=1,NHM) 
BRITE 16,50) (Witi),i=1,NHM) 
7 L=JJ+i1 
BO 5 J=1,1L 
if (JeGT.JJI) GO TG 125 
if (NHN.GT.~1) GO TO 4 
READ (5,301) R(1,J) 
WRITE (6,501) R{1i,J) 
GQ TO 5 
4 READ (5,300) {R( K,yJ),K=1,NHM) 
WRITE (6,500) (R(K,J5),K=-1,NHN) 
® CONTINUE 


OO G 


ROUTINES FCR CALCULATING THE COST COEFFICIENTS TO BE USED 
LIN THE OPTIMIZATION MODEL 


©oqg aa 


(ow, i=14¢ 1,0u pa) Lontene gue 


Oo go G 
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156 


KJ=55+2 
N=KJ* NH 


CALL IN THE SUBBOULTINE LEAEFF 20 DETERMINE THE { WORK ) 
GEPPRPICTIENCY OF THE NEW EMPLOYEES DURING TRAINING /LEARNING 
PERICDS 


CALL LEAEFF(E(20 ),E(16),H(17),6(14),E(15), 
16( 12),8(13),NHP,EHD,LDI,EHDAL, EHDA23 
CALL LEAEPP(E(35 ), E126 ),E(27),8(28),E(29), 


161 24),E 25),NHP,EHI,LIi,EHIAL, EHIA2) 


IF ( KPRESS-eNE.O0) GO TO 70 

WRITE (6,62) 

WRITE (6,602) EHDA1,EHDA2 ,(EHD(1),1=1,LD1 ) 
WRITE (6,621) 

WRITE (6,602) EHIA1,EHIA2 ,( EHI{ 1),i=1,L11) 


C FIXED PORTION OF THE GVERKRHEADL COSTS, CFOR(Y ) 
70 i=0 


ome 


GOOG 


2s 


Vi 


i=iti 
CPOH( 1 = B( 1,1 )+51 2, f +81 3,1 5+B( 4,1 )+ 


1B( 5,4 )*B( 6, 13480 7,1 94+8( S,fi NBC G,1 +B 10,1) 


IF (i1.LTeNGY) GO TO 235 

IF (KPRESS»NE.O) GO TO 71 

WRITE (6,601) 

WRITE (6,602) (CPOH( 1), i=1, NOY) 
OFMHD=0 


COST COEBRFPPICIENT FOR THE COST OF PRODUCTIVTY LOSSES, CPi 


CPTH =CFOHR( 1)/({ GWDtOWT J*FNHP*NPY ) 


MANUFACTURING COST PER UNIT FOR THE PERICD/MONTH PRIOK 
TO THE START OF THIS PLAN 


16 


17 


CSPO=0 

GP (MSOeEQ.1) GO TO 16 
CSPO=NAP*( E( 47)-1 0*( EC 50)* 08D +E(51)* OWL) 
CLDOC=NHP*E( 50 P*kOWD*( 141 £1 46 DF OTDPI )) 
CLIO=NHP*E( 51 )*OW1*(1t0 E( 46 )*OTIPI )) 
OCDM=0. 

OCIM=0. 

JC=0 

JC=JCti 

CDMT=0. 

OFMD( JC) =OP( JC )*C(8,IC) 
OFMHD-OFMHDt+CFMD( JC ) 

CDMT= GPL JC J*C(10,5C) 


9 
(408 ) se salem on aun sees 8, ° 
On ha ADLNONTAL ART meee 2 


080 4 6 


c 


18 


THE 


15% 


OCDM = OCDM +CDMT 

1F( JC.LT.JS) GO TO 17 

OCiN= E( 43)*OCDM 

OCR=(( CLDO+CLIG+CSPO+0C int+( CFOH( 1)/NPY)) 


i/OFMAD 9*+6(11,19t8(12,1) 


OC PUM=0 

JC=0 

JC=JICt+1 

OCPU( JC )=C( 10, J3C )#+C( 5S, IC *¥ CCR 
OCPUM= DMAX1( CCPUM,OCPU(JC )) 

IF( JCeLT»~JJ)}) GO TO 18 

IF (KPRESS.~NE.0) GO TO 721 

WRITE (6,603 ) 

WRITE (6,602) (OCPUC 1 ),1=1,)5) 
WRITE (6,605) 


TOTAL (MAK. }) PRODUCTION CAPACITY OF THE PLANT 


( INCLUDING GVERTIME) WITH EACH ADDITIONAL SHIFT, 
AMHDMX( Q,Y) 


F21 
24 


28 


{=0 

i=i+ti 

J=0 

J=Jti 

AMHDMX( J ,T J=4 B( 13, 1 0445-1 )9*8( 14,1) 9*( 61 44 9+1) 
iP (J21LTMX) GO TO 28 

IF (KPRESS.NEeO) GO TO 72 

WRITE (6,606) I 


© THE COEPPICIENT FCR THE SHIFT START-UP AND SHUT-DGY¥N 
C COST EQUATION, CSSA(Y) 


72 


IF( MXeEQe1) GO TO 119 
HMD2=E£( 44 )* B( 15,1) 

IF (BC 17,1).NE.0) GO TO 22 

BC 17, 1 )=(( OWI*OEL 9-H 2) 9*B( 13,1 )/ 


i( CHD*OCED*NHF) +E(2) 


22 CSSA( 1 )=0 


CSSACT )=( 42 0* BC 17, 1 * BC 11 )* EC 51 )*( HND2/ 
1B8(13,1)))*+0.25*(1-E( 39 ) )*E( 41 )*¥CPTH*( (HND2 


2/NHP J+ B(3)) 


IF (i-L1T-eNOY) GO TO 24 


119 IF (KPRESS.NE.0) GO TO 29 


WRITE (16,1112) 
1=0 


120 i=it1 


‘ “he Spent ay egg (Rly Me a; BB I-Moy 
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C RANGE CONSTRAINTS FOR THE INDEPENDENT VARIABLES 
PRODUCTICN, LCIRECT AND INDIRECT WORKFORCE LEVELS 


Cc 


WRITE (6,1119) CSSA(1I) 
IF (2.LT».NC’Y) GO To 120 


SHORTAGES COST COEFFICIENT, CINSA(J) 
29 JC=0 


3i 


BIRENG AND TRAINING COST COEFFICIENTS, CHTDA AND CHTIA 
73 CHTDAI=E( 21 )+( B( 22 )+( 1-EHDA1I 9*{ CPTH+E( 18 )* 

LEC S2)))*8( 16) 
CHT DA2Z=E( 21 )4+( E{ 22 )+( 1-EHDA2 )*( CPTHT+E( 19 )* 
LE(S2)))* 2017) 
CHTTAI=E{ 32 J)*+( E(33 )+( 1-EHLA1 )*( CPTHT+E( 30 )* 
LE( 53 )))*8(26€) 
CHT LA2=E( 32 J+( E( 33 )+( 1-EHIA2 *( CPTH+E( 31 )* 
LE( 53 )))*E{ 27) 


74 
37 


36 


JCH=JC+1 


CEiNSA( JC) = €(5,J5C ¥C(3,J3C 0401 6,5C )*E( 57) 
L*OCPU( JC o*( 14 E154 94E( 55 )) 400 7,5C )*C( 4,50) 


Fi JCsiT.JJ) GO To 31 

IF (KPRESS.~NE.O) GO TO 73 
WRITE( 6,609) 

WRITE (6,602 )(CINSA( 1), 1=1,J33) 


CHTLA>CHTIA 1* £124 )tCHTIA2Z*E( 25) 
iP (KPRESS.NE.2O) GO TO 74 

WRITE (6,617) CPTH 
CHIDA=CHTDAI*E( 12 D*+CHTDAZ*E( 13) 
WRITE (6,61) 

WRITE (6,1114) CHTDA,CHTIA 


JC=0 

JG=JCt+t1 

KMAX{ JC )=—0. 

KXMIN( JC )=S(1,5C) 

DO 26 I=1,KM 

XMAX( JC J=DMAXK1(S(1,J3C ), XMAX( JC) ) 
AMIN( JC J=DMINI(S(1,J3C),XKMIN{( JC)) 
CONTINUE 

XMAK( JCt+2) =E( 58 )*¥KMAK( JC) 

XMIN( JC#+2) =E( 48 )*XMIN( JC) 

LF( JC.LT~JJ) GO TG 37 

JD=0 

JD=JDt1 

UMseD( JD J=( 8113,3D 04¢ MX—-1 )*B(14,5D))/NHP 
iF( JDeLT.NOY) GO TO 36 
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(i Ate COLLEga) SEERW 
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32 
730 


77 


38 


150 
iSi 


i52 


153 
39 


C HERE 


40 
6000 


Y99 
Cc 


KL=NOY+1 
DO 32 1=1,Ki 

If (1-GT.NOY) GO TO 73vV 
WDMAX=DMAK1( UMHD{ 1), WDMAX) 
CONTINUE 

WDMAXK=1.3* WIMAX 
WIiMAX=WDMAX*( CWL/CHD ) 

GP (KRPRESS.NE.O0) GO TOG 77 
WRITE (6,614) 
XMAX( 1 )=WDMAX 
XMAX{ 2 J= WIMAX 
XMIN( 1 )=0.100 
KMIN( 2 )=0.1C0 
KJ=JI+2 
J=0 
J=Jt1 

pO 39 1=1,KJ 
KJL=( J—-1 )*K JTL 
KMAX(K JI J=KMAXK( I) 
XMIN( KJ J=XMAINC I) 

TFC i-2) 150, 150, 153 

GO TO (151, 152),1 

X( KJ 1 )=WDL J) 
GO TO 393 

X( RJT )=WIt J) 
GO TO 39 
XC KJ1T)=B8( J, 1-2) 
CONTINUE 

IF (J-LTsNH) GO TO 38 

THE VECTOR PROP( NM) IS COMPUTED 
X10=NH-1 

PROPXK=JPROPX 
PROPK=DLOG( PROPZ )/K10 
PROPK=DEXP( PROPK ) 
M1ii=NHt+1 

IP (1.GT.NH) GO TO 6000 
X10=1-1 
PRCOP( I )=PROPK** X10 

IF (KPRESS.»NE.0) GO TO 60 
DO 999 I=1,kKJI 

WRITE (6,615) XMAX(CG ),XMINCT),X( 1) 
CONTINUE 


C HERE THE PERIOD/MONTH LOOP HEGINS 


159 


om fae) 


O 4 G 
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60 NP=0 
42 L=Nit1 
DO 160 M=1,L 
IP (M.»GT.NH) GO TO 161 


ESTABLISH THE FORCAST DATA ON WHICH THE SDK GPTIMIZATION 
FOR THIS PARTICULAR PLANNING HORIZON YILL BE BASED 
J=0 
i156 J=Jt1 
TCM, J )=S(NPt+M, J) 
iF ( JeLT»JJ} GC TO 156 
160 CONTINUE 


HERE NEVAL iS RESET 
161 NEVAL=0 
&PFUINAL=0 
CALL ADAPTIVE PATTERN SEARCH ROUTINE TG 
OPTIMIZE CVYER THIS PLANNING HCRIZON« 
CALL PATS 
KPINAL=1 
CALL FCTIi 


STORE THE RESULTS} PLANNING DATA BASED UPON THE CPTIMAL 

DECISIONS FOR THIS PERICD/MONTH 
NP =NPt+i 
WDC NP J=xX(1) 
WICNP )=X(2) 
J=0 

162 J=Jt1 
RCNP ,J =X( 542) 
CVIT(NP,J)=CVITKIJ) 
VICNP,J)= Vik(J) 
VISP(NP,J)=VISPK(J ) 
IF (JsLTsJJ) GO TO 162 
MS( NP )=MSK 
FCUP( NP )=FCUPK 
CETI( NP )=CETIK 
CHTD(NP )=CHIDK 
CCP( NP )=CCPK 
CDM( NP )=CDMK 
CDW( NP )=CDWK 
CME( NP )=CMEK 
COH( NP )=COHK 
CASE( NP )=CASEK 
CETD( NP )=CETDK 
CHTI( NP )=CHTIK 
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CLPED( NP D=CLFDK 
CLPICNP J=CLEFIK 
COTD(NP )=COTDK 
COTICNP )=COTIK 
CINC( NP J=CINCK 
CINS(NP )= CiNSK 
CSPDCNP )=CSPRK 
CSPIC(NP )=CSFPFIk 
TCOP( NP )=TCCPK 
EDOCNP )=EDOK 
ELO(UNP J=EIOK 
OTDPCNP )=OTIPK 
OTIPCNP )=OTIPEK 
CHECK TO SEE if THE LAST PERIGD/MONTH IN THE 
TIMESPAN (NN) HAS JUST BEEN COMPLETED 
IP (NP.GE2NA) GO TO 50 
iP (N»LEKJ) GO TO 42 
NXIS=N-KJ 
pO 49 I=1,NxK18 
49 X(1)=X( 1TKJ ) 
GO TO 42 


PRINT GUT THE AGGREGATE PLANNING DECISIONS AND RELEVANT 
INFORMATION 
$0 WRITE (16,400 ) 
WRITE (6,418) 
WRITE (6,426 ) 
WRITE (6,422) MSO,0WD,OWI 
DO 43 K=1,NN 
I=(K-1)/NPY +1 
XWD=WD(K )*+0.5 
XWi=WI(K )+0.5 
43 WRITE( 6,421) 1,K,MS(K),XWD,XWI, FCUP( K) 
J=0 
44 j=Jjt1 
WRITE (6,40) 
WRITE (6,401) 
COVI{ J J=CVI(C J P*¥OCPUL J) 
WRITE (6,402) J 
WRITE (6,403) 
WRITE (6,418) OP{J),OVI(J ),COVI(J) 
DO 45 K=1,NN 
i=( K-1 )/NPYt1 
XP=R( K,J)+025 
KViI=VI(K,J)t+0.5 


39 


95 


56 
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WRITE (6,404) Ly Ky SC KyS 0g XP,XVL,CVIT( KyJ ),VISP(K,J) 
CONTINUE 

IF (J.LT.JJ) GO TO 44 

WRITE (6,400) 

WRITE (6,405) 

WRITE (6,406) 

DO 3S K=1,NN 

JD=( K-1)/NPYt1 

WRITE (6,407) JD,K,EDO( K ),EL0( K ),OTDP( K ),OTIP( x ) 
CONTINUE 

WRITE (6,400) 

WRITE (6,408) 

WRITE (6,410) 

J=0 

J=J+1 

Z1{ J )=0 

Z2J )=0 

Z3( J )=0 

Z4( 3 )=0 

Z5( J )=0 

Z6(J )=0 

Z7( J )=0 

Z8( JI )=0 

ZO 5 )=0 

Z10( J )=0 

DO 56 [=i1,NPY 

B=({ J-1 )*NPY+T 

WRITE( 6,411) J,K,CLPD(K ),COTD( K ), CHTD( KD), 


1CETD(K),CSPD( K),CLPA( K),COTI( K),CHTI( K), 
2CETI(CK),CSPI( kK) 


Z14J3) = Zi( J) +CLPD( Kk) 
Z26(5) = Z6(J) +CLPI(K) 
Z2(J) = 22453) +COTD( K) 
ZI(3I) = Z7I) *+COTE(K) 
Z34J5) = Z3(3) + CHTD( K) 
Z8(J5) = Z8( 3) + CHTI(K) 
Z4(5) = 2443) +CETD(K) 
Z9(0 J) = Z29(3) #CETI(K) 
Z54J) = 2543) +CSPD(K) 
410( 3) = Zi0Q({J) +CSPI(K) 


WRITE (6,412) ZLC5 42205), Z35( 5 042445 )9,25( 5), 


12613 ),Z7(5 )yZ8( 5 )4Z9(5),Z10( 5) 


iF (JeLTsNOCY) GO To 55 
WRITE (6,400) 
WRITE (6,413) 
WRITE (6,414) 
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57 


538 


300 
301 
303 
306 
307 
S00 
501 
503 
304 
400 
401 
402 
403 


J=0 
J=Jt+1i 
ZAI )=0 
Z21 J )=0 
Z3(J )=0 
ZAI =0 
Z5(J 9=0 
Z6( J )=0 
ZI J )=0 
Z8( J )=0 
Z3( 5 )=0 


DO 58 i=1,NFY 
K=( J-1 )*¥NPY41 
WRITE (6,415) J,K,yCDM(K ),CDW(K),COH{(K), 


ICIUNC(K ),CINS( K),CCP( K ), CASE(K ), CME( K),TCOP( K 3) 


Z2(143) = Z2(0 3) + CDW(K) 
Z3( 5) Z3 (J) + COH(K) 
ZidJjJ) = Z1 (J) + CDM(K) 
Z4(0 J )=Z415) + CINC(E) 
Z25( 3 9=Z5( 3) + CAINS kK) 
Z6(J3) = Z6(5) + CCP(K) 
Z27¢339=Z7(5) + CASECR) 
Za( J 9=Z8(5) + CMEC(K) 
Z9( J) =Z391 5 )t1TCCP(KR) 
CUMC=CUNCTZE( J) 

WRITE (6,416) 2105 9,220.5 ),23(5 ),Z2405),2515), 


1Z6( 5 ),Z715),Z819),Z91 5) 


iF( JeLT»NCY) GO TO 57 
WRITE( 6,417) CUMC 
WRITE( 6,400 ) 


FORMAT( 6610.0 ) 

FORMAT(G10.0) 

FORMAT( 213 ) 

FORMAT(1015,16) 

FORMAT( 16) 

FORMAT (6F1C.3) 

FORMAT( F10-3) 

FORMAT( 213) 

FORMAT 13 ) 

FORMAT (414 3} 

PORMAT(1X,* AGGREGATE PRODUCTION PLAN? ) 
FORMAT (//*PRGDUCT LINE -*,13) 
FORMAT (//6%,"PERIOD® ,6X%,'SALES',4X, 
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1? PRODUCTICN!,6X," ENDING INVENTORY',6X, 
2'*STOCK-CUIS'/13,' YEAR /', 

3'* MONTH! ,4X, *FORECAST',5X,'LEVEL!,8xX, 
439 LEVEL',6X,'S VYALUE*,6X,'{ PERCENT )'//) 

404 FORMAT (14,16,1X%,4(1%,F12.0),2X,F8.2) 

405 FORMAT (1XK,*WORKFORCE PERFORMANCE ANALYSIS? ) 

406 FORMAT (//21K%,'AVERAGE EMPLOYEE',9X,'E', 
L'MPLOYEE GVERTIME/IDLET INE / 
26X%,' PERIOD! ,6X%,* PERFORMANCE EFPFPICIENCY', 
314X,*{ PERCENT )'/1X,'YE*, 
4° apt,1X%,*/MCNTH! ,4X,*DIke WeFORCE IND. ', 
5°*WeFORCE’,4K,9DIR. WeFORCE INDe WeFORCE'//) 

407 FORMAT( 2X%,12,GXK,13,1%,21 8X, FGH02),2XK, 
12(8X,F6.2)) 

409 FORMAT( 1X,* PRODUCTION WORKFORCE COSTS ) 

410 FORMAT(//124,°{ - - - - —- DIRECT WoRK', 
1*FORCE COSTS - - - — — )8,8K,?{ - - - -! 
2,' INDIRECT WORKFORCE COSTS - - —- - )*/ 
34X%,'*PER.',3X,*REGULABR! ,15X,'HIRING/*,3%K, 
4° EMPLOYEE", 5X,!SHIFT* ,6OX,* REGULAR ,15X, 
S*HIRING/*,3%,*EMPLOYEE® ,5X,*SHIFT*/iX,'YR', 
61X,°/MOe',3X,y*PAYROLL', 3K, * OVERTIME" ,3X, 

FJ! TRAINING, 2X,' TERMINATION! ,2X%,*PREMIUM'’, 
85X,*PAYROLL',3K,* OVERTIME',3X,*TRAINING', 
92K," TERMINATION’ ,2X,' PREMIUM! //) 

411 FORMAT(13,14,1XK,2(1X,5F11.0)) 

412 FORMAT(//* *TOTAL*,1X,2(11X%,5Fi11.0)/) 

413 FORMAT (1iX,*12CTAL COST OF PLANT OPERATIONS") 

414 FORMAT( //15X,*DIRECT* ,6X,' DIRECT" ,16X, 

1* INVENTORY? ,5X,*BACK—!,3X,* PRODUCTION',16X, 
2*MARKETING*/5X,*PERIOD!",3X,'MATERIAL! 
3,5XK,* LABOUR’ ,5X,* OVERHEAD! ,4X,*CARRYING', 
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43X,*ORDERING! ,1X,' FLUCTUATIONS’ ,1X,* ADMINISTRA-',2x, 


S*aAND DISTRI-',4X,*TOTAL'/1X, 
67 YRe*,1X,!/MONTH' ,4X,*COSTS!,7X,'COSTS', 
17X," COSTS! , 7X,*COSTS' ,7X,"COSTS' ,6X,'COSTS' 


8,4K,'*TION COSTS*,3X,*BUTION COSTS',3X,*COSTS!//) 


415 FORMAT( 13,16,1X,9F12.0) 
416 FORMAT(/' * TOTAL’ ,3X,9F12.0/) 
417 FORMAT(//984X,°**TOTAL=' ,F16.0) 


4138 FORMAT(3X,41% ,5X,*O",10X,4---—*',i1X,3( 1k, Fi2.0),6X, 


1!---') 
419 FORMAT (1X,*AGGREGATE MANPOWER PLAN! ) 
420 FORMAT (//6%,*PERIOD NOs OF*,3K,*DIEECT ¥', 


14.FORCE TNDAIRECT W.FORCE CAPACITY UTIL.~'/ 


21X,* YEAR /MCNTH SHiFTS',7X,*LEVEL* ,14x, 
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3'*LEVEL*,10X,*FACTOR (%)*//) 

421 FORMAT (14,16,18,4X,2(F1000,5X%),2X%,F8s.2 ) 

422 FORMAT (3K, %1%,5X,*O",18,4X_2( F1020,8X), 
1625 %=-+*) 

601 FORMAT(//'ANNUAL FIXED CVERHEAD COSTS,CFOH{( YEAR ):*) 

602 PORMAT (6F12.2) 

603 FORMAT( //'MFG.e. COST/UNIIT (PROD. LINE) FOR THE PERIOD! 
1,*/MONTEH PRICK TO THE START OF THIS PLAN:24/*OCPU', 
2°91 FC) gSCHL, 2g GeovevevseseseadJis/ ) 

605 FORMAT( //*TCTAL (MAK. ) PLANT CAPACITY WITH EACH’, 

1* ADDITIONAL SHIFT *'*//5X%,11 SHIPFT*,4X,'2 SHIFIS',4x, 
2139 SHIPTSee vweessscesessece'/) 

606 FORMAT ('YEAR ="',13) 

609 FORMAT(//'SHORTAGES COST COEFF.:CINSA( JC),JC=1,2,020') 

610 FORMAT(//*HIRING. AND TRAINING COST COEFF. :', 
1//6X,'*CHIDAS ,5xX,'CHTIA'/) 

1112 FORMAT(//'*SHIFT START-UP AND SHUT-DOWN COST CUEFFe, ! 
1,'*CSSAC( YEAR Jz? ) 

1113 FORMAT (1X,f12-3) 

1114 FORMAT( 1%, 2F10.2) 

614 FORMAT(// GX, *XMAXK*',DX,*XMIN*',10X%,'XK!/) 

615 FORMAT (Gk, 73Fi2.2) 

620 FORMAT (//7X%,*EHDA1", 7X," EHDA2* , 7X," EHD( K)ocecseses!) 

621 FORMAT (//73%,'EHIAI*t, 7X, EHIA2Z*,7X,* EHIT( K Deoceccecne® ) 

STOP 
END 


PLANNING COST SUBROUTINE ,COMPUTES THE VALUE CF THE 
OBJECTIVE FUNCTICN AND GENERATES NECESSARY PLANNING 
INFORMATION BASED UPON OPTIMAL DECISIGNS 


SUBROUTINE ECTL 

IMPLICIT REAL*8 ( A-B,O-Z) 

DIMENSION K(12,10 ),V1(10),FMD(12,10),9D1 90), 
1Wi(90),EDO( 12),E10( 12), 0TDP(12),CT1P(12),A( 11), 
2PMHD(12),FMDK(10),FMHI(C 12), VI0O(10),CPUK(10), 
BSUTDP( 12 ),UTiP(12 ),CLPD( 12 ),CLPI(12),COoTD( 12), 
4COTi(i2),CHTD( 12),CHTi( 12), CETD( 12), CETI( 12), 
SMS( 12) 

COMMON 1T(12,10),8115,6),C(13,10),8(65),F(13), 
LAMHDMX(4,6),%(150), XKMAK( 150), XMIN(150),0VI(10), 
2V1IL(10,6),ViU(10,6),EHD(30),EHI( 30), 
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BCINSA( 10 ),CSSA(6 ),PROP( 12), VISPK(10), 
4V1iK( 10),SN,CFND( 10), CFOH( 6 ),OCPU(10),CVITK(10), 
5CCPK, CDK, CDWK, CMEK, COHK ,CASEK,CETDK ,CETIK,CPTH, 
6CHTDK ,CHTIK ,CLPDK,CLPIK,COTDK,COTIK, CINCK,CINSK, 
7CSPDK,CSPIK, ICOPK, EDOK, EIOK,OTDPK,OTIPK,OCPUM, 
SFCUPK ,OFMHD ,OWD,OW1,GED,CEI,CHTDA,CHTIA,SJ,SIGMA, 
QN, NP, NH, NHP,NCY,NPY,NWI,JJ,LDI,LIL,LIM,yMSO,MSK, 
+MXK ,NEVAL,KFINAL,KPRESS 

NEVAL=NEVAL+1 


LOAD PRODUCTION (R),DIRECT WORKFORCE (WD) AND INDIRECT 


WORKFORCE (Wi) DECISION VECTORS WiTH TRIAL VALUES FROM 


THE VALUES OF (K)'S SUPPLIED BY THE SEARCH ROUTINE (PATS) 


1=0 
143 I=1+1 
KJ=JI+2 
KJT=( 1-1 )*KJ 
WDC NPTTI )=xX¢( KJ 141 ) 
WIC NPti )=x( KJi+2) 
L=NP+1 
Jc=0 
142 Jc=scti 
RCL ,JC = X(KJI+2+J5C ) 
IF( JCeLT.«JJ) GO TO 142 
IF (IeLT.sNH) GO TO 143 
INITIALIZE THE FENALTY COST COEFRPICIENTS 
11A=0 
11C=0 
1ip=0 
11f=0 
11G=0 
11H=0 
12D=0 
I2F=0 
12cG=0 
12H=0 


SET UP THE (FIRST) PERIOD/MONTH LOOP, TO WORK OUT THE 
WORKFORCE AND SHIFT REQUIREMENTS FOR EACH PERIOD IN 
THE PLANNING HCRIZON 

KB=NHt1 

DG 165 M=1,K5 
IF(M.sGTeNH) GO TG S56 
KT=NP+M 

PMHD( MJ=0 

MSCM )=0 
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JD=( KT-1i )/NEYt1 
IF ( JD.GT.NCY) JD=NOY 
© COMPUTE THE MAN-HOUR EQUIVALENTS GCF THE PRODUCTION LEVELS 
JC=0 
145 JC=JCt1 
FMD{ M,JC J=R( Ny JC )*C(1 5,IC ) 
FMHO( M )=FMHD( M)tPNMD( M,JC ) 
VF( JC.LT.~JJ) GO TO 145 
© DETRERMINE THE NUMBER OF SHAFTS REQUIRED 
IP( MX.EQ.i) GO To 158 
166 iF{ PMHD(M).NEeCFMHD) GO 10 168 
MS(M)=MSO 
GO TO 160 
167 IF (M-1) 252,252,253 
252 IF ( FMHD(M).»NE.FMHDK) GO TO 168 
MS( MO=MSK 
GO TO 160 
253 iF( FMHD( M).EQeFMHD({ M-1)) GO TO 159 
168 DO 157 1=1,NX 
MS(M )=1i 
ITF( FPMHD( M).eLEsAMHDMX( 1,J3JD)) GO TO 160 
157 CONTINUE 
GO TO 160 
158 MSimM)=1 
GO TO 160 
159 MS(M)D=MS(M-1) 
C DETERMINE THE INDIRECT WORKFCRCE REQUIREMENTS 
160 IF(NWI.EQ.1) GO TG 162 
WiIFi=0. 
JC=0 
161 JC=-JCtl 
WiFi =WIFI +(C( 9, J5C D*¥R( M,JC )/NHP ) 
TF( JC.LT.JJ) GO TO 161 
WIF=( E( 7 )t+E( 10) )* WDC KT OTWIPFI 
1+( E( 8 )4+E( 8) )*¥ PMHD(M)/NHP+( MSC M)-1 0*E( 3) 
GP( WiFsLTe£i1)) WiIP=E(1 ) 
GO TO 55 
162 WAF=E( 4 )+E( 5 )*( WDC ET )** E16 )) 
55 FMHiI( M)=WiP*NHP 
GB (KFINAL»EQe1) GO TO 56 
165 CONTINUE 
Cc 
C CALL THE SUBROUTINE OVEFF TO DETERMINE THE AVERAGE ( WORK) 
© EFFICIENCY OF THE (DIRECT/INDIRECT) WORKFORCE EMPLOYEES 
Cc 


ciHVat KONTOVGORS BHT SO 


. | ° 

(A80¥) SOAeEVA ONE 2 
ESHYOIUMA 2 
er 2 


Ooo¢g OG 


eR om erie) 


168 


56 CALL CYEFF( NH, NP,OWD,9D,CEL,EHD,LDI,EDC,KFINAL ) 


CALL OVEFPF( NH ,NP,COWIL,WI,CEI1,EHI,LII,EIO,KFINAL) 
K=NP+1 


CALL THE SUBROUTINE WCOST TC DETERMINE THE WORKFORCE COSTS 
SUCH AS 3: BEGULAK PAYROLL, CVYERTIME, HIRING AND TRAINING 
AND EMPLOYEE TERMINATION COSTS 


CALL WCOST( NH,NP,NHP,E( 50), WD,OWD, EDO, 
iPMHD, E144), £145 ),E( 23), E136 ),E{ 46),CHTDA, 
2E( 37) ~OTDP ,CTDP,CLPD,COTD, 

SCHTD, CETD, NEY,NOY,KPINAL,I1DB,i1F,1iG,114H) 
CALL WCOST( NH,NP,NHP,E( 51), W1,OWI,EIO, 
LPMHI,E(44), 1 45),8( 34), 8138), 46),CHTIA, 
2E( 37 ),CTiIP,TIIP,CLPI,COTI, 

SCHTI ,CETI,NPY ,NOY,KFINAL, 12D, i12F;,12G, 120) 


SET UP THE SECOND PERIGD-LOOP, TC COMPUTE THE INVENTCRY 
COSTS,THE REMAINING CPERATING COSTS,THE VALUE OF THE 
OBJECTIVE FUNCTION AND THE PLANNING DATA 


KC=NHt1 

TCPP=0. 

DO i170 M=i1,KC 

IF (M.GT».NH) GO TC 183 
KT=NP+M 

JD=( ( KT-1 D/ NPY )t1 

iP (JDeGT.NCY) JD=NOY 


INVENTORY LEVELS FOR THIS PERICD, Vit{JC) 


i74 


152 


153 
144 
146 


147 
2350 


CincK=0 

CinNSK=0 

JC=0 

JC=JCt1 

UCP=0-. 

Cis=0. 

Cic=0. 

VISPK{( JC )=0 

iF (NP) 152,152,153 
UCr=OCCPU( JC ) 

GO TO 144 

UCP=CPUK( JC ) 

IP( KT-1) 146,146,147 

Vit JC I=OVIC JC DFR M,ICI-T(M, JC) 
GO TO 148 

4F (M-1) 250,250,251 

Vit JC J=VAK{ JC FRC My JC O- TAM, JC) 
Go TO 148 
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251 Vit JC J=VIOC( JC FRC My IC D- T( My JC) 
148 ¥VS=Vi( Jc )-VIL( JC,JD) 
IF( VS) 149,151,150 

i149 vS=-VS 
VISPK( JC) = VS/T(M,JC) 

CHECK FCR THE CCNSTRAINTS ON INVENTORY SHORTAGES 
Ip ( VISPK( JC )eLE.»C(2,;JC)) GO TO 139 
LIC=L1Ct+VISEK{ JC )—C1 2,30) 

INVENTORY SHCRTAGES COSTS 

139 CIS= VS*CINSA( JC) 

GO TO 151 
150 IF( VIC JC ).GIeVIUi JC,JD)) LiaA=Liati 
INVENTORY CARRYING COSTS. HERE TWO DIFFERENT METHODS HAVE 
BEEN USED TC COMPUTE THE INVENTORY CARRYING COSTS: (1) 
FOR USE IN GBJECTIVE FUNCTICN ,KFINAL=03; AND (2) FOR THE 
PURPOSE GF REPCRTING THE RESULIS, KFINAL=1. 

151i IF (KFINAL.NE.»i) GO TO 126 
CIC=Vij JC PX UC PXC(1,JC) 

CINCK=CINCK#C1iIC 
GO To 127 

HERE IS WHERE JPROPK IS USED 

126 CIC = VS*UCE*C(1,JC) 
PROPL=PROP( 4 ) 

LF( PROPL.»EQ.1i1) GO TO 154 
PROPL=PROPL*( CCPUM/OCPU(JC)) 

154 Cic=( PROPL-1 )*C(1,JC )*UCP*VYSt+CIC 

CUMULATE THE INVENTORY COSTS 
CINCK =CINCKt+CIC 

127 CIinsk = CINSK +CiS5S 
IF ( 3JCelilTeJ3) GO TO 174 

SHIFT START-UP AND SHUT-DOWN COSTS, CSSK 
CSSK=0 
IF (KT-1) 178,175,176 

175 IF (MS(M).ECe-MSOC) GO TO 171 
CSSK=CSSA( JD )*1ABS( MS(MI-MSOC ) 

GO TO 171 

176 iF (M—-1) 207,207,208 

207 IF (MS(M)eEQ.sMSK) GO TO 171 
CSSK=CSSA( JL )*1LABS(MS(M)-MSK) 

Go To 171 

208 iF( MS(M)eEQeNMS(M-1)) GO 10 171 
CSSK=CSSA( JD )* 1 ABS( MS( MI-MS( M-1 )) 

COST OF FLUCTULATIONS IN PRODUCTION LEVELS (CCPX) 

171 CCPK=0 
DMHD=0 
Jc=0 
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172 Jc=Jctl1 
LFCKT—-1) 177,177,175 
177 DMHD=DMED+DABS( FND( M, JC )- OFMD( JC)) 
GO TO 1798 
178 IF (N—-1)209,209,210 
209 DMHD=DMHD+DABS( FMD( My JC I-FMDK( JC )) 
GO TO 179 
210 DMHD=DMHD+DABS( FND( My, JC J-FMD( N-1,JC)) 
179 AF( JCeLT»JJ) GO TO 172 

PMHD=DMHD/F NHD(M) 
wpo=-0 
wio=0 

IF(KT-1) 18€,180,1381 

180 #WbDG=OWD 

wic-=owi 

GO TO 182 
181 wDO = WD(KT-1) 

wio=Wi( KT-1) 

182 DWD=DABS( WDC KT )-WDO) 

DWI=DABS( Wid KT )-Wi0) 

CSPIkK=0 
CSPDK=0 

IF(MS(M).EQ»1) GO TO 173 

SHIFT PREMIUM CCSIS, CSPDK AND CSPIK 

CSPIK=( E( 47 )-12 0 )*¥NHP*( EC 51 )*¥ WIC KT)) 

CSPDK=( E( 47 )-1-0 )¥NHP*¥( E(50 )*WDCKT)) 
STANDBY COSTS OF IDLE PACILITY 
173 CIPK =E( 43 )*UTDP( M)/NPY 

COPK=0.5*(0.5*( DWDtDWI )+PMHD 
1*( WDO+WIO ) )¥(1-E( 40 ) )*B( 42 )#CPTHTCSS& 
24E( S51 )¥E( 11 )*3.0* PMHD* WIC KT )+CIFS 

TOTAL PLANNING COSTS (VARIABLE PORTION OF THE PRODUCTION 
COSTS) 

CPPK=CLPD( M )}+CLPI(M )+COTD(M )+COTICM) 
1+CHTD( N)FCHTI( M)+CETD(M D#CETICN ) 
24+CCPKTC INCKFC INSK+CSPDKF+CSPIBK 

CUMULATIVE PLANNING COSTS FOR THE PLANNING HORIZON 

TCPP=TCPPTtCPPK 

THE FOLLOWING STATEMENT CAUSES A BYPASS (FOR KFINAL=0) OF 
THE FOLLOWING ( PLANNING DATA) ROUTINE DURING THE SEARCH 
PROCESS 
IF (KFINAL.NE«1) GO TO 216 


PREPARE THE PLANNING DATA (FOR BEPORTING) BASED UPON THE 
RESULTS OF THIS STAGE'S OPTIMIZATION 
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CAPACLTY UTILIZATICN FACTOR 
PCU=0 


PCU=FMEHD(M)/(1B(14,J3D +0 (MS(M )—-1 0%58(14,5D))) 


FPCUPK = FCU*100.0 
DIRECT MATERIAL COSTS 
CDMK =0 
JC=0 
200 Jc=JjJCt1 
CDMT=O0. 
CDMT = RCM, JC )*¥C(10,J5C) 
CDAK= CDMK+CDOMT 
TP (JCelTeJS) GO TG 200 
INDIRECT MATERIAL COSTS 
CIiMK = #(49 )*¥CDNK 
DIRECT WORKFORCE COSTS 
CDWK = CLPD(N) + COTD(M) + CSPDK 
INDARECT WCOREFOSCE COSTS 
CiWk = CLPI(N) + COTiI(M) + CSPIK 
OVERREAD COSTS 
COHK = (B(11,3D) + BC 12,35D))*FuMRDi NM) 
i+ Ci#wk + Cixk + CPOH{JD)/NPY 
GVERBEAD RATE ( PER PRODUCTION MAN-HOUR ) 
COHR = COHK/PMHD{ MN) 
BiRECT YVORKFPCRCE COST PER PRODUCTION MAN—-HOUR 
CDWH = CDWK/FMHD( M) 
COST OF MANUFACTURE PER UNIT 
JC=0 
204 JC=JCt1i 
CPUK( JC) = €(10,35C) + C(8,J5C )*( CDWH+COHR ) 
VALUE OF THE ENDING INVENTORIES 
CVATK( JC O=CPURG JC D*¥ VI TC) 
CVITK( JC D=DABS(CVITK( JC )) 
iF (JCsLT.JJ) GO To 204 
TOTAL COST GF PLCDUCTICN 
CPTK=O06 
CPTK = CDWKTCOHKTCDMKTCINCKt+C INSK+CCPK 
ADMINISTRATIVE, SELLING AND MARKETING COSTS 
CMEK=CPTK*E( 55) 
CASEK = CPTK*E(54 ) 
TOTAL COST OF OPERATIONS 
TCOPK= CPTK+CASEKTtCNEK 
PREPARE REMAINING DATA FOR STRORING 
iF (OTDP(1).EQ.0) OTDP( 1 )=-UTDP( 1) 
TF ( OTIP(1)-EQ.0) OTIP( 1 J=-UTIP(1) 
OTDPK= GTDP(1)*100. 
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GOTIPK= CTIP(19*100. 
CLPDK= CLPD(1 ) 
CLPIK= CLPI(1) 
COTDK= COTD(i1) 
COTIK= COTiI(1) 
CHITDK= CHTD(1) 
CHTIK=CHTI({ 1) 
CETDK= CETD(1) 
CETiIK= CETI(1) 
EDOK= EDO( 1) 
ELOK= EI0(1) 
FPMHDK=FMHD{ 1) 
MSK= MS(1) 
J=0 
2ii J=Jti 
FMDK(J)= FMD(C1,J) 
VIK( J )=VI(s ) 
IP (J.sLT.JJ) GO TO 211 
iF (KFINAL. EQC.1)KRETUEN 
2io iF (KPRESS.~NE.0O) GO TO 213 
WRITE (6,402) JD,KT,¥Dt «T),MS(M), Wii KT) 
WRITE 16,403) 
J=U 
217 J=Jt1 
iF (NP) 218,218,218 
2i8S WRITE (6,404) JTC My J), RIM, I),VI1 I), VISPKC J ),OCCPU(J) 
GO TO 220 
220 IF (JeLTeJJ) GO TO 217 
WRITE 16,901) 
WRITE (6,902) KT,CINCK, CINSK, CSSK,CSPDK,CSPIK, 
1CCPK,CPPK 
213 J=0 
2i4 Jj=Jt1 
Vic(J )=Vil J) 
IF (JeLT.JJ)}) GO TO 214 
170 CONTINUE 
183 IF( NEVALeNE«oil) GO TO 117 
C COMPUTE THE RATES FORK CHARGING THE PENALTY COSTS 
DO 116 1=1,10 
ACI J=TCPP*F(TI ) 
116 CONTINUE 
© VALUE OF THE CBJECTIVE FUNCTION 
117 TCPK=TCPPtHtI 1A¥*¥A( 1 FT 1C* AC 2 FTL D* A 3 FTL P*AC 4 ) 
LI+G2D*A(5 JFL 2F*AL 6 PFALIG* ATT JFL 2G%*AC 9 DF ITLH*ALS ) 
2+ i12H*A(10) 
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SJ=TCPP 
SN=TCPK 

402 FORMAT (////1K,"YEAR —*,12,3XK,' MONTH/ PERIC® 
i,7D 34,14,6X,*DIRECT WORKFORCE LEVEL', 

a* 2*,Fi0.0/f/12X%,"NOs»s GF SHIFTS 24,14,6X, 
S'tINDIRECT WCRKFORCE LEVEL +4,F10.0) 

403 FORMAT (//234%, *PRODUCT—',4X,'SALES',6X, 
i*PRODUCTION',1X,* END. INVENT.*, 2X, 
2'STOCK-CUIS* ,4%,* UNIT—-MFG.*t/* ',3X, 

StLineE® ,4X,* FORECAST’, 7K,*LEVEL* , 7X, 
4'* LEVEL!, 7X, *'( PERCENT )*,6X%,'COST'/) 

404 FORMAT (4X%,13,2XK,3( 1X,F12-0),2K,21 F802,6K)) 

901 FORMAT (1 //7%,* INVENTORY! ,3%,* BACK, 3K,*SHIPT*, 
1* START—4,2%,'COST OF SHIFT*,1X,*COST OF', 
2*' CHANG—* ,4X,*TOTAL'/1XK, *PER.*,3K,*CARRYING', 
32X,* CRDERING*,1K,° UP/ SHUTI-DOON® ,4X,'PREMIUMS', 
44X,'ING PROLUCTION',3X,*COST OF'/1X,'MO.', 
54%,'* COSTS! , 5k, 'COSTS* ,6%,*COSTS! ,4K,'D/WF.! 
by, € 1/Wek.? 94K» *LEVELS? 7OX,? PRODUCTION /) 

902 FORMAT (14, 2X%,611X%,F9.0 ),3X,F12.0) 

RETURN 
END 


LEAEFP SUBKCUTIAE, DETERMINES THE EMPLOYER EFFICIENCY 
DURING LEARNING/TRAINING PERIODS 


SUBROUTINE LEAEFPP( RA, HNT1,HNT2, EH1, EH2, 
L1PHLI, PH2,NHP,ET,1T11,HHA1,EHA2 ) 

IMPLICIT REAL*8( A-H,0-Z ) 

DIMENSTON A(2 ),AA(2),8( 2,30 ),BT(1),b8(2), 
AHNT( 2),TNK( 2),TNP( 2), EH( 2), PH( 2), EHA( 2) 
HNT( 1)=HNTI 

HNT(2)= HNT2 

EH(1)= EH1 

EH(2)= EBH2 

PH(1)= PHI 

PH(2)= PH2 

PNHP=DPLOATCNHP ) 

BB=( DLOG( RA ) )/DLOG(2.B0) 

H=-BB/(BHT1 ) 

po 1 J=1,2 

AA(J J=(CBBt1 DX HNT( JS *C ENC J **(1/5)) 

A(J )=( BBt+1)/ AAJ) 

BEHA( J )=( BAC J )*C ACS OK HNTCS ) **B)/( Bt1) 
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TNP( J)=ANT( J) /PNHP 

KK( J )=LPIX( SNGL(TNP(J))) 

IFC KK( J )eLTe1) GO TO 2 

I=0 

3 E=Eti 

E(d.J3 i J=( EH( 3 )*( ACS )** 8B) DX CNP ST De %*( Bt 1 - 
1( NHP*{ I-1 ) )**(Bt1 ))/(NHP*({ Bt+1)) 

GF ( G.1LT.KK({J)) Go TO 3 

LTF(TNP( J ).EQeKK(J))GO TO 1 

2 TNK( J J=HNT( J -C KK( J )*¥NHP ) 

E( J ,KK( 9 )+1 J=(( BR(.5 *C AC 3 **R )*( CHNTC I) *% 
1( B+i ))-( (NHP¥KK( JS ))**( 541 9900/0BT1 D+ NH P- 
2TNK(.5)))/NHP 

K&{ J j9=KAK(J)t1 

1 CONTINUE 

EHA1I=EHWA(1) 

EHA2Z=EHA( 2) 

TP( KK(1).GE.KKi 2))GO To 4 

Lit=Kxk( 2) 

GO TO 5 

4 T1IiI=Kk(1) 
§ iIF( if1eGT.1) GO TO 7 
ET( 1)—=PH(1)* 8611, 19+PH( 2 D*E{ 2,1) 
GC TC 8 
7 po 6 1=1,111 
TFC 1.«GT.»KK{ 2) DEI 2 I)=1 
BTC i J=PH(1)*E( 1,1 PHI 2 * EL 2,1) 
6 CONTINUE 
8 BETURN 
END 


SUBROUTINE OVEFE, DETERMINES THE AVERAGE EFFICIENCY OF 
THE WCREFCRCE 


SUBROUTINE CVEFF( NH, NP, OW, W¥,OEK,E,111,BO,KFINAL) 
IMPLICIT REAL*¥8( A-H,0-Z ) 

DIMENSION E(1 ), E011), ¥W(1),89(50),5790) 

BI=OE 

J=LLItNH 

cWi =Wi1)-cw 

KG=NHT1 

pO 1 M=1,KO 

iF (M.GT»NH) GO TO 18 
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OS=C¥ 

HW( i )=CWi 

$8(i)=wWwi) 

KT=NPtM 

1=0 

i=itl 

K=NPt*NHt1-1 

IF (KeLE»1i) GO TO 26 
pw( K )=W8( K -W( K-15 

S(K )=8(K) 

IF (1.LE.J) GO TO 25 
EO( M)=0 

NT=0 

N=2iit1 

DO 28 1=1,N 

IF (1.GT.111) GO TO 29 
L=1 

NT=KTt1-i 

IF (NT.LE.O) GO TO 2D 
1F( HW(NT).GE.0) GO TO 10 
LP(( NT-L)eGi1e0) GO TO J 
OS=HW( NT +08 

HWC NT )=0 

GO TO 106 

IF (LeLEe( ifi-1)) GO To 4 
SCNT-L )=8¥W{ NT OtS( NT-L) 
HW( NT )=0 

GO TO 10 

HW( NT J=HW( NT I+HW(NT-L) 
LTF( HW( NT ).GE.0) GO TO 9 
HW( NTL )=0 

L=L+1 

GO TO 2 

HW{ NT—L J=HW( NT) 

HW( NT )=0 

BO( M)=EO(M )FHWCNT * ECT) 
CONTINUE ' 

iF (NT-1) 30,930,931 

1Ff ( KT.«GT.2) Ei=1.0 
EO(M )=( BO(M )+OS* EI D/W( KT) 
Go To 1 


EO( M)=( EOC M J+SCNT-1))/WCKT) 


GP ( KFENAL»~EQe1) RETURN 
CONTINUE 
RETURN 
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END 


WCOST SUBRCUTINE, DETERMINES THE WORKFORCE COSTS, 
HIRAING/TERMINATION AND OVERTINE REQUIREEMENTS 


SUBROUTINE WCOST( NH,NP,NHEP, ALB, W, OW, EO 
1,FMHW, POTU, PULIU, PHTU, PTU,ROT,CHIA,PTEP, 


20TP , UTP, CLP ,COT,CHT,CET ,NPY,NOY, 
SKPINAL,I1D, 11F,11G,118) 


IMPLICIT REAL*¥&( A-H,0-Z ) 
 DIMENSICN EC(1),O0TP(1), UTP(1),WTP( 30), 
1PMHWC1),CET(1),CHT( 1),COT(1),CLP(1),9(1)9 

KOQ=NH+1 

DO 60 M=1,KC 

If( M»GT»sNH) GO TO 65 

KT=NPt+M 

LS=8+M 

JD=(( KI-i/NPY)+1 

IF (JD.GT.NCY) JD=NOY 

WOo=0 
WT=0. 

WH=0 
OTH=0 
OTP(M)=0- 

UTP( M)=0. 

UTH=0. 

WTP( LS )=0. 

COT(M )=0 
CHT{ M )=0 

WHP=0 
CET( M)=0 


C REGULAR PAYROLL COSTS: 


CLP( M)=ALK* 4(C KT * NHP 


C LEVEL OF OVERTIME/IDLETIME AND CYERTIME COSTS: 


52 


OTH=( FMHW{ MW )/ EOC M ))-WCUKT D*NHP 
LF( OTHI52,63,53 

UTH=-OTH 

OTH=0 

UTPCM )=UTH/( WCKT D*NHP ) 
DUP=UTP{ Mw I-PUTU 


C CHECK FOR THE WCRKFORCE iDLETINE CONSTRAINTS 


IFC UTP(M).LEsPUTU) GO TO 63 
LiF=11F+*+DoUP 
GO TO 63 
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53 OTP( M)=OTH/(CWCUKT D*NUP ) 
DOP=OTP(M )-(0.75* POTU ) 

C CHECK FOR THE WCRKFORCE OVERTINE CONSTRAINTS 
IFC OTP(M).LE.»POTU) Go TO 74 
Gip=-iibdt+pDcP 

C OVERTINE CCSIS 

74 COT( MN )=ROT* ALR*OTH 
C EMPLOYEE HIRING AND TRAINING AND TERMINATION COS’IS: 
63 IFC KT-1) 61,61,62 
61 FO=oW 
GO TO 66 
62 #0=W(KT-1) 
66 WH=W( KT )-WO 
LPC WH )55,64,58 
C LEVEL GP EMPLCYEE TERMINATIONS 
55 #T=-Wil 
WH=0. 
WIP( LS )=WT/ WO 
DTP=WTPCLS )-PTU 

© CHECK FOR THE LEVEL OF ENPLOYEE TERMINATIONS 
iF (WIP(LS).LEePTU) GO TO 59 
Lik=i1LHt+DOTP 

59 KY=0 

56 KY=KYt1 
PIK=PTRtWIP( LS+1-i ) 

57 CONTINUE 
PTR=PTR/KY 
PLM=PTE-PTRE 

C EMPLOYEE TERMINATION COSTS: 
CET( MI=CHTA*( BWItPLM *89( KT)) 
GO TO 64 

C EMPLOYEE HIRING AND TRAINING COSTS: 

58 WHP= WH/¥C 
DHP=WHP-PHTU 

© CHECK FOR THE CCNSTRAINTS ON THE RATE OF HIRINGS 
If (WHP.LE.»PHTU) GO TG 67 
11G=11G+DHP 

67 CHT( M)=CHTA*DH 
64 IF (KFINAL-EQe1) RETURN 
60 CONTINUE 
65 RETURN 
END 
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SUBROUTINE PATS 
PATTERN SEABCH SUBROUTINE. THIS SUBROUTINE WAS TAKEN FROM 
TeWeSiKES'S UNPUBLISHED PH.».D. DISSERTATION REFERENCE (40) 
WHO OBTAINED LI FROM ANAND KHOKHA OF UCLA IN FEB 1970. 

"IT iS BASED ON THE PATTERN 

SEARCH ROUTINES OF HOOKE-JEEVES, WEISMAN-WOOD-RIVLIN, AND 
TADUBERT. ALTHOUGH THE SUBROUTINE CONTAINS CERTAIN 
INEFFICIENCIES AND REDUNDANCIES (AS DOES THE MAIN 
ROUTINE), THERE ARE NO BUGS WHICH CAUSE ERRCKS. BECAUSE 
OF THE GREATER CCMPLEXITY OF THE LOGIC, A MORE COMPLETE 
ANNOTATION AND GLOSSARY IS PROVIDED FOR THIS SUBROUTINE. 
THE CODE IS BIViIDED INTO 15 CONTINUOUS SUBSETS CALLED 
BOAES, WiTH ENTRANCES AND ERITS FOR EACH BOX EXHAUSTIVELY 
IDENTIFIfD. THIS FORMS A COMPLETE BASIS FOR A FLOW 
CHART, IF THE NEED SHOULD ARISE. THE BEST DESCRIPTION OF 
THE PATTERN SEARCH APPEARS IN WILDE, ‘OPTIMUM SEEKING 
METHODS’. SEE ALSO REFERENCES 38 AND 35." 


THE FOLLOWING STATEMENT CAUSES DOUBLE PRECISION 
ARITHMETIC. 


IMPLICIT REAL*8 (A-H, O-Z) 

DIMENSION D(20),0(720) 

COMMON T412,10),8(118,6),C{13,10),8 65), F113), 
LAMBDMX( 4,6 ),X(150 ), KRAX( 150 ), XKXMIN(1150),0V1T(10), 
2V1L(10,6),V10(10,6),EHD( 30),EHI{ 30), 

SCINSA( 10 ),CESA(6 ),PROP( 12), VISPK(10), 

4V1K( 10),SN,CPMD( 10), CFOH(6 ),O0CPU(10),CVITK( 10), 
5SCCPK,CDMK, CDWK, CMEK, COHK, CASEK, CETDK,CETIK,CPTH, 
6CHTDK, CHTIK,CLPDK,CLPIK,COTDK,COTIK,CINCK,CINSK, 
7TCSPDK ,CSPIK y TCOPK, EDO K, EIOK, OTDPK, OTIPK, OCPUM, 
SPCUPK, OF MAD ,CWD, OWL ,OED,CE1,CHTDA,CHTIA,SJ,SIGMA, 
ON,NP,NH,NHP ,NOY,NPY,NW1I,JJ,LDi,LI1,LIM,MSO,MSK, 
+MX,NEVAL,KFINAL,KPRESS 


BOX 1,CARDS THRE * PURPOGSE- TO INITIALIZE DATA» 


ENTRANCES—- FROM START OF SUBROUTINE AND (DUMMY) FROM BOX 
10 ViA STM 100. EXIT- TO 190 IN BOX 52 COMMENT- VERY 
STRAIGHTFORWARD. 


100 TOL=1.E-5 
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DEL=.1 

BET=.15 

LT=0 

LSN=0 

NPF=0 

iTR=1 

ALP=2.s 

GR=2.2 

LA=1 

K=1 

RY=K 

LT2=0 

LT3=0 

LT4=0 

LT5=0 

LT6=0 

LT7=0 

4b1=0 

Lb2=0 

1D3=0 

ip4=0 

KOUNT=0 

V=0. 

OLDV=O0. 

NEVOLD=0 

DO 180 i-i,s 

OC 1)=x(1) 
NOTE FROM THE NEXT STATEMENT THAT THE INITIAL VALUE OF 
Di K) DEPENDS ON HOW FAR APART THE UPPER AND LOWEK 
BOUNDARIES ABE. 


180 pi 1 )=DEL*( XMAX( I )-XMIN(TI)) 
CALL FCT1 
SNOLD=SN 
GOLDSN=S5N 
SC=SN 
SP=SN 
Mi=1 
a 2=1 
K=1 
KV=K 
KK=1 
L4=1 
iF (KPRESS.NE.0) GO TO 130 
WRITE( 6,999) DEL, BET 
999 FORMAT(40', 20K," PATTERN SEARCH DEL=",F622,' BET', 
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1?=",F6.2) 

WRITE( 6,991 JNEVAL,KOUNT, LT2,LT3,LT4,LT5,L16,L17,K,KK, 
iM1,M2,NPF,LA,1L4,1D1,10D2,1D3,1D4,ySN,SJ_, V, GR 

WRITE( 6,992) KV, (D(i),1=1,N) 

GO TO 190 


BOX 18, CARDS THRU « PURPOSE- CALL OBJECTIVE 
FUNCTION. ENTRANCES-— FROM BOXES 17, 7, AND 14. EXITS- IG 
STMS 100 AND 500 (DUMMY), TC 282 AND 285 IN BOX 2, TO 463 
AND 466 iN BOX 11, AND TO 580 IN BOK 15, TO 510 IN BOX 
13. COMMENT- LA CONTROLS TRAFFIC THROUGH THIS BOX, SEE 
GLOSSARY FOR INTERPRETATION OF VALUES OF LAs 


270 CALL FCTI1 
GO TO (100, 282,463,580 ,285 ,466,510,500),LA 


BOX 2,CARDS THRU ° PURPOCSE- TO EVALUATE A FORWARD 
MOVE. ENTRANCE—- FROM BOX 18 VIA STMS 282, 285-2 EXITS- TO 
360 IN BOX 6 WHEN FWD MOVE FAILS, TO 300 iN BOX 3 WHEN 
FWD MOVE SUCCEEDS. COMMENT- WE REACH THIS BOX WHEN LA=2 
OR Se -- FIRST INCREMENT STATE COUNTER-— 


282 LT2=LT2+1 
GO TO 280 
285 LIS=LT5ti 
—-- THEN TEST FCR SUCCESS (250) OR FAILURE (360) 


280 IFC SN-SP )29C, 360, 360 
290 GO TO (300,2892),L4 

—-- NEXT IF WE ARE IN AN EXPLORATORY MODE FOLLOWING A 
PATTERN MOVE ( ATTEMPT), WE MULTIPLY THIS SUCCESSFUL 
FORWARD STEP SiZE SY ALP. 


292 D(K J=D(K D*ALP 


BOX 3, CARDS THRU » PURPOSE - TO RESET SCONE 
VARIABLES ENTRANCES- FROM BOXES 2 AND 11, BOTH VIA SiM 
300. BXIT- TO 305 IN BOX 4 COMMENT- STRAIGHTFORWARD -- 
SEE GLOSSARY FCR DEFINITIONS. 


300 SP=SN 
M2=1 
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wMi=1i 
NPPF=0 


BOX 4, CARDS THRU 2 PURPOSE - TO DECIDE WHETHER 
TO TEST FOR A BREW BASE POINT. ENTRANCE~- FROM BOXES J3 AND 
12 ViA STN 3905. EXITS- TO 340 IN BOX 7 WHEN A BASE POINT 
TEST iS REQUIRED, -TO 200 IN BGX 5 GTHERWISE. COMMENT—- 
THIS BCX iS THE CNLY ENTRANCE TO THE BASE PCINT TEST, BOX 
Je ALSO NOTE THAT WHEN LT=0, AS IT NORMALLY IS (SEE 
GLOSSARY), IT IS ONLY HERE THAT KK IS INCREMENTED AND 
TESTED. 


JOS A=Kt1 

777 K=1 

320 IF( LT)9330,350,340 
J30 KK=KK+1 


BOX 3S, CARDS THRU a PURPOSE - SET SOME INDICES 

AND INCREMENT THE VALUE OF i(k) IN PREPARATION FOR 
TESTING A PORWARD MOVE IN THE & TH VARIABLE. ENTRKANCES— 
FROM BOX 1 ViA STM 190, BOXES 4 AND 8 ViA Sim 200, BOX 10 
Via STM 180, BOX 15 VIA STH 210-6 EXITS- TO 4590 IN BOX 12 
WHEN VARIABLE 18 FIXED AND CANNOT BE PETURBED, ITO 2350 IN 
BOX 17 (BOUNDARY CHECK) OTHERWISE. COMMENT- HERE WE 
ADJUST X(K) TO TRY A FORWARD NOVE. 


200 GO TO (190,210),L4 
ivv L4=1 
LA=2 
GO TO 220 
210 La=5 
HERE'S WHERE WE CHECK FOR VARIABLES WHICH ABE FIXED. 


220 IF( D(K})9225,490,225 
225 X(K )J=XK(K IFD(K) 
GO TO 230 


BOX 6, CARDS THRU e« PURPCSE- WHEN FWD MOVE FAILS WE 
REACH BOX 6 ANL SET THE VALUE OCF Xik) IN PREPARATION FOR 
ATTEMPTING A REVERSE MOVE. ENTRANCES- FROM BOXED 2 AND 
17, BOTH VIA SIN 360. EXIT- TO 230 IN BOX 17 ( BOUNDARY 
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CHECK) ALWAYS~s STMS 360 THROUGH 350 SHOULD BE "SET, THEN 
TESTED AND RESET if NECESSARY' WHEN REPROGEAMMINAG FOR 
SPEED. 


360 GO TO (370,380),L4 
370 LA=3 
GO TO 3390 
380 LA=6 
390 Xi K)J=X( KB I-2* DUK) 
GO TO 230 


BOX 17, CARDS THRU « PURPGOSE-BOUNDARY CHECK IN 
PREPARATION FOR POSITIVE OK NEGATIVE MOVE GF X11 K)e 
ENTRANCE- FRCM BOXES 5 AND 6 Vid STM 230.2 EXITS- To 270 
GN BOX 18 (FCT1) WHEN BOUNDARY CHECK I5 PASSED, TIO 500 
(DUMMY) IN CERTAIN EVENTS WHEN FAILED, TO 360 IN BOX 6 
WHEN POSITIVE MOVE FAILS THE BCUNDAKY CHECK, TO 480 IN 
BOX 2 WHEN NEGATIVE MOVE FAILS THE BOUNDARY CHECK. 
COMMENT—~ 270 INPLIES BOUNDARY WAS NOT VIOLATED. 


230 IF( X( K O-XMAX(K) 9250, 270, 260 
250 IFC XMIN( K )-X0 KB) 270,270,260 
260 GO TO (500,260,450,500,360,450,500,500),LA 


BOX 7,CAKRDS TERU ® PURPCOSE-TIEST TO SEE WHETHER A NEW 
BASE POINT HAS BEEN ESTABLISHED. ENTRANCE- FROM S5OX 4 VIA 
STM 340. EXITS —- TO 400 IN BOX 8&8 WHEN TEST FOR NEW BASE 
POINT PAILS, TC 270 IN BOX 18 (FCTI) WHEN NEW BASE PCINT 
IS ESTABLISHED. COMMENT~- FikST, NOTE THE TOLERANCE OF 
20001 IN STH 340-- iS THIS THE BEST VALUE-e NEXT, SEE 
NOTE UNDER M2 AND NPF IS GLOSSARY. ONE TROUBLE WITH THE 
EXISTING CODE iS THAT IN SIM 340 £400,400" COULD BE 
£353,353! AND/GR 'S* IN STM 93953 IS TOO HIGH. LOOK [INTO 
THIS SGMETIME. FCR EXAMPLE, IF 1T IS DESIRABLE TO LET NPF 
GROW ABOVE 1, SUCCESSIVE PATTERN MGVE ATTEMPTS (BEFORE 
THIS PATTERN I5 DECLAKED BROKEN) MIGHT TRY INCREASING Tik 
D(K) AND/OR DECREASING THE VALUE OF GR. IN ANY CASE NPF 
PROBABLY SHOULD ALWAYS BE LESS THAN 3 OF 42 AS THE CODE 
NOW STANDS STMS 350 AND 353 ARE BREDUNDANT -- THE BRANCH 
FROM STM 340 TC 350 COULD LEAD DIRECTLY TO STM J352 
INSTEAD, SiNCE UNDER THE CODE'S EXISTING LOGIC NPF CANNOT 
BECOME GREATER THAN 1. 


340 1f( SP+.0001-SC )350,400, 400 
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350 MEANS WE HAVE SUCCESSFULLY IDENTIFIED A NEW BASE 
POINT. iF L4=1 WE SHOULD NCW ATTEMPT THE FPAikRST PATTERN 
MOVE IN THIS SERIES. 


350 GO TO (352,253),L4 
352 LA=7 

Mi=1 

GO TO 270 
353 IF( NPE-5 352,400,400 


BOX &,CARDS TERE 2» PURPGSE- HAVING FAILED TOG 
ESTABLISH A NEW BASE POINT, WE REACH THIS BOX TO DECIDE 
WHERE TO GO NEXT. ENTRANCE~- FROM BOX 7 VIA STM 400. 
BXAITS- TO 200 IN BOX 5 WHEN KKLE.~N APTER BASE POINT TEST 
(-IMPOSSIBLE WiTH LT=0),T0O 410 IN BOX 9 WHEN FAILURE 
OCCURS AT THE CONCLUSTON OF ADJUSTMENTS POLLOWING AN 
ATTEMPTED PATTERN MOVE, 20 440 IN BOX 10 WHEN IN THE 
PROCESS OF SEARCHING POR A NEW PATTERN AND THE VALUE OF 
Mi (SEE GLOSSARY) I5 LESS THAN Nt1, TO 500 IN BOX 16 WHEN 
NO NEW PATTERN IS FOUND (SEE Ni IN GLOSSARY) AND SEARCH 
IS TO BE TERMINATED.» COMMENT- CONCERNING STMS 770 AND 
404, SEE FINAL SENTENCE UNDER KK IN GLOSSARY. STM 420 18S 
REDUNDANT, BUT THE STM NUMBER IS AGT. 


400 IFC(LT)404,404,778 

778 BK=KKt+1 

404 LF( &K-N)200,200,773 

773 GO TO (429,410),14 

423 L4=2 
I1F(M1-—-N 94406 ,440,500 


BOX O,CARDS THEU « PURPCSE- SESTOGRE THE VALUES OF 

X€ 1), APTER UNSUCCESSFUL ATTEMPT TO ESTABLISH NEW BASE 
POINT. ENTRANCE- FROM BOX 8 VIA STM 410-4 EXIT- TO 440 
IN BOX 10. 


410 SP=SC 
420 X(1)=Q11) 
NPF=0 


BOX 10, CARDS THRU e PURPOSE- TO RESET SCNE 
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COUNTESS-- SEE GLOSSARY FOR DEFPANITICNSe. ENTRANCE- FROM 
BOXES 8S AND &, BCTH VIA STM 440. EXIT- TO 190 IN BOX 5- 


440 M2=1 
Mi=1i 
KK=1 
GO TO 130 


BOX 11, CARDS THRU « PURPOSE- TO EVALUATE A KEVERSE 
MOVE- (CORRESPONDS TO BOX 2'S FWD MOVE EVALUATICN). 
ENTRANCE- FRCM BCX 18 (FCT1) ViA STMS 463,446. EXITS TO 
480 IN BOX 12 WHEN REVERSE MOVE FAILS, TO 300 IN BOX 3 
WHEN REVERSE MCVE SUCCEEDS» CGMMENT—- NOTE THAT AS USUAL 
WE SET THE COUNTERS LT3 AND LI6, DEPENDING ON WHETHER WE 
REACHED THIS BCK WITH LA=3 GR 6. 


463 LI3=LTSti 
GC TC 460 
466 LTG=LT671 
460 IF( SN-SP )470,480,480 
SUCCESS—-—- THEREFORE iN THE FUTURE A FORWARD MOVE WILL BE 
WHAT HAS HERETCFORE BEEN A BACKWARD MOVYE-- HENCE, Sim 470 


470 D(KJ=-Di XK) 

GO TO 300 
BOX 12, CARDS THRU » PURPOSE- AFTER REVERSE MOVE 
PAILURE THIS BCX REDUCES STEP SiZE AND ATTEMPTS IC TRY AN 
EXPLORATORY STEP FOR ANOTHER VARIABLEs ENTRANCES—- FROM 
BOX 5S VIA STM 480, FROM BOXES 17 AND i1, BOTH Via STMS 
480.  EXII- TO 305 IN BOX 4. 


480 X¢(K )=XC(K FDI ) 

D(KI=DC KB )*BET 
OBSERVE THIS USE OF THE VARGABLE TOL. WHEN REPROUGKANMIN- 
G, THE DOUBLE PRECISION ABSCLUTE VALUE FUNCTICN SHOULD BE 
REPLACED WITH A COUPLE OF EASY FORTRAN STATEMENTS» 


DX= DABS( X( Kh )/D( K D* TOL) 

IF( 1-pxX )481,48&2,484 
481 D(C Kj)=D( K)*DK 
HOW WAS "128-304 SELECTED? IS IT OPTIMAL WITHIN A FACTOR 
OF 1000? 


482 DX= DABS(1.E-30/D(K)) 
iF( i-DxX )483 ,490,490 
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463 bi kK jI=DCKD*DX 
GO TO 480 
4834 DK= BDABS(1.E-30/D(K)) 
IF( 1-DxX )4835,480,492 
WiTH THIS SCHEME I WOULD SEEM THAT D{ K) MIGHT CSCILLATE, 
BUT FURTHER EMPIRICAL STUDY iS NEEDED HERE. IT iS WORTH 
OBSERVING THAT D(K) IS RARELY AND MODESTLY INCKEASED, BUT 
( IN THIS BGX) SHARPLY- BET= 21- DECREASED. 


455 DUK jJ=D(K )*DE 
490 Mi=Miti 

EN STM 492 THE TRANSFER TO #4895" TRANSFERS CONTROL TO STM 
$905 VIA THE LONG WAY AROUND. [1 SHCULD BE MORE DIRECT. 


492 GO TO (495,493),L4 
495 GO TO 309 
493 M2=M2t1 
780 M2=1 
NPE=NPPT1 
309 GG TO 305 


BOX 13, CARDS THRU « PURPOSE- ADAPTIVE LCGIC To 
COMPUTE Gk, A MULTIPLICATIVE FACTOR THAT GOVERNS THE SiZeE 
OF THE PATTERN MOVE ATTEMPT. ENTRANCE- FROM BOX i8(FCTI) 
Via STM 510. EXiI- TO 883 IN BOX 14. COMMENT-— BY 
OBSERVING THE CUTPUT FROM PATS IN THE HMMS MODEL, iT iS 
SEEN THAT 99 &%® OF THE TiWE GR SIMPLY GROWS BY INCKEMENTS 
OF 21 FROM 2.22 TO 3.5 --THEN IS RESET TO 2.2 FOR ANOTHER 
CYCLE. THE ADAPTIVE LOGIC SEEMS TO OFFER GREAT 
GPPORTUNITIES FOR iMPROVEMENT. SOTE, AMONG CTHER 
THINGS, THAT if AN EXISTING PATTERN IS BROKEN A NEW ONE 
IS FOUND (iF PGSSIBLE) BUT GR iS NCT RESET TO 2e2- ALSO, 
SOMETHING SHCULD BE ADDED TO THE PRINTED OUTPUT OF THIS 
SUBROUTINE TO MARK WHEN AND UNDER WHAT CIRCUMSTANCES AN 
EXASTING PATTERN I5 BROKEN AND ANOTHER ONE FOUND. 


510 LT7J=LT7+1 

KOUNT=KOUNTt1 
THiS FUNCTICN ‘MCD’ IS ALWAYS EQUAL TO ZERO IN THIS CODE, 
SINCE THE SECOND ARGUMENT ALWAYS EQUALS 1. THIS IS A 
FUNCTION WHICH SHOULD BE REDESIGNED TO SERVE A PURPOSE, 
OR BETTER YET, ELIMINATED ALTOGETHER. 


iF( MOD( KOUNT,1))&888, 398, 888 
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JS8 V=( OLDSN-SN )* 100/17 OLDSN*{ NEVAL-NEVOLD )) 
iF{ 26-V)801,781,781 
731 LF( V¥—.9)502 ,782,7852 
782 iF OLDV—-V )803,804, 804 
S0i ib4=1D4ti 
GO TO #815 
802 iDi=ibditi 
Di=.1 
GO TO 810 
803 iD2=1D2t1 
GO TO 815 
804 i1D3=1D3ti 
D1I=.05 
810 GR=GRtDI 
783 GR=2.2 
815 OLDSN=SN 
OLDV=V 
NEVOLD=NEVAL 


BOX 14, CARDS THRU 2 PURPOSE-- SET X(1) TO MAKE 
PATTERN MOVE ATTEMPT.» ENTRANCE- FROM BOX 13 VIA STM 8&8 
EXLIS—- TO 270 iN BOR 18 (FCT1) ALMOST ALWAYS, TC 500 IN 
BOX 16 WHEN NEVAL VICLATES LIMiT CAUSING A RETURN TO MAIN 
ROUTINE. COMMENT—- NOTE THAT A BOUNDARY CHECK I5 CONTAINED 
IN THIS BCX AND iS SLAGHTLY DIFFERENT FROM THE BOUNDARY 
CHECK IN BOX 17 INSOFAR AS THE CONSEQUENCES CF FAILING 
THE CHECK AKE CCNCERNED. ALSG NOTE THAT SiNCE ITR=1 THIS 
MOD FUNCTIGN IN STM &8588 PERFORMS NO USEFUL PRINITCUT 
CONTROL SEKVICE. iT SHOULD BE REMOVED. 


888 GF( MOD( KOCUNT, 1918 ))7,8,7 
8 LF (KPRESS»NE-0) GO TO 7 
®WRITE( 6,991 JNEVAL, KOUNT,LI2,L1I3,LT4,LTS,LTG,LT7T gk yKkK, 
iM1,M2,NPF,LA,L4,1D1,1ED2,1D3,1D4,SN,SJ, Vy GR 
WRITE( 6,992) KV, (D(i),f=1,N) 
784 Kk=1 
IFC LSN )330,530,785 
735 SP=SN 
530 SC=SP 
LA=4 
P=Q( i) 
O(i )=K(1T) 
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RC i )=PtGR*( AC 1)-P ) 

TPC X( 1 I-XMAK(1))550,570, 786 
756 XC i jJ=XMAK( I) 

GO TO 570 
$50 IPC XMIN( 1 )-X%( 1))570,570,787 
787 &£C L)=KMINI I) 
570 CONTINUE 

GO TO 270 


BOX 15, CARDS THRU e PUKPOSE- RESET CERTAIN CCUNTEKS 
(SEE GLOSSA&SY) PRICK TO BEGINNING ADJUSTMENTS OF THE 
PATTERN MOVE ATTEMPT. ENTRANCE- FROM BOK 18 (FCT1) VIA 
STM S80. EXIT- I0 410 iN BOX 5 ALWAYS. 


580 LT4=LT4ti 
L4=2 
SP-SN 
GO TO 210 


BOX 16, CARDS THRU 2 FURPOSE- TO PRINT OUT FINAL 
PARAMETER VALUES AND TERMINATE THE SEARCH. ENTRANCE 
FROM BOXES i7 AND 18 (DUMY), PROM BOXES 8 AND 14, ALL ViA 
STM 500. EAXIT- (RETURNS TO MAIN ROUTINE AL¥AYS) 


S00 LA=8 
990 WRITE( 6, 7534 ) 


7538 FORMAT (1HO) 


1M1,M2,NPF,LA,L4,1D1i,102,1D3,1D4ySNySJ, Vy GR 

WRITE( 6,992) KV, (D(1),T=1,N) 
991 FORMAT (1H ,15, 414, 1X, 314, 413, 312, 1X, 415, 

12F14045F505,F522) 
$92 FORMAT(13,1€F 7-3) 
993 FORMAT(GX,1€F6.0) 

RETURN 
END 
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43300 
2636 
43020 
29220 
550.0 
180.0 
40020 
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2 15 
606.0 
30.0 
70020 
80020 
1200.20 
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3 17 
5000020 
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10000.0 
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10000.90 
5000.0 
200020 
450020 
70000.0 
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4 53 
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600 2b 
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207 


Bee 


4500026 


600020 
8000.0 
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1000020 
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2000.0 
4500-C 
7000026 
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1.0 


600020 
$000.0 
40.6 
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1000.0 
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230 
2J0 
2 2¥ 
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203 
20000 
03 
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219 
20G1 


220 220 
260 210 
2&3 10.0 
2420 1.20 
210 22 
2J0 210 
4205 6205 
220 a2 
207 2015 
2 002 2005 
95.0 9820 
65.20 71.20 
44720 441.0 
6202 750.20 


(END OF SAMPLE 
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DATA DECK) 


240 
$00.00 
220 
205 
280 
80200 
1.2 
&e0 


2015 


90.0 
6520 
40020 
50020 


189 


070 
210 
050 
210 
219 
60.00 
210 
210 


2015 


70.0 
6820 
4U02U 
30020 


UVe 
Uls 
Vie 
Ube 
Ol-s 
UU ald 
O9s 
Uls 


GiVe 


Veve 
U ste 
VU seOP 
U .UUe 


APPENDIX 3 


REGRESSION MODEL SUBROUTINE 


The regression modei subroutine computes the values 
of the coefficients (X, Y and Z) of the regression equation 
(WI* = X + Y (WD) 4) to be used in estimating the indirect 
work force requirements within the aggregate production 
planning model (refer to Section 3.5 for complete details). 
The computer program, written for this subroutine in 
FORTRAN IV, was successfully run on an Amdahl 470V/6 
installation at the University of Alberta Computing Services. 
The input data requirements for using this subroutine are 
as given below. 

Input Data. One requires to furnish the following 
information in the input data: 

(1) The minimum level (X) of the indirect work force 
required to maintain the plant (operations) ata 

zero direct work force level. 

C2) The historical data about the direct work force levels 

(WD) the corresponding indirect work force levels 

(WI) during the past periods/months. 

(3) The number (N) equal to the number of the data sets 

(the values of WD and WI both) included in the input 

data. 


Data Deck. All data cards follow the same general 
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format as shown below: 


Numeric 
Values 


All numeric values appearing on the data cards must 


be integer numbers and should be placed (punched), within 

the respective fields as far to the right as possible. 

Each data card is explained in detail below: 

G)) The first data card will have the number N in the 
field one. 

(2) The second data card will have the value of (X) in 
the field one. 

(3) The third and each subsequent card will have the 
values of (WI) and (WD) in the fields one and two 
respectively. 

A listing of the computer program for the regression 


model subroutine is given on the following page. 


Fy 7 
a. a 
ai 
— — 
cs s Fe 
| a S| 9) 
- 
4 
} | 
OSPF 
: : 
re rg: os 


7] mo 2 tinsogge 2 pa 
ok yougd) beastg. sd bivone bap aiedm 


— a 


ayven | iy Ww DYED? sosb bnopez edt 


_—- 


sno bier salt, 
‘od [Tiw. bted Jasupeedse toes, baa. by fy dt (e) 
ns sho z2biett afd, mt (dW) pat A to 38 
a “Misvitoeqeen " 


sigas oft 46? meteOTG Xedhanes sit 90. atte pa. 


oniwol{o? ona 0 neve ah anh 


= 


‘ 7 
a re 
och I — 


Cc 
Cc 


SUBROUTINE FOR DEVELOPING THE REGRESSION EQUATICN FOR 


FORECASTING THE INDIRECT WORKFORCE REQUIREMENTS 


IMPLICIT REAL*8 ( A-H,O0-Z) 
SP=0 

$Q=0 

SQos=0 

SPQ=0 

READ (5,51) WN 
BREAD (5,51)i1i 

R=1 

READ (5,52) J,K 
p=J-1 

Q=K 

P=DLOG(P ) 
Q=DLOG(Q) 
SP=Sp+p 

SQ=S0O+Q 
SPQ=SPQt( P*C/N) 
SQS=SCQSt{ G*%*2 DIN 
SPA=SP/N 
SQGA=SO/N 
SSQ=SQA**2 

B=( SPO-( SPA*S@QA ))/( SUS-SSO) 
A=SPA-( B¥SQA) 
A=DEXP(A) 
FORMAT( 110) 
FORMAT( 2110 ) 


FORMAT( 5Xy*X=',F1205,4 Y=4,F1205,* 


STCP 
END 


Z=4 eF12-5) 


BS) 


ae 


APPENDIX 4 


METHODE FORT ESTIMATING THE, EFFECTIVE 
LEARNING RATE CONSTANT FOR THE 
COMBINED PLANT-OPERATIONS 


The method, presented herein, for estimating the 
effective learning rate constant for use in the aggregate 
production planning model is an application of the correla- 
tion curve for the learning curve ratios developed by 
Titleman [45]. This correlation curve (shown in the 
Figure A4.1) is based on the empirical data and defines 
the relationship between the job learning rate constants 
and the manual ratios (the proportion of the manual elements 
involved in an operation) for the manufacturing operations. 
For any measured manual ratio this curve gives the value 
of the learning rate constant within a maximum error of 2.5 
percent. The effective learning rate constant for the com- 
bined plant-operations can be determined in the following 
manner: 

First list all the individual operations and deter- 
mine the proportion that each operation occupies in the 
total available production man-hours. Then multiply the 
manual ratio of each operation by its proportion of the 
total production man-hours. This gives a weighted manual 


ratio for each operation. The sum total of these weighted 
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‘manual ratios gives the effective manual ratio for the com- 
bined plant-operations. This effective manual ratio can 
now be used to read corresponding effective learning rate 


constants from the correlation curve (Figure A4.1). 
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